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METHOD OF ANALYZING ELECTROMAGNETIC INTERFERENCE 

Background of the Invention 

1. Field of" the Invention 

5 The present invention relates to a method of 

analyzing electromagnetic interference (EMI) , and more 
particularly, to a method of analyzing electromagnetic 
interference arising in a large-scale, high-speed LSI 
(large-scale integrated circuit) by means of high-speed, 
10 highly-accurate logic simulation. 

2 . Desaription of the related Art 

As is often the case, in the field of computers, LSIs 
find a broadening range of applications, from 
communications devices, such as cellular phones, to 

15 general household products, toys, and automobiles. 

Electromagnetic interference arising in such a product 
induces radio interference noise to arise in a receiver, 
such as a TV set or a radio, or faulty operations of another 
system. In order to prevent this problem, the entirety 

20 of a product is shielded, or filters are provided in a 
product. With a view towards preventing an increase in 
the number of components and cost and difficulty in 
preventing occurrence of electromagnetic interference in 
a product, strong demand exists for suppression of noise 

25 in an LSI alone. 
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Under such a situation, an LSI is ranked as a key 
device for any product which contains an LSI. Demand 
exists for an larger-scale, high-speed LSI for ensuring 
competitiveness of a product. In a situation in which the 
5 cycle of product development becomes shorter, design- 
automation of an LSI is indispensable for satisfying the 
demand. There is growing necessity for adopting 
synchronous circuit as a condition for introducing a 
state-of-the-art design-automation method. In a case 
10 where all circuits of a large-scale, high-speed LSI 
operate synchronously with a reference clock signal, 
instantaneously-changing current becomes very large, 
which induces an increase in electromagnetic 
interference . 

15 The present invention relates to a simulation method 

which enables evaluation of EMI indispensable for reducing 
electromagnetic interference while maintaining a tendency 
toward a larger-scale, higher-speed LSI. 

Noise imposed on another device by an LSI is roughly 

20 classified into two types; radiation noise, and conduction 
noise. Radiation noise emanated directly from an LSI 
includes noise emitted from internal wires of an LSI . 
However, internal wires do not act as an antenna of large 
size. As a matter of course, it is considered that the 

25 noise emitted directly from an LSI will pose a problem in 
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the future, in association with an improvement in the 
operation frequency of an LSI. However, as of now, the 
noise emitted from the inside of an LSI is considered 
trivial . 

5 In contrast, conduction noise affects another device 

mounted on a printed board by way of direct 
interconnections, such as internal wires of an LSI, lead 
frames, or wiring patterns provided on a printed wiring 
board. Noise is emitted from such interconnections while 

10 the interconnections are taken as the source of 

origination or as an antenna. The antenna constituted of 
the interconnections is much larger than that constituted 
by internal wires of an LSI and is a dominant element in 
terms of electromagnetic emission. 

15 A power line and a signal line can act as paths along 

which conduction noise developing in an LSI travels. In 
consideration of an electromagnetic field in the vicinity 
of an LSI, noise which results from variation in an 
electric current of a power source being emitted from a 

20 power line serving as an antenna is considered to be 

dominant. There may be a case where a ringing overshoot 
stemming from variation in a signal poses a problem. 
However, there more frequently arises a case where 
variation in an internal power level of an LSI propagates 

25 as a signal waveform, to thereby presents a problem. Noise 
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emitted from a power line or a signal line is considered 
to have a strong correlation with variation in the electric 
current of a power source (hereinafter referred to as a 
source current") . 
5 A source current of a CMOS circuit will now be 

described by reference to a simple inverter circuit. In 
a case where variation arises in a voltage applied to an 
inverter circuit, there flows a load capacitance 
charge/discharge current, which is the primary source 

10 current of the CMOS circuit. In addition, a short circuit 
current flows together with the load capacitance 
charge/discharge current. In design of such a CMOS 
circuit, all circuits of an LSI are synchronized in 
accordance with constraints on the use of a design- 

15 automation tool. As a result of all circuits being 
synchronized, all circuits of the LSI operate 
simultaneously, and a peak current arises in a power source 
in synchronism with a reference clock signal. Further, 
in order to increase operating speed, or shorten a cycle, 

20 of the LSI, the capacitance of a transistor is increased 
so as to enable a charging/discharging operation to be 
completed within a short period of time. Eventually, a 
peak current increases. As a matter of course, the total 
source current of an LSI is increased even when the scale 

25 of an LSI is increased. Thus , the peak current of the power 
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source is increased, thereby inducing occurrence of an 
abrupt change in a source current. Such an abrupt change 
induces an increase in higher harmonic components, thereby 
resulting in an increase in electromagnetic interference. 
5 Highly-precise simulation of change in a source 

current, which may be said to primarily account for 
electromagnetic interference, is considered to be 
effective in evaluation of electromagnetic interference 
arising in an LSI. 

10 A current simulation method for effecting 

transistor-level current analysis, as will be described 
below, has conventionally been employed. 

FIG. 116 is a block diagram showing the flow of 
processing operations pertaining to a conventional EMI 

15 analysis method. According to this method, on the basis 
of layout information pertaining to an LSI to be analyzed 
through use of a method of analyzing an electric current 
of a transistor, there is performed processing pertaining 
to a layout parameter extraction (hereinafter referred to 

20 simply as an "LPE") step 03 . Subsequently, there is 

performed processing pertaining to a circuit simulation 
step 06 for effecting simulation of a circuit on the basis 
of a switch-level netlist; a source-of -current modeling 
step 08; a source line LPE step OlO; a transient analysis 

25 simulation step 012; and a FFT step 014. 
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Processing pertaining to each of the foregoing 
processing steps will now be described by reference to FIG. 
116. 

First, there are input data Ol pertaining to the 
5 layout of a semiconductor integrated circuit (hereinafter 
referred to simply as ^^layout data Ol") , which circuit is 
to be subjected to EMI analysis; parameters of elements, 
such as transistor elements or various parasitic wiring 
elements <e.g. , resistors and capacitors) ; and an LPE rule 

10 02 for defining a form in which extracted layout parameters 
are to be output. In step 03 , parameters of the respective 
elements included in the layout data Ol are calculated in 
accordance with the LPE rule 02 , to thereby prepare a 
netlist 04 . In step 03 , parameters of parasitic elements 

15 of a power source (and the ground) are not extracted. 

In step 06 are input the netlist 04 prepared in step 
03 and a test pattern 05 for replicating a desired logic 
operation in a circuit serving as an object of analysis. 
There are calculated a load capacitance charge/discharge 

20 current and a short circuit current, both responding to 
the operating state of an internal circuit, wherewith 
information Ol about the waveform of an electric current 
of each transistor (hereinafter referred to simply as 
"current waveform information Ol") is produced. The 

25 first operation of the processing pertaining to step 06 
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is effected on the assumption that the potential of a power 
source (and that of ground) is a variation-free, ideal 
potential . 

In step 08, the current waveform information pieces 
5 07 concerning the respective transistors prepared in step 
06 are entered. Each of the thus-entered current waveform 
information 07 pieces is modeled into a mode which can be 
applied to subsequent step 012 , wherewith information 09 
concerning a source-of -current element model (hereinafter 
10 referred to simply as "source-of -current element model 
information 09") is prepared. According to a common 
method, in order to alleviate the load associated with 
processing pertaining to step 012 and subsequent steps, 
each functional circuit block consisting of a plurality 
15 of transistors is modeled as a source-of -current element. 

Processing pertaining to step OlO differs from 
processing pertaining to step 03 , only in that parameters 
of parasitic elements of a power source and parameters of 
those of parasitic elements of ground wiring (e.g., 
20 resistors, decoupling capacitance , and like elements) are 
extracted and subjected to EMI analysis in lieu of 
extracted parameters of transistor elements and extracted 
parameters of various parasitic wiring elements. Hence, 
repeated explanation of processing pertaining to step OlO 
25 is omitted. In step 010, a power source (and ground) 




wiring netlist Oil is prepared. 

In step 012 , there are entered the source-of -current 
element model information 09 prepared in step 08 , the power 
source (and ground) wiring netlist Oil prepared in step 
5 OlO, and an impedance 016 (including, resistance, 

capacitance, and inductance) of a wire or lead frame. 
Through analysis of these input data carried out by a 
transient analysis simulator typified by SPICE, variation 
in the line voltage of a circuit to be analyzed is 
10 calculated, wherewith a line voltage drop result 016 is 
prepared . 

Subsequently, processing pertaining to step 06 is 
performed again. In contrast with the first operation of 
the processing pertaining to step 06 having been effected 

15 on the assumption that the potential of the power source 
(and the ground) is a variation-free, ideal potential, the 
line voltage drop result 016 prepared in step 012 is 
entered, and the current waveform information pieces 07 
concerning the respective transistors are prepared again 

20 in consideration of the line voltage drop. Similarly, 
processing pertaining to steps 08 and 012 is performed 
again . 

Processing pertaining to steps 06, 08 , and 012 is 
effected several times in a looped manner, wherewith there 
25 is produced a current waveform result 013 which duplicates 
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variation in a line voltage with high precision. 

In step 014 , the current waveform result 013 prepared 
in step 012 is entered and subjected to fast Fourier 
transformation (hereinafter abbreviated FFT) , to thereby 
5 enable frequency spectrum analysis. There is obtained an 
EMI analysis result 015, 

In the conventional example, the precision of 
verification varies greatly according to combination of 
processing pertaining to the LPE step 03 , processing 

10 pertaining to the source line LPE step 010, and processing 
pertaining to the source-of -current modeling step 08 . 
However, a certain level of accuracy of analysis can be 
expected. A transient analysis simulator typified by 
SPICE is used for transistor-level analysis of an electric 

15 current. Hence, a limitation is imposed on the level of 
a circuit to be analyzed, and an enormous amount of 
processing time is required. The level of a semiconductor 
integrated circuit has increased recently, and 
establishment of an EMI analysis method which enables 

20 high-speed analysis of an electric current on a larger 
scale than a transistor level by dealing with the analysis 
in abstract more than a transistor level dealing is 
desired . 

A gate-level current analysis method (hereinafter 
25 referred to as a ^^gate-level power consumption analysis 
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method") has conventionally been proposed as a current 
analysis method which can be made faster . This gate-level 
current analysis method is used for analyzing power 
consumption . 

5 The gate-level power consumption analysis method for 

estimating power consumption on a gate level comprises the 
following steps: namely, a step of preparing beforehand 
a library from the total amount of electric current 
required per change when a logic signal appearing at an 

10 output port of a circuit element is changed from a logic 
value of 0 to 1 or from a logic value of 1 to 0; and a step 
of adding, to a total amount of current of the entire 
circuit at each time, a total amount of electric current 
required per change when a logic signal appearing at an 

15 output port of an instance in a circuit, the circuit being 
an object of estimation of power consumption, is changed 
from a logic value of 0 to 1 or a logic value of 1 to 0, 
to thereby estimate the total power consumed by the overall 
circuit (i.e., a value obtained by means of multiplying 

20 the total amount of current by an external voltage) . FIG. 
109 shows the configuration of an example system for 
effecting the conventional, popular gate-level power 
consumption analysis method. The gate-level power 
consumption analysis system shown in the drawing is 

25 operated through use of an amount-of -current analyzer 
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comprising circuit connection information storage means 
1001; signal pattern storage means 1002; means 1003 for 
storing the amount of electric current flowing through an 
element (hereinafter referred to simply as ^^element 
5 current storage means 1003") ; total-amount-of-current 
storage means 1004; and total-amount-of-current 
calculation means 1005, 

Individual elements constituting the gate-level 
current analysis system shown in FIG, 109 will now be 

10 described. Procedures for calculating the amount of 

current shown in FIG. 113 will be described, by reference 
to the gate-level circuit connection information shown in 
FIG. 110, a pattern of signal change arising at each 
terminal shown in FIG. Ill, and the amount of current 

15 flowing through respective elements shown in FIG. 112. 

The circuit connection information storage means 
1001 stores circuit connection information concerning 
connection of circuits which are objects of analysis . 
Circuit connection information, such as that shown in FIG. 

20 110, is stored in the circuit connection information 

storage means 1001 in advance. The circuit connection 
information is generally called a netlist. 

The circuit connection information comprises an 
external input port D, an external clock signal input port 

25 CLK, an external output port Q, flip-flops FFl and FF2 , 
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a buffer element BUFl , and wires Dl and Ql for 
interconnecting elements . Each of the flip-flops FFl and 
FF2 comprises a data input pin D, a clock signal input pin 
CK, and a data output pin Q. A logic value appearing at 
5 the data input pin D immediately before arrival of a 
leading edge of a clock signal input to the clock input 
pin CK is transmitted to the data output pin Q. 

The buffer element BUFl has a data input pin A and 
a data output pin Y. A logic value appearing at the data 

10 input pin A is output to the data output pin Y. 

The above elements have a lag, such as an intra- 
element lag, and permit propagation of a signal with a time 
lag of 1 ns . The external input port D is connected 
directly to the input pin D of the flip-flop FFl. The 

15 external clock input port CLK is a terminal for receiving 
variation in a clock signal from the outside and is 
connected directly to the input pin CK of the flip-flop 
FFl and to the input pin CK of the flip-flop FF2 . The 
external output port Q is connected directly to the output 

20 pin Q of the flip-flop FF2 and outputs a signal value of 
the flip-flop FF2 . A wire Dl interconnects the output pin 
Y of the buffer BUFl and the input pin D of the flip-flop 
FF2 . A wire Ql interconnects the output pin Q of the 
flip-flop FFl and the input pin A of the buffer BUFl . 

25 The signal pattern storage means 1002 stores a 
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pattern of signal change applied to information concerning 
connections of a circuit to be analyzed. As shown in FIG. 
Ill;, a pattern of change in a signal input to the input 
port CLK and a pattern of change in a signal input to the 
5 input port D shown in FIG. 110 are stored in the signal 
pattern storage means 1002 in advance. The horizontal 
axis in FIG. 110 represents the time at which change arises 
in a signal pattern. Upon each change, the signal pattern 
is changed between a logic value of 0 and a logic value 
10 of 1. 

The element current storage means 1003 stores 
information concerning the amount of electric current 
flowing through an element; that is, information 
concerning the total amount of electric current flowing 

15 when single change arises in an external terminal of each 
element. Information concerning the amount of electric 
current flowing through an element, such as that shown in 
FIG. 112, is stored beforehand in the element current 
storage means 1003. Reference numeral 1301 designates 

20 the amount of electric current flowing through a flip- 
flop FF when change arises in a signal input to the clock 
signal input terminal CK of the flip-flop FF . Reference 
numeral 1302 designates the amount of electric current 
flowing through the flip-flop FF when change arises in a 

25 signal output from the output terminal of the flip-flop 
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FF . Reference numeral 1303 designates the amount of 
electric current flowing through the buffer BUF when 
change arises in a signal output from the output terminal 
Y of the buffer BUF. 
5 The total -amount-of -current calculation means 1005 

is given a flowchart such as that shown in FIG. 114. 

In step 1501, the circuit connection information 
which is stored in the circuit connection information 
storage means 1001 and is shown in FIG. 110 is read. In 

10 step 1502, a signal pattern stored in the signal pattern 
storage means 1002 and shown in FIG. Ill is read. In step 
1503, information concerning the amount of electric 
current (hereinafter referred to simply as ''^element 
current information") stored in the element current 

15 storage means 1003 and shown in FIG. 112 is read. 

In step 1504, the signal pattern shown in FIG. Ill 
is delivered to the circuit connection information shown 
in FIG. 110, to thereby enable propagation of a signal. 
When change arises in a signal appearing at the external 

20 terminal of a circuit element written in the element 
current information, the amount of electric current 
flowing through the external terminal is added to the 
information concerning the total amount of electric 
current at each time. A signal propagation method used 

25 herein is identical with a signal propagation method 
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employed in a timing simulator. More specifically, the 
method is to calculate, in time sequence, change arising 
in a real circuit in accordance with a given signal pattern 
by means of tracking circuits which are affected by the 
5 signal pattern. The signal propagation method employed 
in a timing simulator is generally embodied in the form 
of a plurality of commercial tools. Hence, detailed 
description of the signal propagation method is omitted. 
FIG. 112 shows change in a logic value of a signal appearing 

10 at each pin, which has been caused as a result of 
propagation of the signal pattern in each step. 

FIG. 112 shows change in a logic value in the same 
form as that in which the signal pattern is shown in FIG. 
111. After completion of propagation of the signal, in 

15 step 1505 the information concerning the total mount of 
electric current calculated at respective times shown in 
FIG. 113 is stored in the total -amount-of -current storage 
means 1004. The total-amount-of -current storage means 
1004 stores the information concerning the total amount 

20 of electric current detected at respective times according 
to the foregoing procedures , in the form shown in FIG. 113. 

In a case where information concerning the total 
amount of electric current is calculated in greater detail , 
a load dependence coefficient of a current which is 

25 dependent on an output load capacitance (or on a charge 




current) and the total amount of electric current which 
is not dependent on an output load capacitance (or on a 
short circuit current) are stored in an element current 
storage section. Thereafter, information concerning the 
5 capacitance of each wiring (i.e., capacitance 

information) is read, the sum of the load-independent and 
the product of the capacitance information and the load 
dependence coefficient is taken as the total amount of 
electric current. As described above, according to the 

10 conventional gate-level power consumption calculation 
method shown in FIG. 110, the amount of electric current 
flowing upon a single change arising in a signal appearing 
at an input or output pin is taken as a unit or as flowing 
momentarily. In other words, the conventional method 

15 determines only the total amount of electric current. In 
terms of power consumption, the accuracy of the thus- 
calculated total amount of current is sufficient. 
However, EMI analysis requires information concerning 
chronological change in electric current, and the total 

20 amount of electric current is not sufficient in terms of 
EMI analysis. 

The conventional example using the transistor-level 
current analysis method can be expected to yield a certain 
level of accuracy. However, a transient analysis 

25 simulator typified by SPICE is used for transistor-level 
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current analysis. A limitation is imposed on the level 
of a circuit to be analyzed, and an enormous amount of 
processing time is required. The scale of a semiconductor 
integrated circuit has increased recently, and 
5 establishment of an EMI analysis method which enables 
high-speed analysis of an electric current at a larger 
level than a transistor level is desired. 

In the conventional example using the gate-level 
current analysis method, high-speed analysis of a current 

10 is feasible. However, only the total amount of electric 
current is determined. The gate-level current analysis 
method is sufficient in terms of power consumption but in 
sufficient in terms of EMI analysis. 
Summary of the Invention The present invention has 

15 been conceived to solve the drawbacks of the conventional 
methods and is aimed at providing a method for enabling 
valuation of electromagnetic interference in an LSI 
through simulation within a realistic period of time. 

To this end, an EMI analysis method according to the 

20 present invention is characterized by means of adding a 
function of highly precisely analyzing a change in power 
current to a high-speed gate-level simulation. 

The present invention provides an electromagnetic 
interference analysis method for analyzing the amount of 

25 electromagnetic interference arising in an LSI by means 




of performing a logic simulation, the method comprising: 

an instantaneous current calculation step of 
calculating the amount of instantaneous electric current 
from event information, the information being produced 
5 when a change arises in a signal and including the instance 
name of each cell in which the change has arisen, the name 
of the signal, a time at which the change has arisen, and 
transition information; 

a modeling step of modeling the instantaneous 

10 electric current according to a predetermined rule; and 
an FFT processing step of subjecting to fast Fourier 
processing (hereinafter referred to as '^^FFT 
processing" ) the information concerning a change in 
electric current, the information being calculated 

15 through a modeling step. Here, the amount of 

instantaneous current means the amount of electric current 
induced by an event. The amount of electric current 
flowing through a cell when a change arises in an output 
of the cell as an event is taken as the amount of 

20 instantaneous electric current. 

Preferably, the modeling step includes an averaging 
step of averaging the instantaneous current over a certain 
discrete width, and the FFT processing step includes a step 
of subjecting to FFT processing information concerning a 

25 change in current, the information being produced by the 
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averaging step. 

Preferably, the modeling step includes a rectangular 
waveform modeling step of modeling the instantaneous 
current as a rectangular waveform whose height is 
5 calculated for information for each event such that the 
area of the rectangular waveform becomes equal to the 
amount of electric current of each event, and the FFT 
processing step includes a step of subjecting to FFT 
processing information concerning a change in current, the 

10 information being calculated in the rectangular waveform 
modeling step. 

Preferably, the modeling step includes a 
geometrically-similar rectangular waveform modeling step 
of modeling the instantaneous current as a 

15 geometrically-similar rectangular waveform whose height 
and width are calculated such that the area of the 
rectangular waveform becomes equal to the amount of 
electric current of each event, and the FFT processing step 
includes a step of subjecting to FFT processing 

20 information concerning a change in current, the 

information being calculated in the geometrically-similar 
rectangular waveform modeling step. 

Preferably, the modeling step includes a rectangular 
waveform modeling step of calculating the instantaneous 

25 electric current for information for each event, and a step 




of modeling the instantaneous current as a rectangular 
waveform through use of the amount of electric current and 
a table representing the relationship between the width 
and height of a rectangular waveform^ to thereby subject 
5 to FFT processing the information concerning a change in 
electric current calculated in the rectangular waveform 
modeling step. 

Preferably, the modeling step includes a step of 
calculating the instantaneous electric current for 

10 information for each event, and a rectangular waveform 
modeling step of modeling the instantaneous current as a 
rectangular waveform through use of a slew in input 
waveform and a table representing the relationship between 
the width and height of a rectangular waveform, to thereby 

15 subject to FFT processing the information concerning a 
change in electric current calculated in the rectangular 
waveform modeling step. 

Preferably, the modeling step includes a step of 
calculating the instantaneous electric current for 

20 information for each event, and a rectangular waveform 
modeling step of modeling the instantaneous current as a 
rectangular waveform through use of an output load 
capacitance and a table representing the relationship 
between the width and height of a rectangular waveform, 

25 to thereby subject to FFT processing the information 




concerning a change in electric current calculated in the 
rectangular waveform modeling step. 

Preferably, the modeling step includes a step of 
calculating a drop in voltage from the amount of electric 
5 current flowing in each cell and the resistance of a power 
and correcting the amount of instantaneous electric 
current of each cell for each event, on the basis of the 
relationship between the drop in voltage and the amount 
of instantaneous electric current. 

10 Preferably, the modeling step includes a step of 

calculating a drop in voltage from the amount of electric 
current flowing in each cell, the resistance of a power 
line, and the capacitance of an on-chip capacitor, and 
correcting the amount of instantaneous electric current 

15 of each cell for each event, on the basis of the 

relationship between the drop in voltage and the amount 
of instantaneous electric current. 

Preferably, the modeling step includes a step of 
transiently-analyzing a power RC of each cell and a cell 

20 power source, accurately calculating a drop in voltage, 
and a correction step of correcting the amount of 
instantaneous electric current of each cell for each event, 
on the basis of the relationship between the drop in 
voltage and the amount of instantaneous electric current. 

25 Preferably, the modeling step includes a triangular 
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waveform modeling step of modeling the instantaneous 
current as a triangular waveform which has a given width 
and whose height is calculated for information for each 
event such that the amount of instantaneous electric 
5 current becomes equal to the area of the triangular 

waveform, and the FFT processing step includes a step of 
subjecting to FFT processing information concerning a 
change in current, the information being calculated in the 
triangular waveform modeling step. 

10 Preferably, the modeling step includes a multi- 

order-function waveform modeling step of modeling the 
instantaneous current as a multi -order-function waveform, 
and the FFT processing step includes a step of subjecting 
to FFT processing information concerning a change in 

15 current, the information being calculated in the 
multi-order-function waveform modeling step. 

Preferably, the modeling step includes a exponential 
function waveform modeling step of modeling the 
instantaneous current as an exponential-function waveform, 

20 and the FFT processing step includes a step of subjecting 
to FFT processing information concerning a change in 
current, the information being calculated in the 
exponential-function waveform modeling step. 

Preferably, the modeling step includes a step of 

25 modeling the amount of instantaneous electric current 
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while separating the same into a short circuit electric 
current component and a charge current component. 

Preferably, the modeling step includes a calculation 
step of calculating the height of a rectangular waveform 
5 from a library in which peak currents of cells are 
characterized according to the type of cell, and a 
rectangular waveform modeling step of modeling the amount 
of instantaneous electric current as a rectangular 
waveform, _the peak current calculated in the calculation 

10 step being taken as the height of the rectangular waveform 
and the area of the triangular waveform being equal to the 
amount of electric current of each event, and the FFT 
processing step includes a step of subjecting to FFT 
processing information concerning a change in current, the 

15 information being calculated in the rectangular waveform 
model i ng s t ep . 

Preferably, the calculation step includes a step of 
calculating a peak current for each cell from information 
concerning a slew in the cell, by reference to a library 

20 in which the relationship between a slew in input waveform 
and a peak current is characterized in the form of a table 
according to the type of cell. 

Preferably, the calculation step includes a step of 
calculating a peak current for each cell from information 

25 concerning a load capacitance of a cell, by reference to 
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a library in which the relationship between a load 
capacitance and a peak current is characterized in the form 
of a table according to the type of cell. 

Preferably, the calculation step includes a step of 
5 setting a plurality of peak currents for a composite cell 
and calculating the heights of a plurality of rectangular 
waveforms through use of a characterized library, and the 
rectangular waveform modeling step corresponds to a step 
of modeling the amount of electric current into a plurality 
10 of rectangular waveforms. 

Preferably, the calculation step includes a step of 
setting a plurality of peak currents for each of the rise 
and fall of a flip-flop (FF) cell and calculating the 
heights of a plurality of rectangular waveforms through 
15 use of a characterized library, and the rectangular 

waveform modeling step corresponds to a step of modeling 
the amount of electric current into a plurality of 
rectangular waveforms. 

Preferably, the calculation step includes a step of 
20 calculating the height of a rectangular waveform through 
use of a library in which peak currents are characterized, 
in consideration of the state of an input signal. 

Preferably, the modeling step includes 

a step of calculating a drop in voltage from the 
25 amount of electric current determined according to the 
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type of cell and from the resistance of a power line; and 
a correction step of characterizing, for each cell, 
the relationship between a drop in voltage and the amount 
of instantaneous electric current in the form of a table, 
5 to thereby correct the amount of instantaneous electric 
current for each event of the cell. 

Preferably, the modeling step includes 
a step of calculating a drop in voltage from the 
amount of electric current determined according to the 
10 type of cell, the resistance of a power line, and the 
capacitance of an on-chip capacitor; and 

a correction step of characterizing, for each cell, 
the relationship between a drop in voltage and the amount 
of instantaneous electric current in the form of a table, 
15 to thereby correct the amount of instantaneous electric 
current for each event of the cell . 

Preferably, the correction step includes a step of 
iterating several times calculation of a drop in voltage 
and correction of a current waveform. 
20 Preferably, the electromagnetic interference 

analysis method as defined in any one of claims 15 through 
23, wherein the calculation step includes a step of 
modeling the amount of instantaneous electric current 
while separating the same into a short circuit electric 
25 current component and a charge current component. 
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Preferably, the modeling step includes 
a triangular waveform modeling step of modeling the 
instantaneous current as a triangular waveform whose width 
is calculated for each event information in consideration 
5 of slew information (i.e., an output slew) for an output 
terminal of a cell for each event information such that 
the area of the triangular waveform becomes equal to the 
amount of electric current of each event, the height of 
the triangular waveform being calculated on the basis of 
10 the width, and 

the FFT processing step includes a step of subjecting 
to FFT processing information concerning a change in 
current, the information being calculated in the 
triangular waveform modeling step. 
15 Preferably / the modeling step includes 

a triangular height calculation step of calculating 
the height of a triangular waveform such that the area of 
the triangular waveform becomes equal to the amount of 
electric current of each event, by means of multiplying 
20 the amount of instantaneous electric current by a 

coefficient corresponding to the state of an event of a 
cell, in consideration of whether the event of the cell 
is in a rising state or a falling state. 

Preferably, the modeling step includes 
25 a step of calculating the amount of instantaneous 
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electric current for each event information in the case 
of a composite cell; and 

a triangular waveform modeling step of modeling the 
amount of instantaneous electric current as a plurality 
5 of triangular waveforms which are equal in number to the 
stages provided in the composite cell, through use of a 
table representing the relationship between the width and 
height of a triangular waveform; and 

the FFT processing step includes a step of subjecting 
10 to FFT processing information concerning a change in 
current, the information being calculated in the 
triangular waveform modeling step. 

Further, the present invention provides an 
electromagnetic interference analysis system for 
15 analyzing the amount of electromagnetic interference 

arising in an LSI by means of performing a logic simulation , 
the system comprising: 

a logic simulator; 

computation means which is connected to the logic 
20 simulator and calculates the amount of instantaneous 

electric current from event information, the information 
being produced when a change arises in a signal and 
including the instance name of each cell in which the 
change has arisen, the name of the signal, a time at which 
25 the change has arisen, and transition information; 
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modeling means for modeling the instantaneous 
electric current according to a predetermined rule; and 
fast Fourier (FFT) conversion means for subjecting 
to fast Fourier processing the information concerning a 
5 change in electric current, the information being 

calculated through a modeling step , thereby analyzing the 
amount of electromagnetic interference arising in an LSI 
on the basis of a signal output from the FFT conversion 
means . 

10 

Brief Description of the Drawings 
FIG. 1 is a block diagram showing the overall flow 
of processing pertaining to an EMI analysis method 
according to a first embodiment of the present invention; 
15 FIG. 2 is a block diagram showing a current waveform 

calculation section according to the first embodiment; 

FIG. 3 is a graph showing a calculation result of 
an instantaneous current; 

FIG. 4 is a flowchart showing processing for 
20 averaging a current waveform; 

FIG. 5 is a graph showing a calculation result of 
a current waveform (according to the first embodiment) ; 

FIG. 6 is a graph showing a result of an FFT 
operation ; 

25 FIG. 7 shows a rectangular waveform model (according 
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to a second embodiment of the present invention) ; 

FIG. 8 is a block diagram showing calculation of a 
current waveform (according to the second embodiment) ; 

FIG. 9 is a flowchart showing a first rectangular 
5 waveform- shaping operation; 

FIG. 10 is a graph showing a calculation result of 
a current waveform (according to the second embodiment) ; 

FIG. 11 shows a rectangular waveform model 
(according to a third embodiment of the present 
10 invention) ; 

FIG. 12 is a block diagram showing calculation of 
a current waveform (according to the third embodiment) ; 

FIG. 13 is a flowchart showing a second rectangular 
waveform- shaping operation; 
15 FIG. 14 is a graph showing a calculation result of 

a current waveform (according to the third embodiment) ; 

FIG. 15 shows a rectangular waveform model 
(according to fourth, fifth, and sixth embodiments of the 
present invention) ; 
20 FIG. 16 shows an i-h table; 

FIG. 17 is a block diagram showing calculation of 
a current waveform (according to a fourth embodiment of 
the present invention) ; 

FIG. 18 is a flowchart showing a third rectangular 
25 waveform- shaping operation; 



29 




FIG. 19 is a graph showing a calculation result of 
a current waveform (according to the fourth embodiment) ; 
FIG. 20 shows an i-s-h table; 

FIG. 21 is a block diagram showing calculation of 
5 a current waveform (according to a fifth embodiment of the 
present invention) ; 

FIG. 22 is a flowchart showing a fourth rectangular 
waveform- shaping operation; 

FIG. 23 is a graph showing a calculation result of 
10 a current waveform (according to the fifth embodiment) ; 
FIG. 24 shows an i-c-h table; 

FIG. 25 is a block diagram showing calculation of 
a current waveform (according to a sixth embodiment of the 
present invention) ; 
15 FIG. 26 is a flowchart showing a fifth rectangular 

waveform-shaping operation; 

FIG. 27 is a graph showing a calculation result of 
a current waveform (according to the sixth embodiment) ; 

FIG. 28 is a flowchart for describing an EMI analysis 
20 method according to a seventh embodiment of the present 
invention ; 

FIG. 29 is a flowchart for describing the flow of 
processing operations of a line voltage drop 
calculation/instantaneous current correction section 
25 according to the seventh embodiment; 
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FIG. 30 is a schematic representation showing layout 
data which are entered into the line voltage drop 
calculation/instantaneous current correction section 
according to the seventh embodiment; 
5 FIG. 31 is a schematic representation showing 

information about segments in a power line, which 
information is produced in the line voltage drop 
calculation/instantaneous current correction section 
according to the seventh embodiment; 

10 FIG. 32 is a schematic representation showing 

information about the resistance of each of segments, 
which information is produced in the line voltage drop 
calculation/instantaneous current correction section 
according to the seventh embodiment; 

15 FIG. 33 is a schematic representation showing 

information about an electric current flowing in each of 
segments, which information is produced in the line 
voltage drop calculation/instantaneous current 
correction section according to the seventh embodiment; 

20 FIG. 34 is a schematic representation showing 

information about a drop in the voltage of each of segments , 
which information is produced in the line voltage drop 
calculation/instantaneous current correction section 
according to the seventh embodiment; 

25 FIG. 35 is a schematic representation showing 
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information about a drop in the voltage of each instance, 
which information is produced in the line voltage drop 
calculation/instantaneous current correction section 
according to the seventh embodiment; 
5 FIG. 36 is a flowchart for describing an EMI analysis 

method according to an eighth embodiment of the present 
invention ; 

FIG. 37 is a flowchart for describing the flow of 
processing operations of a line voltage drop 
10 calculation/instantaneous current correction section 
according to the eighth embodiment; 

FIG. 38 is a schematic representation showing 
information about a capacitor provided between segments, 
which information is produced in the line voltage drop 
15 calculation/instantaneous current correction section 
according to the eighth embodiment; 

FIG. 39 is a schematic representation showing 
information about a drop in voltage of each of segments, 
which information is produced in the line voltage drop 
20 calculation/instantaneous current correction section 
according to the eighth embodiment; 

FIG. 40 is a flowchart for describing an EMI analysis 
method according to a ninth embodiment of the present 
invention ; 

25 FIG. 41 a flowchart for describing the flow of 
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processing operations of a line voltage drop 
calculation/instantaneous current correction section 
according to a ninth embodiment of the present invention; 

FIG. 42 shows a triangular waveform model (according 
5 to a tenth embodiment of the present invention) ; 

FIG. 43 is a block diagram showing a current waveform 
calculation section (according to a tenth embodiment of 
the present invention) ; 

FIG. 44 is a flowchart of a triangular waveform- 
10 shaping operation (according to the tenth embodiment) ; 

FIG. 45 is a graph showing a calculation result of 
a current waveform (according to the tenth embodiment) / 

FIG. 4 6 is a multi-order-function model (according 
to an eleventh embodiment of the present invention) ; 
15 FIG. 47 is a block diagram showing a current waveform 

calculation section (according to a tenth embodiment of 
the present invention) ; 

FIG. 48 is a flowchart of a multi -order-function 
waveform-shaping operation; 
20 FIG. 4 9 is a graph showing a calculation result of 

a current waveform (according to an eleventh embodiment 
of the present invention) ; 

FIG. 50 is a schematic representation showing a 
Gaussian model (according to twelfth and twenty-fifth 
25 embodiments of the present invention) ; 
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FIG. 51 is an i-h table; 
FIG. 52 is an i-w table; 

FIG. 53 is a block diagram showing a current waveform 
calculation section (according to a twelfth embodiment of 
5 the present invention) ; 

FIG. 54 is a flowchart of a Gaussian function 
waveform-shaping operation; 

FIG. 55 is a graph showing a calculation result of 
a current waveform (according to a twelfth embodiment of 
10 the present invention) ; 

FIG. 56 is a block diagram showing the configuration 
of a system for effecting an EMI analysis method according 
to a thirteenth embodiment of the present invention; 

FIG. 57 is an illustration showing example data 
15 stored in means for storing basic information to be used 
for estimating an electric current according to the 
thirteenth embodiment; 

FIG. 58 is a graph showing example data stored in 
total current storage means according to the thirteenth 
20 embodiment; 

FIG. 59 is a graph showing example data stored in 
element current storage means according to the thirteenth 
embodiment ; 

FIG. 60 is a flowchart showing processing operations 
25 of the total current calculation means according to the 
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thirteenth embodiment; 

FIG. 61 is a block diagram showing a current waveform 
calculation section (according to a fourteenth embodiment 
of the present invention) / 
5 FIG. 62 is a flowchart of a sixth rectangular 

waveform- shaping operation; 

FIG. 63 shows an s-h table; 

FIG. 64 is a block diagram showing a current waveform 
calculation section (according to a fifteenth embodiment 
10 of the present invention) ; 

FIG. 65 is a flowchart of a seventh rectangular 
waveform-shaping operation; 

FIG. 66 shows a c-h table; 

FIG. 67 is a block diagram showing a current waveform 
15 calculation section (according to a sixteenth embodiment 
of the present invention) ; 

FIG. 68 is a flowchart of an eighth rectangular 
waveform-shaping operation; 

FIG. 69 is a graph showing example data stored in 
20 a pattern-of -signal -change storage means according to the 
seventeenthembodimentof the present invention; 

FIG. 7 0 is a graph showing an example current 
waveform arising in a real circuit according to the 
seventeenth embodiment; 
25 FIG. 71 is a graph showing example data stored in 
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conventional total current storage means according to the 
thirteenth embodiment; 

FIG. 72 is a block diagram showing the configuration 
of a system for effecting a gate-level EMI analysis method 
5 according to the seventeenth embodiment and an eighteenth 
embodiment of the present invention; 

FIG. 73 is an illustration showing example data 
stored in circuit connection information storage means 
according to the seventeenth embodiment; 
10 FIG. 74 is a graph showing example data stored in 

total current storage means according to the seventeenth 
embodiment ; 

FIG. 75 is a graph showing example data stored in 
element current storage means according to the seventeenth 
15 embodiment; 

FIG. 76 is a flowchart showing processing operations 
of the total current calculation means according to the 
seventeenth embodiment; 

FIG. 77 is an illustration showing example data 
20 stored in circuit connection information storage means 
according to the eighteenth embodiment; 

FIG. 78 is a graph showing example data stored in 
a pattern-of -signal -change storage means according to the 
eighteenth embodiment; 
25 FIG. 7 9 is a graph showing an example current 
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waveform arising in a real circuit according to the 
eighteenth embodiment; 

FIG. 80 is a graph showing example data stored in 
conventional total current storage means according to the 
5 eighteenth embodiment; 

FIG. 81 is a graph showing example data stored in 
element current storage means according to the eighteenth 
embodiment; 

FIG. 82 is a graph showing example data stored in 
10 total current storage means according to the eighteenth 
embodiment ; 

FIG. 83 is a flowchart showing processing operations 
of the total current calculation means according to the 
eighteenth embodiment; 
15 FIG. 84 is an illustration showing example data 

stored in circuit connection information storage means 
according to a nineteenth embodiment of the present 
invention ; 

FIG. 85 is a graph showing example data stored in 
20 a pattern-of -signal -change storage means according to the 
nineteenth embodiment; 

FIG. 86 is a graph showing an example current 
waveform arising in a real circuit according to the 
nineteenth embodiment; 
25 FIG. 87 is a graph showing example data stored in 
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conventional total current storage means according to the 
nineteenth embodiment; 

FIG. 88 is a graph showing example data stored in 
element current storage means according to the nineteenth 
5 embodiment; 

FIG . 89 is a flowchart showing processing operations 
of the total current calculation means according to the 
nineteenth embodiment; 

FIG. 90 is a graph showing example data stored in 
10 total current storage means according to the nineteenth 
embodiment ; 

FIG. 91 is a flowchart for describing an EMI analysis 
method according to a twentieth embodiment of the present 
invention; 

15 FIG. 92 is a flowchart for describing the flow of 

processing operations of a line voltage drop 
calculation/instantaneous current correction section 
according to the twentieth embodiment; 

FIG. 93 is a flowchart for describing an EMI analysis 

20 method according to a twenty-first embodiment of the 
present invention; 

FIG. 94 is a flowchart for describing the flow of 
processing operations of a line voltage drop 
calculation/instantaneous current correction section 

25 according to the twenty-first embodiment; 
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FIG. 95 is a flowchart for describing an EMI analysis 
method according to a twenty- second embodiment of the 
present invention; 

FIG. 96 shows a triangular waveform model (according 
5 to a twenty-third embodiment of the present invention) ; 

FIG. 97 is a block diagram showing a current waveform 
calculation section (according to the twenty- third) ; 

FIG. 98 is a flowchart of a triangular waveform- 
shaping operation; 
10 FIG. 99 is a multi-order-function model (according 

to a twenty-fourth embodiment of the present invention) ; 

FIG. 100 is a block diagram showing a current 
waveform calculation section (according to the twenty- 
fourth embodiment) ; 
15 FIG. 101 is a flowchart of a mul ti -order-function 

waveform-shaping operation; 

FIG. 102 is a block diagram showing a current 
waveform calculation section (according to a twenty-fifth 
embodiment of the present invention) ; 
20 FIG. 103 is a flowchart of a Gaussian function 

waveform- shaping operation; 

FIG. 104 is a graph showing an example current 
waveform arising in a real circuit according to a 
twenty-sixth embodiment of the present invention; 
25 FIG. 105 is a graph showing example data stored in 
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element current storage means according to the twenty- 
sixth embodiment; 

FIG. 106 is a graph showing example data stored in 
conventional total current storage means according to the 
5 twenty-sixth embodiment; 

FIG. 107 is a flowchart showing processing 
operations of the total current calculation means 
according to the twenty- sixth embodiment; 

FIG. 108 is a diagram showing example data concerning 
10 internal connections stored in the circuit connection 
information storage means according to the seventeenth 
embodiment shown in FIG. 73; 

FIG. 109 is a block diagram showing the configuration 
of a system for effecting an EMI analysis method employed 
15 in an example conventional EMI analysis method; 

FIG. 110 is an illustration showing example data 
stored in circuit connection information storage means 
employed in an example conventional EMI analysis method 
and in a twenty-sixth embodiment; 
20 FIG. Ill is a graph showing example data stored in 

a pattern-of-signal-change storage means employed in an 
example conventional EMI analysis method and in the 
twenty-sixth embodiment; 

FIG. 112 is a graph showing example data stored in 
25 element current storage means employed in an example 
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conventional EMI analysis method; 

FIG. 113 is a graph showing example data stored in 
total current storage means employed in an example 
conventional EMI analysis method; 
5 FIG, 114 is a flowchart showing processing 

operations of the total current calculation means employed 
in an example conventional EMI analysis method; 

FIG. 115 is a diagram of changes in signals flowing 
through an internal net of a circuit, the changes being 
10 detected according to an example conventional EMI analysis 
method; 

FIG. 116 is a flowchart for describing an example 
conventional EMI analysis method; 

FIG. 117 is a triangular waveform model according 
15 to the present invention (a twenty-seventh embodiment) ; 

FIG. 118 is a block diagram showing calculation of 
a current waveform (according to the twenty-seventh 
embodiment) ; 

FIG. 119 is a flowchart showing shaping of a 
20 triangular waveform (according to the twenty- seventh 
embodiment) ; 

FIG. 120 shows a triangular waveform model according 
to the present invention (a twenty-eighth embodiment) ; 

FIG. 121 is a block diagram showing calculation of 
25 a current waveform (according to the twenty-eighth 
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embodiment) ; 

FIG. 122 shows a triangular waveform model according 
to the present invention (a twenty-ninth embodiment) ; and 

FIG. 123 is a block diagram showing calculation of 
5 a current waveform (according to the twenty-ninth 
embodiment) . 

Description of the preferred Embodiments (First 
Embodiment) 

10 An electromagnetic interference (EMI) analysis 

method according to a first embodiment of the present 
invention will be described hereinbelow. The EMI 
analysis method is characterized by calculating a current 
waveform. 

15 The method comprises the steps of: 

extracting event information on a per-cell basis 
from a logic simulator; 

calculating from the event information a 
momentarily-changing current of each cell; 
20 averaging the instantaneously-changing current by 

means of dividing by a certain duration of time, to thereby 
calculate a more realistic change in current; and 

subjecting the calculation result of a current 
waveform to fast Fourier transformation (FFT) , to thereby 
25 determine a frequency characteristic of electromagnetic 
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interference (EMI) emitted from a circuit of interest, to 
thereby effect EMI analysis. 

FIG. 1 is a block diagram showing the overall flow 
of processing pertaining to the EMI analysis method of the 
5 first embodiment. In the drawing, there are described a 
netlist 1 representing data pertaining to a circuit to be 
subjected to EMI analysis; a test vector 2 for activating 
the netlist 1; a library 3 having described therein 
information concerning individual cells constituting the 

10 netlist 1; a logic simulator 4 for performing a simulation 
on the basis of the information from the library 3; and 
a current waveform calculation section 7 for calculating 
a current waveform from event information 5, the 
information having been calculated on a per-cell basis by 

15 means of the logic simulator 4 . An FFT processing section 
9 performs an FFT processing operation on the basis of a 
current waveform calculation result 8 calculated by the 
current waveform calculation section 7, to thereby 
calculate an EMI component. Accordingly, there is 

20 obtained an EMI analysis result 10 of the circuit of 

interest. Thus, the EMI analysis method is characterized 
by use of the current waveform calculation section 7 . The 
logic of a cell, delay information concerning the cell, 
and momentary-current-value information concerning the 

25 cell are described in the library 3 according to the types 
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of individual cells . Every time an event arises , the logic 
simulator 4 outputs the event information 5. The term 
^^evenf used herein designates occurrence of any change 
in a signal within the logic simulator 4 . The event 
5 information 5 contains the instance name of an object of 
EMI analysis which accounts for outputting of the event 
information 5, the name of a signal in which change has 
arisen, the time at which the change has arisen, and 
information concerning transition of the signal. The 

10 expression "instance name" corresponds to a proper name 
assigned to each of the cells constituting the netlist 1. 
Reference numeral 6 designates a simulation result output 
from the logic simulator 4. 

As can be seen from a detail block diagram shown in 

15 FIG. 2, the current waveform calculation section 7 
comprises an instantaneously-changing current 
calculation section 11 and calculates an 
instantaneously-changing current from information 
concerning an instantaneously-changing current of an 

20 instance calculated by the instantaneously-changing 

current calculation section 11, the instance serving as 
an object of EMI analysis. On the basis of a calculation 
result 12 of an instantaneously-changing current 
(hereinafter referred to as an "instantaneously-changing 

25 current calculation result 12") , the current waveform 
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calculation section 7 averages a current waveform by means 
of dividing an instantaneously-changing current by a 
predetermined duration of time. The instantaneously- 
changing current calculation section 11 calculates a 
5 momentary change in the current of each cell from the event 
information 5 concerning each cell, the information being 
output from the logic simulator 4. Calculation of 
momentary change in current is technically analogous to 
calculation of power consumption through a conventional 
10 logic simulation. Upon receipt of the event information 
5 output from the logic simulator 4, the 

instantaneously-changing current calculation section 11 
extracts, from the library 3, information concerning an 
instantaneously-changing current of an instance, the 

15 instance being an object of EMI analysis, to thereby 
calculate an instantaneously-changing current. 
Reference numeral 12 designates an instantaneously- 
changing current calculation result determined by the 
instantaneously-changing current calculation section 11. 

20 FIG. 3 is a plot showing an example instantaneously- 
changing current calculation result 12. In the drawing, 
the vertical axis represents the amount of current ^^it" 
(As), and the horizontal axis represents time "t." 
Reference symbol T designates a duration of time for 

25 averaging a current waveform by means of dividing an 
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instantaneously-changing current by a certain duration of 
time and corresponds to a predetermined constant. At 
designates a minimum time unit to be used for analysis and 
differs from the minimum unit time of the logic simulator 
5 4 . Reference symbol ''^t" denotes a current processing time 
and corresponds to a simulation time of the logic simulator 
4 . Reference symbol ^^it" denotes the value of an 
instantaneously-changing current. As can be seen from 
FIG. 3, according to the instantaneously-changing current 

10 calculation result, the waveform of an electric current 
assumes a 6 function. Therefore, an EMI component cannot 
be analyzed through an FFT operation. In order to enable 
analysis of an EMI component through FFT, a current 
waveform averaging section 13 averages the 

15 instantaneously-changing current calculation result 12 by 
means of dividing by a certain duration of time. 

FIG. 4 shows a flowchart of steps relating to 
averaging of a current waveform. The following 
processing operations are repeated until an averaging 

20 operation is completed; namely, until ''^t" exceeds an end 
time of simulation (step 1010) . First, a current buffer 
"lb" for averaging purpose is set to 0 (step 1011) . Next, 
i(t) is added to "lb" until the current time "t" becomes 
"t+T" (steps 1012 and 1013). Here, i (t) denotes an 

25 instantaneously-changing current at time "t" (step 1014) . 
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Next, lave is assigned to I (x) until a variable ^^x" shifts 
from -T to ^^t," thereby adding At to the variable "x" (steps 
1015 and 1016) . Here, I (x) denotes a current at time ^^x." 
The current waveform calculation result 8 can be 
5 determined by means of the flowchart of averaging 

operations. FIG. 5 shows a current waveform calculation 
result 8 obtained at time T=5 provided in the 
instantaneously-changing current calculation result 8 
shown in FIG. 3. The FFT processing section 9 subjects 

10 the thus-calculated waveform to an FFT processing 
operation, thereby determining the frequency 
characteristic of an EMI component. In this way, there 
can be obtained a result 10 of EMI analysis (hereinafter 
referred to as an "EMI analysis result 10") . FIG. 6 shows 

15 a sample of the EMI analysis result 10. From the EMI 

analysis result 10, the frequency characteristic of an EMI 
component arising in the circuit of interest can be 
determined, and the designer of the circuit can study a 
method of addressing EMI. 

20 The EMI analysis method of the present embodiment 

models a current waveform by means of an averaging 
operation and subjecting the thus-modeled waveform to an 
FFT operation, to thereby enable highly-precise analysis 
of EMI. Further, the EMI analysis method of the present 

25 embodiment can analyze an EMI component of a circuit within 
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a shorter period of time than can a transistor-level EMI 
analysis method. 
(Second Embodiment) 

An electromagnetic interference (EMI) analysis 
5 method according to a second embodiment of the present 
invention will be described hereinbelow . The EMI 
analysis method comprises the steps of: 

calculating a momentarily-changing current of each 
cell from the event information 5 output from the logic 
10 simulator 4 ; 

modeling the thus-calculated momentarily-changing 
current a rectangular waveform having a certain duration 
of time, the peak value of the waveform having been 
calculated such that an instantaneously-changing current 
15 coincides with the area of the rectangular waveform. 
Accordingly, a more realistic change in current is 
calculated, and the calculation result of a current 
waveform is subjected to FFT operation, wherewith a 
frequency characteristic of an EMI component arising in 
20 a circuit of interest is determined and EMI analysis is 
effected. 

FIG. 7 shows an instantaneously-changing current 
value and a modeled rectangular waveform according to the 
second embodiment. In the drawing, ^^i" represents an 
25 instantaneously-changing current value, and "t" denotes 
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a current processing time. Reference symbol "W" denotes 
the width of a modeled rectangular waveform. The width 
of a rectangular waveforms corresponds to a predetermined 
constant. Further, reference symbol "^^h" denotes the 
5 height of a rectangular waveform, which is to be calculated 
from ^^i" and "W." 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 

10 embodiment (see FIG. 1) . 

FIG. 8 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 

15 calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 

20 conventional logic simulation. Reference numeral 12 

designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
^^instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 

25 calculation section 11. FIG. 3 shows an example 
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instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
5 Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with that of the current wave calculation 
operation of the first embodiment. The EMI analysis 
method according to the present embodiment comprises a 

10 first rectangular waveform height calculation section 20 
for calculating the height of a rectangular waveform such 
that the area of a rectangular waveform serving as a model 
of the waveform of the instantaneously-changing current 
coincides with the instantaneously-changing current 

15 value; and a first rectangular waveform-shaping section 
21 for modeling the instantaneously-changing current 
calculation result 12 as a rectangular waveform having a 
certain duration of time, thereby calculating the waveform 
of an electric current. 

20 FIG. 9 shows a flow of rectangular waveform- shaping 

operations including an operation for calculating the 
height of a rectangular waveform. The following 
operation is repeated on a per-event basis until a current 
wave calculation operation is completed (step 1020) . The 

25 area of a rectangular waveform is defined as ^^W x h" and 
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assumes a value of ^^i." The height ^^h" of the rectangular 
waveform is calculated as i/W (step 1021), Symbol "i" 
denotes an instantaneously-changing current value 
concerning a cell of an event, the event being an object 
5 of EMI analysis. This processing corresponds to 

processing pertaining to the first rectangular waveform 
height calculation section 20. Next, ^^h" is added to 
I (t+x) until a variable ^^x" changes from -H/l to W/2 , and 
At is added to "x" (steps 1022 and 1023). I (t+x) 

10 corresponds to the value of electric current flowing 
through all the cells at time t+x. This processing 
corresponds to processing pertaining to the first 
rectangular waveform- shaping section 21. 

The current wave calculation result 8 can be obtained 

15 by means of the flow of a rectangular waveform height 
calculation operation and a rectangular waveform-shaping 
operation. FIG. 10 shows an example current wave 
calculation result 8 at W=5 in the instantaneously- 
changing current calculation result shown in FIG. 3. The 

20 FFT processing section 9 subjects the waveform to an FFT 
operation, to thereby determine the frequency 
characteristic of an EMI component. As a result, there 
can be obtained an EMI analysis result 10. 

According to the EMI analysis method of the present 

25 embodiment, the waveform of a current is modeled as a 
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rectangular waveform, and the thus -modeled rectangular 
waveform is subjected to an FFT operation, thereby 
enabling more accurate analysis of EMI. The EMI analysis 
method according to the present embodiment enables 
5 processing within a more shorter period of time than does 
the transistor-level EMI analysis method. 
(Third Embodiment) 

An EMI analysis method according to a third 
embodiment of the present invention will now be described. 

10 Under the EMI analysis method according to the present 
embodiment, a momentarily-changing current of each cell 
is calculated from event information concerning to each 
cell , the information being output from a logic simulator . 
The momentarily-changing current is modeled as a 

15 rectangular waveform, wherein a ratio of a duration of time 
to a peak value is made constant such that an 
instantaneously-changing current value becomes equal to 
the area of the rectangular waveform. As a result, there 
is calculated a more realistic change in current, and the 

20 result of calculation of a current waveform is subjected 
to an FFT operation, to thereby determine the frequency 
characteristic of an EMI component of a circuit of interest 
and enable EMI analysis. The EMI analysis method will now 
be described. 

25 FIG. 11 shows a momentary current value and a modeled 
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rectangular waveform according to the third embodiment. 
Reference symbol '^i" denotes an instantaneously-changing 
current value, and "t" denotes a current processing time. 
Reference symbol '^^K" denotes a ratio of a peak value of 
5 a rectangular waveform to a duration of time. In the 
present embodiment, the ratio is expressed as '^^peak- 
value/width . " The peak- value/width ratio is constant, 
and an appropriate value has been determined beforehand 
as the value of the peak-value/width ratio. Reference 

10 symbol '^^h"^ denotes the height of a rectangular waveform, 
and the height ^^h" is determined from the value of "i" and 
the value of ^'^K." Further, reference symbol "*^w^' denotes 
the width of a rectangular waveform, and the width ^^w" of 
a rectangular waveform is calculated from the value of ^^h , " 

15 which has been calculated previously, and the value of ^^i . " 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

20 FIG, 12 is a detailed block diagram showing the 

current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 

25 from the event information 5 concerning each cell, the 
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information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
5 designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 

10 instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 6 function. 
Therefore, an EMI component cannot be analyzed through an 

15 FFT operation. Thus far, processing has been identical 
with the current wave calculation operation of the first 
embodiment . 

Reference numeral 30 designates a second rectangular 
waveform height calculation section for calculating the 

20 height of a rectangular waveform such that a ratio of a 
peak value of a rectangular waveform to be modeled to a 
duration of time is constant and such that the area of a 
rectangular waveform becomes equal to an 
instantaneously-changing current value. Reference 

25 numeral 31 designates a second rectangular waveform width 
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calculation section for calculating the width of a 
rectangular waveform such that a ratio of a peak value of 
a rectangular waveform serving as a model of the waveform 
of the instantaneously-changing current to a duration of 
5 time is constant and such that the area of a rectangular 
waveform becomes equal to an instantaneously-changing 
current value. Reference numeral 32 designates a second 
rectangular waveform-shaping section for modeling the 
instantaneously-changing current calculation result 12 as 
10 a rectangular waveform having a constant peak -value/width 



FIG. 13 shows the flowchart of processing pertaining 
to the second rectangular waveform-shaping section, 
including processing pertaining to the second rectangular 

15 waveform height calculation section and processing 
pertaining to the second rectangular waveform width 
calculation section. The following processing is 
iterated on a per-event basis until calculation of a 
current waveform is completed (step 1030) . The height of 

20 a rectangular waveform is expressed by Eq . 1 provided below , 
from the area of a rectangular waveform and "K = h/w." 



h^^iixK)... (1) 



The height ^^h" of a rectangular waveform is calculated 
through use of Eq.l (step 1031) . Reference symbol ''^i" 
25 denotes the value of an instantaneously-changing current 



ratio . 
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of a cell of an event, the event being an object of EMI 
analysis. Reference symbol ^^K" denotes a ratio of a peak 
value of a rectangular waveform to the width of the same 
and assumes a predetermined constant. This processing 
5 corresponds to processing pertaining to the second 

rectangular waveform height calculation section 30. The 
area of a rectangular waveform is defined by ^^w x h." The 
width ^^w" of a rectangular waveform is calculated on the 
basis of "w = i/h" <1032) , Reference symbol ^^h"' denotes 

10 the height "h" of a rectangular waveform calculated 
through processing pertaining to step 1031. This 
processing corresponds to processing pertaining to the 
second rectangular waveform width calculation section 31. 
Next, '''h" is added to I (t+x) until a variable "^x" changes 

15 from -W/2 to W/2 , and At is added to ^^x" (steps 1033 and 
1034). I (t+x) corresponds to the value of electric 
current flowing through all the cells at time t+x. This 
processing corresponds to processing pertaining to the 
second rectangular waveform-shaping section 32. 

20 The current wave calculation result 8 can be obtained 

by means of the flow of a rectangular waveform height 
calculation operation, a rectangular waveform width 
calculation operation, and a rectangular waveform-shaping 
operation. FIG. 14 shows an example current wave 

25 calculation result 8 in the instantaneously-changing 
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current calculation result shown in FIG. 3. The FFT 
processing section 9 subjects the waveform to an FFT 
operation, to thereby determine the frequency 
characteristic of an EMI component. As a result, there 
5 can be obtained an EMI analysis result 10. 

According to the EMI analysis method of the present 
embodiment, the waveform of a current is modeled as a 
rectangular waveform having a constant peak-value/width 
ratio, and the thus-modeled rectangular waveform is 

10 subjected to an FFT operation, thereby enabling more 

accurate analysis of EMI . In an LSI circuit involving flow 
of a large amount of electric current, flow of an electric 
current consumes a long time. The EMI analysis method 
according to the present embodiment enables realization 

15 of a more realistic model and hence enables more accurate 
processing. Further, the EMI analysis method enables EMI 
analysis within a much shorter period of time than the 
transistor-level EMI analysis method. 
(Fourth Embodiment) 

20 A fourth embodiment of the present invention will 

now be described. The present embodiment describes a 
method of analyzing EMI by the steps of: 

calculating an instantaneously-changing current of 
each current from event information concerning each cell, 

25 the information being output from a logic simulator; 
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modeling the instantaneously-changing current as a 
rectangular waveform whose area becomes equal to the 
instantaneously-changing current, through use of a global 
table representing the relationship between an 
5 instantaneously-changing current and the height of a 
rectangular waveform, to thereby calculate a more 
realistic change in current; and 

subjecting a calculation result of the current 
waveform to an FFT operation, to thereby determine the 

10 frequency characteristic of an EMI component of a circuit 
of interest. The global table signifies a table which is 
generally prepared without regard to the type of a cell. 

FIG. 15 shows an instantaneously-changing current 
value and a modeled rectangular waveform according to the 

15 present embodiment. Reference symbol "i'' denotes an 

instantaneously-changing current; ^^t" denotes a current 
processing time; and '*^h" denotes the height of a 
rectangular waveform serving as a model of the waveform 
of the instantaneously-changing current. The height of 

20 a rectangular waveform is calculated, through use of a 
global table representing the relationship between an 
instantaneously-changing current value and the height of 
a rectangular waveform. Reference symbol ^^w" denotes the 
width of a rectangular waveform, and the width ^^w" of a 

25 rectangular waveform is calculated from the value of "^^h," 
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which has been calculated previously, and the value of . " 

FIG. 16 is a graph showing an example global i-h table 
representing the relationship between an 

instantaneously-changing current value and the height of 
5 a rectangular waveform. The table shows the value of "h" 
at each of certain points '^^i." The table indicates that 
"h" assumes a value of 0 at i=0, a value of 5 at i=50, and 
a value of 20 at i=100. In reality, the value of "i" 
changes continuously, and hence the table is expressed as 

10 a graph plotting points provided in the table. The value 
of ^^h" [i.e. , h(i) ] at a certain value of ^^i" is expressed 
by Eq.2 provided below, through linear interpolation of 
the values of two points closest to "i" provided in the 
i-h table: one point has lower values (il, hi) , and the 

15 other point has higher values (12, h2 ) . 

/; = :^i:i^(/_/l) + M ... (2) 
il - i\ 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
20 embodiment (see FIG. 1). 

FIG. 17 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
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calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
5 analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"^^instantaneously-changing current calculation result 

10 12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 

15 waveform of an electric current assumes a 5 function. 

Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. 

20 In the present embodiment, an i-h table 40 has 

described therein the previously-described global table 
representing the relationship between an 

instantaneously-changing current value and the height of 
a rectangular waveform (see FIG. 16) . A third rectangular 
25 waveform height calculation section 41 calculates a peak 
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value of a rectangular waveform serving as a model of an 
instantaneously-changing current, from an 
instantaneously-changing current value and the i-h table 

40. A third rectangular waveform width calculation 

5 section 42 calculates the width of a rectangular waveform 
from a peak value of a rectangular waveform calculated by 
the third rectangular waveform height calculation section 

41, and the instantaneously-changing current value. A 
third rectangular waveform- shaping section 43 models the 

10 instantaneously-changing current calculation result 12 as 
a rectangular waveform, the height of the rectangular 
waveform corresponding to the peak value of the 
rectangular waveform calculated by the third rectangular 
waveform height calculation section 41, and the width of 

15 the rectangular waveform corresponding to the width 
calculated by the third rectangular waveform width 
calculation section 42. 

FIG. 18 is a flowchart showing processing pertaining 
to the third rectangular waveform- shaping section 43, 

20 including processing pertaining to the rectangular 

waveform height calculation section 41 and processing 
pertaining to the rectangular waveform width calculation 
section 42. First, the i-h table 40 is read (step 1040) . 
Next, the subsequent processing operations are iterated 

25 on a per-event basis until the current waveform 
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calculation operation is completed (step 1041) . Two 
points closest to "*^i," (il, hi) and (i2, h2) , are extracted 
from the i-h table 40. Reference symbol "i" is the value 
of an instantaneously-changing current of a cell of an 
5 event, the event being an object of EMI analysis. In the 
table shown in FIG. 16, the two points closest to i=70 are 
(il, hi) = (50, 5) and (12, h2) = (100, 20) (step 1042) . 

The height "*^h" of a rectangular waveform is 
calculated according to Eq.2 (step 1043) . In connection 

10 with the table shown in FIG. 16, in the case where i=70, 
h(70)=ll. This processing; i.e., processing pertaining 
to steps 1042 and 1043, corresponds to processing 
pertaining to the third rectangular waveform height 
calculation section 41. The area of a rectangular 

15 waveform is defined as ^^w x h." The width ^^w" of a 

rectangular waveform is calculated from w=i/h (step 1044) . 
Reference symbol "i" denotes the value of an 
instantaneously-changing current of a cell of an event, 
the event being an object of EMI analysis. Reference 

20 symbol ^^h" denotes the height of a rectangular waveform 
calculated by the third rectangular waveform height 
calculation section 41. In the table shown in FIG, 16, 
when i=70, w=6 . 3 . This processing corresponds to 
processing pertaining to the third rectangular waveform 

25 width calculation section 42 . Next , "h" i s added to I ( t+x) 
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until variable "x" changes from -'W/2 to W/2 , and At is added 
to variable ^^x" (steps 1045 and 1046). I (t+x) denotes 
electric currents flowing through all the cells at time 
^^t+x." This processing corresponds to processing 
5 pertaining to the third rectangular waveform-shaping 
section 43. 

The current waveform calculation result 8 can be 
determined by means of the flow of a rectangular waveform 
height calculation operation, a rectangular waveform 

10 width calculation operation, and a rectangular 

waveform-shaping operation. FIG. 19 shows an example 
current waveform calculation result 8 obtained in a case 
where the table shown in FIG. 16 is used as the i-h table 
40 and where the instantaneously-changing current value 

15 is 70. The rectangular waveform shown in FIG. 19 is 

subjected to an FFT operation by the FFT processing section 
9, to thereby determine the frequency characteristic of 
an EMI component. Thus, there can be obtained an EMI 
analysis result 10. 

20 (Fifth Embodiment) 

An EMI analysis method according to a fifth 
embodiment of the present invention will now be described. 
The present embodiment describes the EMI analysis method 
comprising the steps of: 

25 calculating an instantaneously-changing current of 
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each cell from event information concerning each cell, the 
event information being output from a logic simulators- 
modeling the instantaneously-changing current into 
a rectangular waveform whose area becomes equal to the 
5 instantaneously-changing current, through use of a global 
table representing the relationship among a slew in an 
input waveform, an instantaneously-changing current value , 
and the height of a rectangular waveform, to thereby 
calculate a change in current in consideration of the 

10 more-realistic slew in an input waveform; and 

subjecting the result of calculation of the current 
waveform to an FFT operation, thereby determining the 
frequency characteristic of an EMI component arising in 
a circuit of interest. 

15 The relationship between the instantaneously- 

changing current value and a modeled rectangular waveform 
according to the present embodiment is the same as in the 
case of the rectangular waveform model according to the 
fourth embodiment (shown in FIG. 15) , except that ^^w" and 

20 ^^h" are calculated in consideration of the influence of 
slew in an input waveform. FIG. 20 shows an example global 
i-s-h table representing the relationship among an 
instantaneously-changing current value, the slew in an 
input waveform, and the height of a rectangular waveform, 

25 and a graph plotting the i-s-h table. The table shows the 



64 




value of ^^h^' at a certain point (i , s) , where '*^s" denotes 
the slew (or slew) of an input waveform. In connection 
with this table, '^^h" assumes a value of 0 at i=0 and s=0; 
a value of 5 at i=50 and s=0 ; a value of 14 at i=100 and 
5 s=0 ; a value of 0 at i=0 and s=2 / a value of 6 at i=50 and 
s=2 ; and a value of 20 at i=100 and s=2 . In reality, the 
value of ^^i" and the value of ^^s" change continuously, and 
hence the table is represented by a three-dimensional 
graph plotting points appearing in the table. The value 

10 of h(s,i) at a certain value of ^^s" and a certain value 
of ''^i" is expressed by Eq.3 provided below, through use 
of four points in combination, the four points comprising 
a combination of a low value (il) and a high value (12) 
closest to ^^i" within the i-s-h table and a low value (si) 

15 and a high value (s2) closest to "s"; namely, h(sl, il) , 
h(sl, 12), h(s2, il) , and h(s2, 12). 

[ h(sljl)(s2-s) h{s2Jl)(s - sVj Y i2-i ^ [ h{sl,i2Xs2- s) h{s2,i2){s - sVi Y ~\ 
[ (52-51) ^ (s2-5l) jt/2-/lJ'^[ (52-51) ^ (52-51) JI/2-/1J 

... (3) 

The overall flow of processing pertaining to the EMI 
20 analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 21 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
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present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
5 information being output from the logic simulator 4, 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 

10 changing current (hereinafter referred to as an 

^^instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 

15 As can be seen from FIG. 3, according to the 

instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 

20 identical with the current wave calculation operation of 
the first embodiment. 

The i-s-h table 50 employed in the present embodiment 
is the above-described global table representing the 
relationship among an instantaneously-changing current 

25 value, the slew in an input waveform, and the height of 
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a rectangular waveform (FIG. 20) . A fourth rectangular 
waveform height calculation section 52 calculates a peak 
value of a rectangular waveform serving as a model of an 
instantaneously-changing current, from an 
5 instantaneously-changing current value of an instance 
serving as an object of EMI analysis, the slew in an input 
waveform, and the i-s-h table 50, all of which are read 
from Slew data 51. The Slew data 51 are provided in the 
netlist 1 and describe information concerning the slew in 

10 a waveform input to each instance. Reference numeral 53 
designates a fourth rectangular waveform width 
calculation section for calculating the width of a 
rectangular waveform, through use of the peak value 
calculated by the fourth rectangular waveform height 

15 calculation section 52 and the instantaneously-changing 
current value such that the area of the rectangular 
waveform becomes equal to the instantaneously-changing 
current value. Reference numeral 54 designates a fourth 
rectangular waveform- shaping section for modeling the 

20 instantaneously-changing current calculation result 12 as 
a rectangular waveform, the height of the rectangular 
waveform corresponding to the peak value calculated by the 
fourth rectangular waveform height calculation section 52 , 
and the width of the rectangular waveform corresponding 

25 to the width calculated by the fourth rectangular waveform 
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width calculation section 53. 

FIG. 22 shows the flowchart of processing pertaining 
to the fourth rectangular waveform- shaping section 54, 
including processing pertaining to the rectangular 
5 waveform height calculation section 52 and processing 
pertaining to the rectangular waveform width calculation 
section 53. First, the i-s-h table 50 is read (step 1050) , 
and the Slew data 51 are read (step 1051) . The processing 
provided below is iterated on a per-event basis until 

10 calculation of a current waveform is completed (step 1052) . 
From the i-s-h table 50, there are extracted four points 
which are combinations of a low value closest to '^^i," a 
high value closest to "i,''' a low value closest to ^^s," and 
a high value closes t to ^^s" ; for example, h(sl, il) , h(sl, 

15 12) , h (s2 , il) , and h (s2 , 12) . Here, reference symbol "i" 
denotes the value of an instantaneously-changing current 
of a cell of an event, the event being an object of EMI 
analysis. Further, reference symbol "s" denotes the slew 
in an input waveform pertaining to the cell. The four 

20 points are extracted on the basis of the name of the 
instance to be processed and the Slew data 51. 

In connection with the table shown in FIG. 20, in 
the case of i=70 and s=l , there are extracted four values; 
25 h(0, 50)=5, h(0, 100)=14, h(2, 50)=6, and h(2, 100)=20 
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(step 1053) . The height "h" of the rectangular waveform 
is calculated by Eq. 3 (step 1054) . In connection with the 
table shown in FIG. 20, in the case where i=70 and s=l , 
h(l, 70) =10.1. Processing pertaining to steps 1053 and 
5 1054 correspond to processing pertaining to the fourth 
rectangular waveform height calculation section 52. 

The area of the rectangular waveform is defined by 
''^w X h." The duration (w) of a rectangular waveform is 
calculated on the basis of w=i/h (step 1055) . Here, ^^i" 
10 denotes the value of an instantaneously-changing current 
of a cell of an event, the event being an object of EMI 
analysis, and «h" denotes the height of a rectangular 
waveform calculated by the fourth rectangular waveform 
height calculation section 52 . 

15 

In connection with the table shown in FIG. 20, in 
the case of i=70 and s=l , w=6 . 9 . This processing 
corresponds to processing pertaining to the fourth 
rectangular waveform width calculation section 53. Next, 

20 ^^h" is added to I (t+x) until variable ^^x" changes from -W/2 
to W/2 , and At is added to ^^x" (steps 1056 and 1057) . Here, 
I (t+x) represents total current flowing through all the 
cells at time ^^t+x." This processing corresponds to 
processing pertaining to the fourth rectangular 

25 waveform-shaping section 54 . 
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The current waveform calculation result 8 can be 
determined by means of the flow of a rectangular waveform 
height calculation operation^ a rectangular waveform 
width calculation operation, and a rectangular 
5 waveform-shaping operation. FIG. 23 shows an example 
current waveform calculation result 8 obtained in a case 
where the i-s-h table 50 corresponds to the table shown 
in FIG. 20, the instantaneously-changing current value is 
70, and the slew in an input waveform assumes a value of 

10 1. The rectangular waveform is subjected to an FFT 

operation by the FFT processing section 9, to thereby 
determine the frequency characteristic of an EMI component. 
Thus, there can be obtained an EMI analysis result 10. 
(Sixth Embodiment) 

15 An EMI analysis method according to a sixth 

embodiment of the present invention will now be described. 
The present embodiment describes the EMI analysis method 
comprising the steps of: 

calculating an instantaneously-changing current of 

20 each cell from event information concerning each cell, the 
event information being output from a logic simulators- 
modeling the instantaneously-changing current as a 
rectangular waveform whose area becomes equal to the 
instantaneously-changing current, through use of a global 

25 table representing the relationship among an output load 
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capacitance, an instantaneously-changing current value, 
and the height of a rectangular waveform, to thereby 
calculate a change in current in consideration of the 
more-realistic output load capacitance; and 
5 subjecting the result of calculation of the current 

waveform to an FFT operation, thereby determining the 
frequency characteristic of an EMI component arising in 
a circuit of interest. 

The relationship between the ins tantaneously- 

10 changing current value and a modeled rectangular waveform 
according to the present embodiment is the same as in the 
case of the rectangular waveform model according to the 
fourth embodiment (shown in FIG. 15) , except that "w" and 
"h" are calculated in consideration of the influence of 

15 slew in an input waveform. FIG. 24 shows an example global 
i-c-h table representing the relationship among an 
instantaneously-changing current value, an output load 
capacitance, and the height of a rectangular waveform, and 
a graph plotting the i-c-h table. The table shows the 

20 value of "h" at a certain point (i , c) , where ''^c" denotes 
an output load capacitance . In connection with this table , 
^^h" assumes a value of 0 at i=0 and c=0 ; a value of 5 at 
i=50 and c=0 ; a value of 20 at i = 100 and c=0 ; a value of 
0 at i=0 and c=20; a value of 8 at i=50 and c=20; and a 

25 value of 25 at i=100 and c=20. In reality, the value of 
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"i" and the value of "c'' change continuously, and hence 
the table is represented by a three-dimensional graph 
plotting points appearing in the table. The value of 
h{c,i) at a certain value of "c" and a certain value of 
5 ^^i" is expressed by Eq,4 provided below, through use of 
four points in combination, the four points comprising a 
combination of a low value (il) and a high value (i2) 
closest to ^^i" within the i-c-h table and a low value (cl) 
and a high value (a2) closest to ^^c"; namely, h(cl, il) , 
10 h(cl, ±2), h(c2, il) , andh(c2, 12). 

y--.l2L!lx+h ... (4) 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 

15 embodiment (see FIG. 1) . 

FIG. 25 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 

20 calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 




conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously-changing current calculation result 
5 12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 

10 waveform of an electric current assumes a 5 function. 

Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. 

15 The i-c-h table 60 employed in the present embodiment 

is the above-described global table representing the 
relationship between an instantaneously-changing current 
value and the height of a rectangular waveform (FIG. 24) . 
Reference numeral 62 designates a fifth rectangular 

20 waveform height calculation section which calculates a 
peak value of a rectangular waveform serving as a model 
of an instantaneously-changing current, from an 
instantaneously-changing current value of an instance 
serving as an object of EMI analysis and an output load 

25 capacitance. Reference numeral 63 designates a fifth 
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rectangular waveform width calculation section for 
calculating the width of a rectangular waveform, through 
use of the peak value calculated by the fifth rectangular 
waveform height calculation section 62 and the 
5 instantaneously-changing current value and such that the 
area of the rectangular waveform becomes equal to the 
instantaneously-changing current value. Reference 
numeral 64 designates a fifth rectangular waveform- 
shaping section for modeling the instantaneously-changing 

10 current calculation result 12 as a rectangular waveform, 
the height of the rectangular waveform corresponding to 
the peak value calculated by the fifth rectangular 
waveform height calculation section 62, and the width of 
the rectangular waveform corresponding to the width 

15 calculated by the fifth rectangular waveform width 
calculation section 63. 

FIG. 2 6 shows the flowchart of processing pertaining 
to the fifth rectangular waveform-shaping section 64 , 
including processing pertaining to the rectangular 

20 waveform height calculation section 62 and processing 
pertaining to the rectangular waveform width calculation 
section 63. First, the i-c-h table 60 is read (step 1060) , 
and capacitance data 61 are read (step 1061) . The 
processing provided below is iterated on a per-event basis 

25 until calculation of a current waveform is completed (step 
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1062) . From the i-c-h table 60, there are extracted four 
points which are combinations of a low value closest to 
^^i," a high value closest to "i," a low value closest to 
"c," and a high value closest to ^^c" ; for example, h(cl, 
5 il), h(cl, i2), h(c2, il), and h (c2 , i2) . Here , reference 
symbol ^^i" denotes the value of an instantaneously- 
changing current of a cell of an event, the event being 
an object of EMI analysis. Further, reference symbol ^^c" 
denotes an output load capacitance pertaining to a cell 

10 of an event, the event being an object of EMI analysis. 
The four points are extracted on the basis of the name of 
instance to be processed and the capacitance data 61. 

In connection with the table shown in FIG. 24, in 
the case where i=70 and c=10, there are extracted four 

15 values; h(0, 50)=5, h(0, 100)=20, h(20, 50)=8, andh(20, 
100)=25 (step 1063). The height "h" of the rectangular 
waveform is calculated by Eq . 3 (step 1064) . In connection 
with the table shown in FIG. 24, in the case where i=70 
andc=10, h(10, 70) =12. 9. Processing pertaining to steps 

20 1063 and 1064 corresponds to processing pertaining to the 
fifth rectangular waveform height calculation section 62 . 
The area of the rectangular waveform is defined by "w x 
h." The duration (w) of a rectangular waveform is 
calculated on the basis of w=i/h (step 1065) . Here, ^^i" 

25 denotes the value of an instantaneously-changing current 
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of a cell of an event, the event being an object of EMI 
analysis, and ^^h" denotes the height of a rectangular 
waveform calculated by the fifth rectangular waveform 
height calculation section 62 . 
5 In connection with the table shown in FIG. 24, in 

the case where i=70 and c=10, ws=5 . 4 . This processing 
corresponds to processing pertaining to the fifth 
rectangular waveform width calculation section 53, Next, 
"h" is added to I (t+x) until variable '^^x" changes from -W/2 

10 to W/2 , and At is added to '^x'Msteps 1066 and 1067) . Here, 
I (t+x) represents total current flowing through all the 
cells at time ^^t+x." This processing corresponds to 
processing pertaining to the fifth rectangular 
waveform- shaping section 64 . 

15 The current waveform calculation result 8 can be 

determined by means of the flow of a rectangular waveform 
height calculation operation, a rectangular waveform 
width calculation operation, and a rectangular 
waveform-shaping operation. FIG. 27 shows an example 

20 current waveform calculation result 8 obtained in a case 
where the i-c-h table 60 corresponds to the table shown 
in FIG. 24, the instantaneously-changing current value is 
70, and an output load capacitance assumes a value of 10. 
The rectangular waveform is subjected to an FFT operation 

25 by the FFT processing section 9, to thereby determine the 
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frequency characteristic of an EMI component. Thus, 
there can be obtained an EMI analysis result 10. 

The EMI analysis method of the present embodiment 
forms a rectangular waveform model from a global table 
5 defining the relationship between an output load 

capacitance, an instantaneously-changing current value, 
and the height of a rectangular waveform, thereby enabling 
preparation of a more realistic model and more accurate 
processing . 
10 (Seventh Embodiment) 

An EMI analysis method according to a seventh 
embodiment of the present invention will now be described. 
FIG. 28 is a block diagram showing the overall flowchart 
of the EMI analysis method according to the present 
15 embodiment. The EMI analysis method according to the 
present embodiment is identical with the EMI analysis 
method according to the first embodiment described in 
connection with FIG. 1, except for addition of data A9 and 
AlO, a voltage drop calculation/instantaneously-changing 
20 current correction section All, and an 

instantaneously-changing current calculation result A6 
output in step All. Those portions of the EMI analysis 
method according to the present embodiment which are 
identical with those of the EMI analysis method according 
25 to the first embodiment are assigned the same reference 
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numerals, and their repeated explanations are omitted. 
A9 designates data concerning the layout of a 
semiconductor integrated circuit which is an object of EMI 
analysis. AlO designates an I-V table file (or an I-V 
5 function) representing the relationship between the 

amount of electric current and the degree of a drop in line 
voltage. All designates a section of calculating a drop 
in line voltage and correcting the result of an 
instantaneously-changing current calculated in step 11 
10 (hereinafter referred to as a ^^voltage drop 

calculation/instantaneously-changing current correction 
section") . A6 designates the result of an 
instantaneously-changing current being calculated by the 
voltage drop calculation/instantaneously-changing 
15 current correction section All in consideration of a drop 
in line voltage. 

FIG. 2 9 is a block diagram showing the flow of 
processing of the voltage drop 

calculation/instantaneously-changing current correction 
20 section All according to the seventh embodiment of the 
present invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 
section All will now be described by reference to FIG. 29. 
Since data 12 have already been described in detail in 
25 connection with the first embodiment, the same reference 
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numeral is assigned to the data 12, and repeated 
explanation is omitted. The voltage drop 
calculation/instantaneously-changing current correction 
section All comprises a power net extraction section B2 , 
5 a section B6 for calculating the resistance of each segment 
(hereinafter referred to as a "^^segment resistance 
calculation section B6") , a section B3 for calculating an 
electric current flowing through each segment 
(hereinafter referred to as a "segment current calculation 

10 section B3") , a section B7 for c^culating a drop in line 
voltage of each segment (hereinafter referred to as a 
^^segment voltage drop calculation section B7") , a section 
BIO for calculating a drop in line voltage for each 
instance (hereinafter referred to as an "instance voltage 

15 drop calculation section BIO") , and a current correction 
section B12 . 

These steps will now be described. In step B2 , 
layout data A9 are entered, andr the geometry of a power 
line is determined in consideration of information 

20 concerning electrical interconnections. On the basis of 
the thus -determined geometry of a power line, nodes, 
change in line width, bends/turns, and switching to 
another interconnection are identified, and the power line 
is branched in accordance with identification results, 

25 thereby forming segment information B4 pertaining to 
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respective segments of the power line. The thus- 
constructed database for the segment information B4 
includes the length and width of a line in each segment. 
FIG. 30 is a schematic illustration shoving the 
5 layout data A9 which is input to be entered in step B2 . 
FIG. 31 is a schematic illustration showing for the data 
base the segment information B4 produced in step B2 . 
Segments 1 through 12 formed as a result of branching of 
an interconnection are represented by Sgl through Sgl2. 

10 The segment information B4 concerning respective segme'nts 
of the power line, the information being produced in ste^p 
B2 , is entered in step B6, and resistance information B8 
concerning the resistance of a line in each segment is 
produced from information concerning the length and width 

15 of a line in each segment and a predetermined resistance 
parameter assigned to each interconnection layer. 

FIG. 32 is a schematic representation showing a data 
base for the resistance information B8 produced in step 
B6 . In this drawing, the resistance of segment Sgl of a 

20 power line (VDD) is indicated by Rl ; the resistance of 
segment Sg2 of the power line is indicated by R2 ; the 
resistance of segment Sg3 of the power line is indicated 
by R3 ; the resistance of segment Sg4 of the power line is 
indicated by R4 / the resistance of segment Sg5 of the power 

25 line is indicated by R5 ; and the resistance of segment Sg6 
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of the power line is indicated by R6 , In step B3 , there 
are entered information Bl concerning an 
instantaneously-changing current of each instance 
produced on the basis of the instantaneously-changing 
5 current calculation result 12 (hereinafter referred to 
simply as ^^instantaneously-changing current information 
Bl") and the segment information B4 concerning the 
respective segments of a power line produced in step B2 . 
The electric current of each instance is distributed into 

10 current components which flow through the respective 

segments of the power line . The current components of each 
instance, which components have been distributed into the 
respective segments of the power line, are summed, to 
thereby produce current information B5 for each segment. 

15 FIG. 33 is a schematic representation showing the 

data base for the current information B5 produced in step 
B3 . In the drawing, electric currents flowing through 
instances 1 through 3 are indicated by II through 13, and 
electrical currents flowing through segments Sgl through 

20 Sg6 are indicated by Isl through Is6. In this case, 

131=11+12+13, Is2=12+I3, Is3=I43, Is4=Il , Is5=I2, and 
Is6=I3. In step B7 , there are entered the current 
information B5 prepared each of the segments in step B3 
and the resistance information B8 prepared for each of the 

25 segments in step B6 . Information B9 concerning a drop in 
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line voltage in each segment (hereinafter referred to 
simply as ^^voltage drop information B9") is prepared 
according to Ohm's low. 

FIG. 34 is a schematic representation showing the 
5 data base for the voltage drop information B9 prepared in 
step B7 . In the drawing, voltage drops in respective 
segments Sgl through Sg6 are indicated by AVsl through AVs6. 
In this case, AVsl=RlxIsl; AVs2=R2xIs2; AVs3=R3xIs3 ; 
AVs4=R4xIs4; AVs5=R5xIs5; and AVs6=R6xIs6. In step BIO, 

10 there is entered the voltage drop information B9 prepared 
for each of the segments in step B7 . Voltage drops in 
respective segments of a current flow path of each instance 
are summed, to thereby produce voltage drop information 
Bll for each instance. 

15 FIG. 35 is a schematic representation showing the 

data base for the voltage drop information Bll prepared 
in step BIO. In the drawing, voltage drops arising in 
instances 1 through 3 are indicated by AVl through AV3. 
In this case, AVl=AVsl+AVs4 ; AV2=AVsl+AVs2+AVs5 ; and 

20 AV3=AVsl+AVs2+AVs3+AVs6 . Here, Avi is expressed as 
AVl=Rlx (I1+I2+I3) +R4xll ; AV2 is expressed as 
AV2=Rlx (I1+I2+I3) +R2x (12+13) +R5xI2 ; and AV3 is expressed 
as AV3=Rlx (I1+I2+I3) +R2x (12+13) +R3xI3+R6xl3 . 

In step B12 , there are entered the voltage drop 

25 information Bll prepared for each instance in step BIO and 
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an I-V table file AlO (or an I-V function) representing 
the relationship between the amount of electric current 
and the degree of drop in line voltage, wherewith the 
instantaneously-changing current calculation result 12 is 
5 corrected. As a result, there is produced an 

instantaneously-changing current calculation result A6 
for each instance, the calculation result taking into 
consideration a drop in line voltage which is dependent 
on the resistance of a power line. 

10 As mentioned above, the EMI analysis method 

according to the seventh embodiment takes into 
consideration a drop in line voltage which is dependent 
on the resistance of a power line, to thereby correct the 
instantaneously-changing current calculation result 12 

15 obtained in the first embodiment in consideration of a 
physical phenomena, such as fluctuations in a line voltage 
of a real semiconductor integrated circuit. Therefore, 
the EMI analysis method of the present embodiment can 
produce a more accurate model of a current waveform, 

20 thereby reducing undesired higher harmonic components 
included in the frequency characteristic of a circuit to 
be analyzed, the characteristic being determined by means 
of subjecting the current waveform to FFT processing. 

Although the seventh embodiment has referred to the 

25 first embodiment, similar results are yielded even when 
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the technical concept of the present embodiment is applied 
to any of the second through sixth embodiments. 
(Eighth Embodiment) 

An EMI analysis method according to an eighth 
5 embodiment of the present invention will now be described. 
FIG. 36 is a block diagram showing the overall flowchart 
of the EMI analysis method according to the present 
embodiment. The EMI analysis method according to the 
present embodiment is identical with the EMI analysis 
10 method according to the seventh embodiment described in 
connection with FIG. 28, except for addition of data CI 
and replacement of the voltage drop 

calculation/instantaneously-changing current correction 
section All by a voltage drop 
15 calculation/instantaneously-changing current correction 
section C2 . Those portions of the EMI analysis method 
according to the present embodiment which are identical 
with those of the EMI analysis method according to the 
first embodiment are assigned the same reference numerals , 
20 and their repeated explanations are omitted. 

CI designates a C-V table file (or a C-V function) 
representing the relationship between decoupling 
capacitance existing between a power and the ground and 
the degree of a drop in line voltage. C2 designates a 
25 section for calculating a drop in line voltage and 
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correcting the result of an instantaneously-changing 
current being calculated in step 11 (hereinafter referred 
to as a ''Woltage drop calculation/ins tantaneously- 
changing current correction section") . FIG. 37 is a block 
5 diagram showing the flow of processing of the voltage drop 
calculation/instantaneously-changing current correction 
section C2 according to the eighth embodiment of the 
present invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 
10 section C2 will now be described by reference to the 
accompanying drawings. Those portions of the EMI 
analysis method according to the present embodiment which 
are identical with those of the EMI analysis method 
according to the seventh embodiment are assigned the same 
15 reference numerals, and their repeated explanations are 
omitted. Since the data 12 have already been described 
in detail in connection with the first embodiment, the same 
reference numeral is assigned to the data 12, and its 
repeated explanation is omitted. 
20 The voltage drop calculation/instantaneously- 

changing current correction section C2 is configured such 
that a section Dl for calculating the capacitance of each 
segment in a power- to-ground line (hereinafter referred 
to simply as a ^'capacitance calculation section Dl") is 
25 added to the voltage drop calculation/ins tantaneously- 
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changing current correction section All, and the segment 
voltage drop calculation section B7 is changed to D3 . 
Those sections and their corresponding steps will now be 
described. In step Dl, there is entered the segment 
5 information B4 concerning a power line prepared in step 
B2 in the seventh embodiment. Information D2 concerning 
the capacitance of each segment in the power- to-ground 
line is formed from the information concerning the area 
of an overlap existing between the segments in the 
10 power-to-ground line, the information being included in 
the segment information B4 , as well as from a predetermined 
capacitance parameter assigned to each interconnection 
layer . 

FIG. 38 is a schematic illustration showing the data 
15 base for the capacitance information D2 prepared in step 
Dl . In the drawing, a capacitor provided between segments 
Sgl and SgV is indicated by CI; a capacitor provided 
between segments Sg2 and Sg8 is indicated by C2 ; and a 
capacitor provided between segments SgS and Sg9 is 
20 indicated by C3 . In step D3 , there are entered the current 
information B5 prepared for each segment in step B3 in the 
seventh embodiment, the resistance information BB 
prepared for each segment in step B6 in the seventh 
embodiment, the C-V table file CI (or C-V function) 
25 representing the relationship between the decoupling 
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capacitance existing between the power and the ground and 
the degree of drop in line voltage, and the capacitance 
information D2 prepared in step Dl with regard to the 
capacitance of each of the segments in the power-to-ground 
5 line. Those information pieces are corrected according 
to Ohm's law and the C-V table, to thereby produce the 
voltage drop information B9 for each segment so as to 
reflect decoupling capacitance as well as the resistance 
of a power line. 

10 FIG. 39 is a schematic illustration showing the 

voltage drop information B9 prepared in step D3 . In the 
drawing, a voltage drop arising in segment Sgl is indicated 
by AVsll; a voltage drop arising in segment Sg2 is 
indicated by AVs21; a voltage drop arising in segment Sg3 

15 is indicated by AVs31; a voltage drop arising in segment 
Sg4 is indicated by AVs41; a voltage drop arising in 
segment Sg5 is indicated by AVs51; and a voltage drop 
arising in segment Sg is indicated by AVs61 . In this case, 
there are obtained AVsll=f (CI, AVsl) , AVs21=f (C2, AVs2) , 

20 AVs31=f(C3, AVs3) , AVs41=AVs4, AVs51=AVs5, and 
AVs61=AVs6 . 

As has been described, the EMI analysis method 
according to the eighth embodiment newly adopts correction 
of a current value in response to fluctuations in a line 
25 voltage, which would be induced by decoupling capacitance 
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existing between the power and the ground, in addition to 
correction of the current value in consideration of a drop 
in the line voltage, which would be induced by the 
resistance of a power line. The instantaneously-changing 
5 current result 12 obtained for each instance in the first 
embodiment can be corrected in consideration of physical 
phenomena, such as fluctuations in a line voltage of a real 
semiconductor integrated circuit, thereby enabling 
realization of more accurate modeling of a current 
10 waveform. Therefore, the EMI analysis method of the 

present embodiment can produce a more accurate model of 
current waveform, thereby reducing undesired higher 
harmonic components included in the frequency 
characteristic of a circuit to be analyzed, the 
15 characteristic being determined by means of subjecting the 
current waveform to FFT processing. 

Although the eighth embodiment has referred to the 
first embodiment, similar results are yielded even when 
the technical concept of the present embodiment is applied 
20 to any of the second through sixth embodiments. 
(Ninth Embodiment) 

An EMI analysis method according to a ninth 
embodiment of the present invention will now be described. 
FIG. 40 is a block diagram showing the overall flowchart 
25 of the EMI analysis method according to the present 



88 




embodiment. The EMI analysis method according to the 
present embodiment is identical with the EMI analysis 
method according to the seventh embodiment described in 
connection with FIG. 28, except that the voltage drop 
5 calculation/instantaneously-changing current correction 
section All is replaced by a voltage drop 
calculation/instantaneously-changing current correction 
section El. Those portions of the EMI analysis method 
according to the present embodiment which are identical 

10 with those of the EMI analysis method according to the 
first embodiment are assigned the same reference numerals , 
and their repeated explanations are omitted. El 
designates a section for calculating a drop in line voltage 
and correcting the result of an instantaneously-changing 

15 current calculation performed in step 11 (hereinafter 
referred to as a ^^voltage drop 

calculation/instantaneously-changing current correction 
section") . 

FIG. 41 is a block diagram showing the flow of 

20 processing of the voltage drop 

calculation/instantaneously-changing current correction 
section El according to the ninth embodiment of the present 
invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 

25 section El will now be described by reference to the 
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accompanying drawings . Those portions of the EMI 
analysis method according to the present embodiment which 
are identical with those of the EMI analysis methods 
according to the seventh and eighth embodiments are 
5 assigned the same reference numerals, and their repeated 
explanations are omitted. 

Since the data 12 have already been described in 
detail in connection with the first embodiment, the same 
reference numeral is assigned to the data 12, and repeated 
10 explanation is omitted. The voltage drop 

calculation/instantaneously-changing current correction 
section El is configured by means of removing, from the 
voltage drop calculation/instantaneously-changing 
current correction section C2 employed in the eighth 
15 embodiment, the segment current calculation section B3 for 
calculating an electric current flowing through each 
segment, the voltage drop calculation section D3 for 
calculating a drop in the voltage of each segment, and the 
instance voltage drop calculation section BIO for 
20 calculating a drop in the voltage of each instance; and 
newly adding a transient analysis section Fl to the same. 

These sections and their corresponding steps will 
now be described. In step Fl , there are entered the 
instantaneously-changing current information Bl prepared 
25 for each instance from the instantaneously-changing 
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current calculation result 12 in step 11, the resistance 
information B8 prepared for each segment in step B6 in the 
eighth embodiment, and the capacitance information D2 
concerning a capacitor provided between segments prepared 
5 in step Dl in the eight embodiment. Through a simulation 
performed by a transient analysis simulator typified by 
SPICE, there is produced voltage drop information Bll 
concerning a voltage drop of each instance. 

As has been described, under the EMI analysis method 
10 according to the ninth embodiment, the resistance of a 
power line and fluctuations in a line voltage, which would 
be induced by decoupling capacitance existing between the 
power line and ground, are determined through use of a 
transient analysis simulator. Physical phenomena, such 
15 as fluctuations in a line voltage of a real semiconductor 
integrated circuit, can be faithfully reflected n the 
instantaneously-changing current result 12 obtained for 
each instance in the first embodiment. Therefore, the EMI 
analysis method of the present embodiment enables 
20 materialization of highly-accurate EMI analysis as 

described in the background art. Further, there can be 
realized high-speed EMI analysis which can be applied to 
a large-level circuit. 

Although the ninth embodiment has referred to the 
25 first embodiment, similar results are yielded even when 
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the technical idea of the present embodiment is applied 
to any of the second through sixth embodiments . 
(Tenth Embodiment) 

An EMI analysis method according to a tenth 
5 embodiment of the present invention will now be described. 
In the second through ninth embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 
the information being output from a logic simulator, and 
10 the thus-calculated instantaneously-changing current is 
modeled as a rectangular waveform, thereby determining 
change in electric current. However, so long as the 
waveform of an instantaneous electric current is modeled 
as a triangular wave, a more realistic change in electric 
15 current can be determined. By way of example, there will 
now be described an EMI analysis method comprising the 
steps of: 

calculating a peak value such that an 
instantaneously-changing electric current becomes equal 
20 to the area of a triangular wave ; 

modeling an instantaneously-changing electric 
current as a triangular wave having the peak value and a 
certain width, to thereby calculate a more realistic 
change in electric current; and 
25 subjecting the result of calculation of an electric 
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waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
of interest. 

FIG. 42 shows an instantaneous electric current and 
5 a modeled triangular wave according to the present 
embodiment. In the drawing, ^^i" denotes an 
instantaneously-changing current; and '^'t"^ denotes a 
current processing time. Further, ''^W" denotes the width 
of a triangular wave serving as a model of an 

10 instantaneously-changing current and corresponds to a 
constant. An appropriate value has been determined as a 
constant beforehand. Reference symbol "h" denotes the 
height of a triangular wave and is calculated on the basis 
of ^^i" and "*^W . " An electric current at time t+x is 

15 expressed by Eq.5 provided below. 



... (5) 

Further, an electric current at time t-x is also expressed 
by Eq.5. 



method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 43 is a detailed block diagram showing the 
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The overall flow of processing of the EMI analysis 



93 




current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
5 from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 

10 designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
^^instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 

15 instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 

20 FFT operation. Thus far, the processing has been 

identical with the current wave calculation operation of 
the first embodiment. Reference numeral 100 denotes a 
triangular waveform height calculation section for 
calculating the height of a triangular wave such that the 

25 area of a triangular wave serving as a model of an 
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instantaneously-changing current becomes equal to an 
instantaneously-changing current. Reference numeral 101 
denotes a triangular-waveform-shaping section for 
modeling the instantaneously-changing current result 12 
5 into a triangular wave having a constant width. 

FIG. 44 is a flowchart of processing pertaining to 
the triangular-waveform-shaping section 101, including 
processing pertaining to the triangular waveform height 
calculation section 100. The following processing is 
10 iterated on a per-event basis until calculation of a 

current waveform is completed (step 1100) . The area of 
a triangular wave is Wxh/2 and assumes a value of an 
instantaneously-changing current ^^i." The height ^^h'' of 
a triangular wave is calculated from 2xi/W (step 1101) . 
15 Reference symbol "i" denotes an instantaneously-changing 
current of a cell of an event, the event being an object 
of EMI analysis. This processing corresponds to 
processing pertaining to the triangular waveform height 
calculation section 100. Until variable ^^x" changes from 
20 0 toW/2, h(c, i) expressed by Eq. 5 is added to I (t+x) and 
I(t-x) . Further, At is added to variable ''x" (steps 1102 
and 1103). Here, I (t+x) denotes total electric current 
flowing through all the cells at time t+x, and I (t-x) 
denotes total electric current flowing through all the 
25 cells at time t-x. This processing corresponds to 
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processing pertaining to the triangular-waveform-shaping 
section 101. 

The current wave calculation result 8 can be obtained 
by means of the flow of a triangular waveform height 
5 calculation operation and a triangular-waveform-shaping 
operation, FIG. 45 shows an example current wave 
calculation result 8 at time W=5 in the 
instantaneously-changing current calculation result 
shown in FIG. 3. The FFT processing section 9 subjects 
10 the waveform to an FFT operation, to thereby determine the 
frequency characteristic of an EMI component. As a result, 
there can be obtained an EMI analysis result 10, 
(Eleventh Embodiment) 

An EMI analysis method according to an eleventh 
15 embodiment of the present invention will now be described. 
In the second through ninth embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 
the information being output from a logic simulator, and 
20 the thus-calculated instantaneously-changing current is 
modeled as a rectangular waveform, thereby determining 
change in electric current. However, so long as the 
waveform of an instantaneous electric current is modeled 
as a multi-order-function wave, a more realistic change 
25 in electric current can be determined. By way of example, 
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there will now be described an EMI analysis method 
comprising the steps of: 

calculating a peak value such that an 
instantaneously-changing electric current becomes equal 
5 to the area of a triangular wave; 

modeling an instantaneously-changing electric 
current as a multi-order-function having the peak value 
and a given width, to thereby calculate a more realistic 
change in electric current; and 
10 subjecting the result of calculation of an electric 

waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
of interest. 

[0135] 

15 FIG. 4 6 shows an instantaneous electric current and 

a modeled multi -order-function according to the present 
embodiment. In the drawing, "i" denotes an 
instantaneously-changing current; and "t" denotes a 
current processing time. Further, "W" denotes the width 

20 of a multi-order-function wave serving as a model of an 
instantaneously-changing current and corresponds to a 
constant. An appropriate value has been determined as a 
constant beforehand. Reference symbol "h" denotes the 
height of a multi-order-function wave and is calculated 

25 on the basis of '^^i" and '''^W . " A multi -order-function wave 
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serving as a model in the present embodiment is expressed 
by Eq.6 provided below. 

y = a(x-Wyix + W)\.. (6) . 
The width of a multi-order-function wave model is taken 
5 as '^2 X W." 

According to Eq.6, the height ^^h" of a multi- 
order-function wave model is determined when x=0 and is 
expressed by Eq.7 provided below. 
h = axW^ ... (7) . 

10 The area of the multi-order-function wave model is 

expressed by Eq.8 provided below. 

S^^axW'... (8) . 

Since the thus-defined area becomes equal to the 
instantaneously-changing current value "i," "a" is 
15 expressed though use of ^^W" and ^^i" as Eq. 9 provided below. 

a = — (9) . 

\6xW' 

As a result of calculation of ^^a" through use of Eq. 9 , 
the height ''^h" of a mul ti -order-function wave serving as 
a model is determined. Let W=2 . 5 and i=100, then a=0.96 
20 is derived from Eq.9. Eq.7 consequently gives h=37.5. 
Further, let W=2 . 5 and i=50 , then a=0.48 is derived from 
Eq.9. Eq.7 consequently gives h=18.75. 

The overall flow of processing of the EMI analysis 
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method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 47 is a detailed block diagram showing the 
5 current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
10 information being output from the logic simulator 4. 

Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
15 changing current (hereinafter referred to as an 

"instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. 

FIG. 3 shows an example instantaneously-changing 
20 current calculation result 12. As can be seen from FIG. 
3, according to the instantaneously-changing current 
calculation result, the waveform of an electric current 
assumes a 6 function. Therefore, an EMI component cannot 
be analyzed through an FFT operation. Thus far, the 
25 processing has been identical with the current wave 
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calculation operation of the first embodiment. 

Reference numeral 110 denotes a multi-order- 
function waveform height calculation section for 
calculating the height of a multi-order-function wave such 
5 that the area of a multi -order-function wave serving as 
a model of an instantaneously-changing current becomes 
equal to an instantaneously-changing current. Reference 
numeral 111 denotes a multi -order-function waveform- 
shaping section for modeling the instantaneously-changing 
10 current result 12 as a multi-order-function wave having 
a constant width. 
[0142] 

FIG. 4 8 is a flowchart of processing pertaining to 
the triangular-waveform-shaping section 101^ including 

15 processing pertaining to the triangular waveform height 
calculation section 100. The following processing is 
iterated on a per-event basis until calculation of a 
current waveform is completed {step 1110) . A peak- 
value-dependent coef f icient ''^a" of a multi-order-function 

20 is cal culated f rom Eq . 9 (step 1111). Reference symbol "i" 
denotes an instantaneously-changing current of a cell of 
an event, the event being an object of EMI analysis, and 
■^^W" is a constant which is dependent on the width of a wave. 
This processing corresponds to processing pertaining to 

25 the multi-order-function waveform height calculation 
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section 110. 

Until variable '*^x" changes from -W to W, multi- 
order-function '^y" expressed by Eq.6 is added to I (t+x) . 
Further, At is added to variable ^^x" (steps 1112 and 1113) . 
5 Here, I (t+x) denotes total electric current flowing 
through all the cells at time t+x. This processing 
corresponds to processing pertaining to the multi- 
order-f unction waveform- shaping section 111. 

The current wave calculation result 8 can be obtained 

10 by means of the flow of a multi-order-function waveform 
height calculation operation and a multi-order-function 
waveform- shaping operation. FIG. 4 9 shows an example 
current wave calculation result 8 at time W=2 . 5 in the 
instantaneously-changing current calculation result 

15 shown in FIG. 3. The FFT processing section 9 subjects 
the waveform to an FFT operation, to thereby determine the 
frequency characteristic of an EMI component . As a result, 
there can be obtained an EMI analysis result 10. 
(Twelfth Embodiment) 

20 An EMI analysis method according to a twelfth 

embodiment of the present invention will now be described. 
In the second through ninth embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 

25 the information being output from a logic simulator, and 
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the thus-calculated instantaneously-changing current is 
modeled as a rectangular waveform, thereby determining 
change in electric current. However, so long as the 
waveform of an instantaneous electric current is modeled 
5 as a Gauss function wave, a more realistic change in 
electric current can be determined. 

By way of example, there will now be described an 
EMI analysis method comprising the steps of: 

calculating the height and width of a Gauss function 
10 wave through use of a global table representing the 
relationship between an instantaneously-changing 
electric current and a Gauss function serving as a model, 
and another table representing the relationship between 
the instantaneously-changing current and the duration of 
15 the Gauss function wave; 

modeling the instantaneously-changing electric 
current as a Gauss function wave, to thereby calculate a 
more realistic change in electric current; and 

subjecting the result of calculation of an electric 
20 waveform to FFT processing, thereby determining the 

frequency characteristic of an EMI component of a circuit 
of interest. 

FIG. 50 shows an instantaneous electric current and 
a modeled Gauss function according to the present 
25 embodiment. In the drawing, ^"^i" denotes an 
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instantaneously-changing current; "t" denotes a current 
processing time; and ^^h" denotes the height of a Gauss 
function wave serving as a model and calculated from a 
global table representing the relationship between an 
5 instantaneously-changing current and the height of a Gauss 
function serving as a model . 

FIG. 51 is a graph showing an example global i-h table 
representing the relationship between an 

instantaneously-changing current value and the height of 

10 a Gauss function waveform. The table shows the value of 
^^h" at each of certain points ^^i." The table indicates 
that ^^h" assumes a value of 0 at i=0, a value of 10 at i=50, 
and a value of 40 at i=100. In reality, the value of '^i" 
changes continuously, and hence the table is expressed as 

15 a graph plotting points provided in the table. The value 
of ^^h" [i.e., h(i)] at a certain value of "i"^ is expressed 
by Eq.2 provided below, through linear interpolation of 
the values of two points closest to "i" provided in the 
i-h table: one point has lower values (il, hi), and the 

20 other point has higher values (12, h2) . Here, ^^w" 

designates the width of a Gauss function serving as a model , 
and the width "w" of a Gauss function is calculated through 
use of a global table representing the relationship 
between an instantaneously-changing current and the width 

25 of a Gauss function serving as a model . 
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FIG. 52 is a graph showing an example global i-w table 
representing the relationship between an 

instantaneously-changing current value and the width "w" 
of a rectangular waveform. The table shows the value of 
5 "w" at each of certain points "i . " The table indicates 
that ''^w" assumes a value of 0 at i=0, a value of 1 at x=20 , 
and a value of 1.5 at i=100. In reality, the value of ''"^i" 
changes continuously, and hence the table is expressed as 
a graph plotting points provided in the table. 
10 The value of ^'w" [i.e. , w(i) ] at a certain value of 

^^i" is expressed by Eq.lO provided below, through linear 
interpolation of the values of two points closest to ''"i" 
provided in the i-w table: one point has lower values (il, 
wl) , and the other point has higher values (i2, w2) . 

15 w{i)= '^^~^^ (i-n) + w\ ... (10) 

z2 - il 

In the present embodiment, the Gauss function 
serving as a model is expressed by Eq.lO. This function 
is called a Gaussian function. Let x=0 , Eq.lO provides 
"h," which represents the height of a Gaussian function 
20 wave. From Eq.lO, let x=3xW, there is derived a value of 
O.OOOlxh (i.e., 0.01% of the height of a Gaussian function 
wave) . This value is taken as a limit on the width of a 
Gaussian function serving as a model, and the width of the 
Gaussian function is taken as ^"*6xw." 
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The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
5 embodiment (see FIG. 1). 

FIG. 53 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
10 calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
15 conventional logic simulation. Reference numeral 12 

designates a calculation result of an instantaneously- 
changing current {hereinafter referred to as an 
^Unstantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
20 calculation section 11. 

FIG. 3 shows an example instantaneously-changing 
current calculation result 12. As can be seen from FIG. 
3, according to the instantaneously-changing current 
calculation result, the waveform of an electric current 
25 assumes a 6 function. Therefore, an EMI component cannot 
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be analyzed through an FFT operation. Thus far, the 
processing has been identical with the current wave 
calculation operation of the first embodiment. Reference 
numeral 4 0 designates a global i-h table representing the 
5 relationship between an instantaneously-changing current 
value and the height of a Gaussian function wave (see FIG. 
51) . Reference numeral 120 designates a Gaussian 
function waveform height calculation section for 
calculating, from an instantaneously-changing current 
10 value and the table i-h 40, a peak value of a Gaussian 
function serving as a model . 

Reference numeral 121 designates a global i-w table 
representing the relationship between an 

instantaneously-changing current value and the width of 
15 a Gaussian function wave (see FIG. 52) . Reference numeral 
122 designates a Gaussian function waveform width 
calculation section for calculating the width of a 
Gaussian function wave from the instantaneously-changing 
current value and the i-w table 121. Reference numeral 
20 123 denotes a Gaussian function waveform-shaping section 
for modeling the instantaneously-changing current result 
12 as a Gaussian function wave having a peak value 
calculated by the Gaussian function waveform height 
calculation section 120 and a width calculated by the 
25 Gaussian function waveform width calculation section 122 . 
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FIG. 54 shows the flowchart of processing pertaining 
to the Gaussian function waveform-shaping section 123, 
including processing pertaining to the Gaussian function 
5 waveform height calculation section 120 and processing 
pertaining to the Gaussian function waveform width 
calculation section 122, First, the i-h table 40 is read 
(step 1120) , and the i-w table 121 is read (step 1121) . 
The processing provided below is iterated on a per-event 

10 basis until calculation of a current waveform is completed 
(step 1122) . There are extracted two points closest to 
"*^i" provided in the i-h table 40 : one point has lower values 
(il, hi) , and the other point has higher values (12, h2) 
(step 1123) . Here, reference symbol "i" denotes the value 

15 of an instantaneously-changing current of a cell of an 
event, the event being an object of EMI analysis. 

In connection with the table shown in FIG. 51, in 
the case of i=50, there are obtained (il, hl) = (0, 0) and 
(12 , h2) = (50 , 5) . In the case of i=100 , there are obtained 

20 (il, hl) = (50, 5) and (12, h2) = (100, 20) . FromEq. 2, there 
is calculated the height '^^h" of the Gaussian function wave 
(step 1124) . In connection with the table shown in FIG. 
51, in the case of 1=50, there is obtained h(50)=10. In 
the case of 1=100, there is obtained h(100)=40. This 

25 processing corresponding to processing pertaining to the 
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Gaussian function waveform height calculation section 



There are extracted two points closest to "^^i" 
provided in the i-w table 121 : one point has lower values 
5 (il , wl) , and the other point has higher values (i2, w2) 
(step 1125) . Here, reference symbol ^^i" denotes the value 
of an instantaneously-changing current of a cell of an 
event, the event being an object of EMI analysis. In 
connection with the table shown in FIG. 52, in the case 

10 of i=50, there are obtained (il, wl)=(0, 0) and (i2, 

w2) = (20 , 1) . In the case of i=100 , there are obtained (il , 
wl)=(20, 1) and (i2 , w2)=(100, 1.5). 

Eq. 10 provides the width "^^w" of the Gaussian function 
wave (step 1126) . In connection with the table shown in 

15 FIG. 52, in the case of i=5 0 , there is obtained w (50) =1 . 3 . 
In the case of i = 100, there is obtained w (100) =1.5. The 
processing pertaining to steps 1125 and 1126 corresponds 
to processing pertaining to the Gaussian function waveform 
width calculation section 122 . 



A function expressed by Eq. 11 provided below is added 
to I(t+x) until variable ^^x" changes from -3xW to 3xW , and 
At is added to "x" (steps 1127 and 1128) . 



120. 



20 



[0158] 



y-hxe 




(11) . 
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Here, I(t+x) represents total current flowing through all 
the cells at time «t+x." This processing corresponds to 
processing pertaining to the Gaussian function 
waveform- shaping section 123. 
5 [0159] 

The current wave calculation result 8 can be obtained 
by means of the flow of a Gaussian function waveform height 
calculation operation, a Gaussian function waveform width 
calculation operation, and a Gaussian function 
10 waveform-shaping operation. FIG. 55 shows an example 
current waveform calculation result 8 as the result of 
calculation of an instantaneously-changing current shown 
in FIG. 3, the calculation result being obtained in a case 
where the table shown in FIG. 51 is used as the i-h table 
15 40 and the table shown in FIG. 52 is used as the i-w table 
121. The current wave calculation result 8 is subjected 
to an FFT operation by the FFT processing section 9, to 
thereby determine the frequency characteristic of an EMI 
component. Thus, there can be obtained an EMI analysis 
20 result 10. 

(Thirteenth Embodiment) 

An electromagnetic interference (EMI) analysis 
method according to a thirteenth embodiment of the present 
invention will now be described. The EMI analysis method 
25 according to the present embodiment is characterized by 
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modeling an instantaneously-changing current while 
information concerning estimation of a current waveform 
(hereinafter referred to simply as "current waveform 
estimation information") is separated into an element 
5 dependent on a short circuit current (also called "short 
circui t-current-dependent element") and another element 
dependent on a charge current (also called "charge- 
current-dependent element") . The expression "short 
circuit current" used herein means an electric current 

10 which flows from the source terminal to the drain terminal 
in each of a p-channel MOS transistor (hereinafter 
abbreviated "p-MOS transistor") and an n-channel MOS 
transistor (hereinafter abbreviated "n-MOS transistor") 
in a case where the p-MOS and n-MOS transistors are 

15 instantaneously and simultaneously brought into an active 
state as a result of a change arising in a signal appearing 
at the gate terminal of a CMOS transistor constituting a 
circuit element. The expression "charge current" means 
an electric current which flows when electric charge is 

20 stored in or discharged from an output load capacitor of 
a circuit element. A short circuit current arises 
substantially simultaneous with occurrence of a change in 
an input signal, and a charge current flows together with 
a change in an output signal at a slight time lag from the 

25 time when change has arisen in an input signal. 
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Accordingly, estimation of a more accurate current 
waveform becomes feasible by means of separate handling 
of the short circuit-current-dependent element of the 
current waveform estimation information and the 
5 charge-current-dependent element. 

FIG. 56 shows the configuration of the EMI analysis 
method according to the present embodiment using a 
technique of averaging a short circuit current and a charge 
current by means of dividing by a discrete width. 

10 A system used with the EMI analysis method shown in 

FIG. 56 has means for analyzing the amount of current 
(hereinafter referred to as ^^amount-of -current analysis 
means") , the amount-of -current analysis means comprising 
circuit-connection information storage means 7301 for 

15 storing information concerning circui t interconnections / 
pattern-of -signal-change storage means 7302 for storing 
a pattern of signal change; element-current storage means 
7303 for storing the amount of electric current flowing 
through an element; basi c-inf ormation-f or-es timating- 

20 current storage means 7306 for storing basic information 
concerning estimation of an electric current; total- 
amount-of -current storage means 7304 for storing the total 
amount of current flowing through an element; and 
total-amount-of-current calculation means 7305 for 

25 calculating the total amount of electric current. 
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These means constituting the EMI analysis method 
shown in FIG. 56 will now be described. There will now 
be described procedures for calculating the total amount 
of electric current shown in FIG. 58 from circuit 
5 connection information shown in FIG. 110, a pattern of 
signal change shown in FIG. Ill, and basic information 
concerning estimation of an electric current shown in FIG. 
59. The amount-of-current analysis means comprising 
circuit-connection information storage means 7301 is 

10 identical with the circuit connection information storage 
means 1001 described in connection with the background art, 
and the pattern-of- signal -change storage means 7302 is 
identical with the signal pattern storage means 1002 
described in connection with the background art- Hence, 

15 their repeated explanations are omitted for brevity. 

The element-current storage means 7303 stores 
information concerning the amount of electric current 
flowing through an element; that is, information 
concerning the total amount of electric current flowing 

20 when single change arises in an external terminal of each 
element, while the information is divided into a short 
circuit current component and a charge current component. 
The information concerning the amount of electric current 
flowing through an element, such as that shown in FIG. 57, 

25 has been stored in the element-current storage means 7303 
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beforehand. Reference numeral 74 01 denotes the amount of 
short circuit current which flows when a change arises in 
a signal appearing at a clock input terminal CK of a 
flip-flop FF. Reference numeral 7402 denotes the amount 
5 of charge current which flows when a change arises in a 
signal appearing at a clock input terminal CK of the 
flip-flop FF. Reference numeral 7403 denotes the amount 
of short circuit current which flows when a change arises 
in a signal appearing at an output terminal Q of the 

10 flip-flop FF. Reference numeral 7404 denotes the amount 
of charge current which flows when a change arises in a 
signal appearing at the output terminal Q of the flip- 
flop FF . Reference numeral 74 05 denotes the amount of 
short circuit current which flows when a change arises in 

15 a signal appearing at an output terminal Y of a buffer BUF . 
Reference numeral 74 06 denotes the amount of charge 
current which flows when a change arises in a signal 
appearing at the output terminal Y of the buffer BUF. The 
basic-inf ormation-f or-estimating-current storage means 

20 7306 stores the discrete width of a short circuit current 
and the discrete width of a charge current in the form shown 
in FIG. 59. The basic information concerning estimation 
of an electric current has a discrete width to be used for 
averaging the total amount of electric current which flows 

25 in response to occurrence of a change in a signal, as in 
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the case of the information employed in the first 
embodiment. However, in contrast with the information 
employed in the first embodiment, the basic information 
concerning estimation of an electric current comprises two 
5 pieces of information: that is, information 7601 

concerning the discrete width of a short circuit current, 
and information 7602 concerning the discrete width of a 
charge current. Information 7601 represents the discrete 
width of a short circuit current, and information 7602 

10 represents the discrete width of a charge current. 

The total -amount-of -current calculation means 7305 
operates according to a flowchart such as that shown in 
FIG. 60. In contrast with the conventional method, the 
flowchart shown in FIG. 60 of the EMI analysis method of 

15 the present embodiment is additionally provided with a 
processing operation for averaging the calculated short 
circuit current by means of dividing by the discrete width 
thereof and averaging the calculated charge current by 
means of dividing by the discrete width thereof (step 7706) . 

20 In step 7701, there is read the circuit connection 

information, which is shown in FIG. 110 and is stored in 
the circuit-connection information storage means 7301. 
In step 7702, there is read the pattern of a change in a 
signal, the pattern being stored in the pattern-of- 

25 signal-change storage means 7302. In step 7703, there is 
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read information concerning the amount of electric current 
flowing through an element, the information being stored 
in the element-current storage means 7303 and shown in FIG. 
57. 

5 In step 7704, there is read the information 

concerning estimation of a current, the information being 
stored in the basic-information-for-estima ting-current 
storage means 7306 and shown in FIG. 59. In step 7705, 
the pattern of a change in a signal shown in FIG. Ill is 

10 imparted to the circuit connection information shown in 
FIG. 110, to thereby enable propagation of a signal. In 
a case where a change arises in a signal appearing at an 
external terminal of a circuit element, the element being 
described in the information concerning the amount of 

15 electric current flowing through an element, a short 

circuit current component and a charge current component 
included in the electric current flowing through the 
external terminal are added, at respective points in time, 
to the short circuit current component and the charge 

20 current component included in the information concerning 
the total amount of electric current (simply called as 
"total current information"), respectively. 
[0166] 

In step 7706, the result of an electric current being 
25 averaged by means of dividing by the discrete width stored 
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in the information concerning estimation of a current 
shown in FIG. 59 is added to the short circuit current 
component included in the total current information. The 
result of an electric current being averaged by means of 
5 dividing by the discrete width stored in the information 
concerning estimation of a current shown in FIG. 59 is 
added to the charge current component included in the total 
current information. In step 7707, after propagation of 
a signal, information concerning the total amount of 

10 electric current calculated at each points in time shown 
in FIG. 58 is stored in the total-amount-of -current 
storage means 7304. As a result, there is obtained a 
current waveform, such as that shown in FIG. 58, which is 
more realistic than the current waveform obtained by the 

15 conventional method shown in FIG. 14. 
[0167] 

(Fourteenth Embodiment) 

An EMI analysis method according to a fourteenth 
embodiment of the present invention will now be described. 
20 The present embodiment describes an EMI analysis method 
comprising the steps of: 

characterizing a peak current of a cell according 
to the type of cell and the amount of instantaneously- 
changing current of each event, to thereby calculate an 
25 instantaneously-changing current of each cell from event 
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information, the information being output from a logic 
simulator ; 

modeling the instantaneously-changing current as a 
rectangular waveform whose width is calculated such that 
5 a peak value corresponds to the peak current of the 
rectangular waveform and such that the 

instantaneously-changing current becomes equal to the 
area of the rectangular waveform, to thereby calculate a 
more realistic change in an electric current; and 

10 subjecting the result of calculation of a current 

waveform to FFT processing, to thereby determine the 
frequency characteristic of an EMI component arising in 
a circuit of interest. 

FIG. 15 shows an instantaneously changing current 

15 and a modeled rectangular waveform according to the 
present embodiment. Here, ^"^i" denotes an 
instantaneously-changing current, and "t" denotes a 
current processing time. 

The overall flow of processing of the EMI analysis 

20 method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1). 

FIG. 61 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 

25 present embodiment. Reference numeral 11 designates an 
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instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
5 Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 

10 ^^instantaneously-changing current calculation result 

12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 

15 instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 

20 the first embodiment. 

Reference numeral 130 denotes a peak current library. 
A peak current of each cell is characterized according to 
the type of cell, and the thus-characterized value is 
stored in the peak current library 130 . Reference numeral 

25 131 designates a sixth rectangular waveform height 
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calculation section for extracting, from the peak current 
library 130 , a peak value of a rectangular waveform serving 
as a model. Reference numeral 132 designates a sixth 
rectangular waveform width calculation section for 
5 calculating the width of a rectangular waveform, from the 
peak value extracted by the sixth rectangular waveform 
height calculation section 131 and the 

instantaneously-changing current value, such that the 
instantaneously-changing current becomes equal to the 

10 area of the rectangular waveform. Reference numeral 133 
denotes a sixth rectangular waveform-shaping section for 
modeling the instantaneously-changing current 
calculation result 12 as a rectangular waveform, the 
height of the rectangular waveform corresponding to the 

15 peak value extracted by the sixth rectangular waveform 
height calculation section 131, and the width of the 
rectangular waveform corresponding to the width 
calculated by the sixth rectangular waveform width 
calculation section 132. 

20 FIG. 62 shows the flowchart of processing pertaining 

to the sixth rectangular waveform-shaping section 133, 
including processing pertaining to the sixth rectangular 
waveform height calculation section 131 and processing 
pertaining to the sixth rectangular waveform width 

25 calculation section 132 . First, the peak current library 
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130 is read (step 1130) . The following processing is 
iterated on a per-event basis until calculation of a 
current waveform is completed (step 1130) . Apeak current 
of an instance which is an object of EMI analysis is 
5 extracted from the peak current library 130, and the peak 
current is taken as the height "h'' of a rectangular 
waveform (step 1132). This processing corresponds to 
processing pertaining to the sixth rectangular waveform 
height calculation section 131. 

10 The area of the rectangular waveform is defined by 

"w X h." The width "w" of a rectangular waveform is 
calculated on the basis of "w = i/h" (step 1133) . 
Reference symbol ^^i" denotes an instantaneously-changing 
current of a cell of an instance, the instance being an 

15 object of EMI analysis. Reference symbol "h" denotes the 
height ^^h" of a rectangular waveform calculated through 
processing pertaining to the sixth rectangular waveform 
height calculation section 131. This processing 
corresponds to processing pertaining to the sixth 

20 rectangular waveform width calculation section 132 . Next, 
"h" is added to I(t+x) until a variable ^^x" changes from 
-W/2 to W/2, and At is added to ^^x" (steps 1134 and 1135) . 
I(t+x) corresponds to total electric current flowing 
through all the cells at time t+x . This processing 

25 corresponds to processing pertaining to the sixth 
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rectangular waveform- shaping section 133. 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 
calculation operation, a rectangular waveform width 
5 calculation operation, and a rectangular waveform- shaping 
operation. The FFT processing section 9 subjects the 
current waveform calculation result 8 to an FFT operation, 
to thereby determine the frequency characteristic of an 
EMI component. As a result, there can be obtained an EMI 
10 analysis result 10. 

(Fifteenth Embodiment) 

An EMI analysis method according to a fifteenth 
embodiment of the present invention will now be described. 
The present embodiment describes an EMI analysis method 
15 comprising the steps of: 

characterizing the relationship between a slew in 
a waveform input to each cell and a peak current of the 
cell, according to the type of cell, and the amount of 
instantaneously-changing current for each event, to 
20 thereby calculate an instantaneously-changing current of 
each cell from event information concerning each cell, the 
information being output from a logic simulator; 

modeling the instantaneously-changing current as a 
rectangular waveform whose width is calculated such that 
25 a peak current taking into consideration a slew in the 
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waveform input to the cell corresponds to a peak current 
of the waveform and such that the instantaneously-changing 
current becomes equal to the area of the rectangular 
waveform, to thereby calculate a more realistic change in 
5 an electric current; and 

subjecting the result of calculation of a current 
waveform to FFT processing, to thereby determine the 
frequency characteristic of an EMI component arising in 
a circuit of interest. 

10 FIG. 15 shows an instantaneously changing current 

and a modeled rectangular waveform according to the 
present enibodiment. Here, ^^i" denotes an 
instantaneously-changing current, and "t" denotes a 
current processing time. Further, "h" denotes the height 

15 of a rectangular waveform. The height "h" of a rectangular 
waveform is calculated from a library which characterizes 
the relationship between a slew in the waveform input to 
each cell and the height of a rectangular waveform. 
Reference symbol ^^w" denotes the width of a rectangular 

20 waveform, and the width ^^w" of a rectangular waveform is 
calculated from the height ^^h" calculated previously and 
the instantaneously-changing current "i." The library 
characterizing the relationship between a slew in an input 
waveform and the height of a rectangular waveform is 

25 prepared according to the type of cell in the form of an 
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s-h table. 

FIG. 63 shows an example s-h table and a graph 
plotting the table. Here, ^^s" denotes a slew in a waveform 
input to an instance which is an object of EMI analysis. 
5 The table shows the value of "h" at each of certain points 
^^s." The table indicates that ^^h" assumes a value of 0 
at s=0 , a value of 5 at s=l , and a value of 2 0 at s=2 . In 
reality, the value of ^^s" changes continuously, and hence 
the table is expressed as a graph plotting points provided 
10 in the table. The value of ^^h" [i.e., h(s)] at a certain 
value of ^^s" is expressed by Eq . 12 provided below, through 
linear interpolation of the values of two points closest 
to ^^s" provided in the s-h table : one point has lower values 
(si, hi) , and the other point has higher values (s2, h2) . 

15 h(s) = ^^~^^ (s-sl) + hl ... (12) 

52-51 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 
20 FIG. 64 is a detailed block diagram showing the 

current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
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from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
5 conventional logic simulation. Reference numeral 12 

designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
^^instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 

10 calculation section 11. FIG. 3 shows an example 

instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 

15 Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. 

In the present embodiment, an s-h table library 14 0 

20 characterizes the relationship between a slew in a 

waveform input to each cell and the height of a rectangular 
waveform (see FIG. 63) . Reference numeral 51 designates 
Slew data describing a slew in a waveform input to each 
instance. Reference numeral 141 designates a seventh 

25 rectangular waveform height calculation section for 
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extracting a slew in a waveform input to each instance, 
the instance being an object of EMI analysis, from the Slew 
data 51, and calculating a peak value of a rectangular 
waveform serving as a model, from the s-h table library 
5 140 and the thus -extracted slew. Reference numeral 142 
designates a seventh rectangular waveform width 
calculation section for calculating the width of a 
rectangular waveform from the peak value calculated by the 
seventh rectangular waveform height calculation section 

10 141 and from the instantaneously-changing current, such 
that the area of the rectangular waveform becomes equal 
to the instantaneously-changing electric current. A 
seventh rectangular waveform- shaping section 14 3 models 
the instantaneously-changing current calculation result 

15 12 as a rectangular waveform, the height of the rectangular 
waveform corresponding to the peak value of the 
rectangular waveform calculated by the seventh 
rectangular waveform height calculation section 141, and 
the width of the rectangular waveform corresponding to the 

20 width calculated by the seventh rectangular waveform width 
calculation section 142. 

FIG. 65 is a flowchart showing processing pertaining 
to the seventh rectangular waveform-shaping section 143, 
including processing pertaining to the seventh 

25 rectangular waveform height calculation section 141 and 
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processing pertaining to the seventh rectangular waveform 
width calculation section 142. First, the s-h table 
library 140 is read (step 1140) . The Slew data 51 are read 
{step 1141) . Next, the subsequent processing operations 
5 are iterated on a per-event basis until the current 

waveform calculation operation is completed (step 1142) . 
Two points closest to ^^s"; (si, hi) and (s2, h2), are 
extracted from the s-h library table 140. Reference 
symbol ^^s" denotes a slew in an input waveform, the slew 
10 being extracted from the Slew data 51 concerning an 
instance which is an object of EMI analysis. 

In connection with the table shown in FIG. 63, in 
the case where s=l . 2 , the two points closest to s=l . 2 are 
(si, hi) = (1, 5) and (s2, h2) = (2, 20) (step 1143) . The 
15 height "h" of a rectangular waveform is calculated 

accordingtoEq.il (step 1144). In the table shown in FIG . 
63, in the case where s=l . 2 , h(1.2)=8. This processing; 
i.e., processing pertaining to steps 1143 and 1144, 
corresponds to processing pertaining to the seventh 
20 rectangular waveform height calculation section 141. 

The area of a rectangular waveform is defined as "'w 
X h." The width "w" of a rectangular waveform is 
calculated from w=i/h (step 1145) . Reference symbol ^^i" 
denotes the value of an instantaneously-changing current 
25 of a cell of an event, the event being an object of EMI 
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analysis. Reference symbol "h" denotes the height of a 
rectangular waveform calculated by the seventh 
rectangular waveform height calculation section 141. 
This processing corresponds to processing pertaining to 
5 the seventh rectangular waveform width calculation 

section 142. Next, ^^h" is added to I (t+x) until variable 
"x" changes from -W/2 to W/2 , and At is added to variable 
"X" (steps 1146 and 1147) . I (t+x) denotes total electric 
current flowing through all the cells at time 
10 "t+x." This processing corresponds to processing 

pertaining to the seventh rectangular waveform-shaping 
section 143 . 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 

15 calculation operation, a rectangular waveform width 

calculation operation, and a rectangular waveform-shaping 
operation. The FFT processing section 9 subjects the 
current waveform calculation result 8 to an FFT operation, 
to thereby determine the frequency characteristic of an 

20 EMI component. As a result, there can be obtained an EMI 
analysis result 10. 
(Sixteenth Embodiment) 

An EMI analysis method according to a sixteenth 
embodiment of the present invention will now be described. 

25 The present embodiment describes an EMI analysis method 
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comprising the steps of: 

characterizing the relationship between an output 
load capacitance of each cell and a peak current of the 
cell, according to the type of cell, and the amount of 
5 instantaneously-changing current for each event, to 

thereby calculate an instantaneously-changing current of 
each cell from event information concerning each cell, the 
information being output from a logic simulator; 

modeling the instantaneously-changing current as a 

10 rectangular waveform whose width is calculated such that 
a peak current taking into consideration of the output load 
capacitance of the cell corresponds to a peak current of 
the waveform and such that the instantaneously-changing 
current becomes equal to the area of the rectangular 

15 waveform, to thereby calculate a more realistic change in 
an electric current; and 

subjecting the result of calculation of a current 
waveform to FFT processing, to thereby determine the 
frequency characteristic of an EMI component arising in 

20 a circuit of interest. 

FIG. 15 shows an instantaneously changing current 
and a modeled rectangular waveform according to the 
present embodiment. Here, "i" denotes an 
instantaneously-changing current, and ^'*t" denotes a 

25 current processing time. Further, '^h" denotes the height 
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of a rectangular waveform. The height "*'h" of a rectangular 
waveform is calculated from a library which characterizes 
the relationship between the output load capacitance of 
each cell and the height of a rectangular waveform. 
5 Reference symbol ^^w" denotes the width of a rectangular 
waveform, and the width ^^w" of a rectangular waveform is 
calculated from the height ''^h" calculated previously and 
the instantaneously-changing current ^^i . " The library 
characterizing the relationship between the output load 
10 capacitance and the height of a rectangular waveform is 
prepared according to the type of cell in the form of a 
c-h table. 

FIG. 6 6 shows an example c-h table and a graph 
plotting the table. Here, "s" denotes an output load 

15 capacitance of an instance which is an object of EMI 

analysis. The table shows the value of "*^h" at each of 
certain points ^^c." The table indicates that "h" assumes 
a value of 0 at c=0 , a value of 5 at c=10, and a value of 
20 at c=20. In reality, the value of ^^c" changes 

20 continuously, and hence the table is expressed as a graph 
plotting points provided in the table. The value of ^^h" 
[i.e., h (c) ] at a certain value of "c" is expressed by Eq . 13 
provided below, through linear interpolation of the values 
of two points closest to "c" provided in the s-h table: 

25 one point has lower values {cl, hi) , and the other point 
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has higher values (c2, h2) . 

/j(c) = :^^Il^(c-cl) + M ... (13) 
c2-cl 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
5 by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 67 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 

10 instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 

15 analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously-changing current calculation result 

20 12") determined by the instantaneously-changing current 
calculation section 11. 

FIG. 3 shows an example instantaneously-changing 
current calculation result 12. As can be seen from FIG. 
3, according to the instantaneously-changing current 
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calculation result, the waveform of an electric current 
assumes a 5 function. Therefore, an EMI component cannot 
be analyzed through an FFT operation. Thus far, the 
processing has been identical with the current wave 
5 calculation operation of the first embodiment . Reference 
numeral 150 designates a table library. The c-h table 
library 150 is made up of a c-h table library (see FIG. 
66) characterizing the relationship between the 
previously-described output load capacitance of each cell 

10 and the height of a rectangular waveform. Reference 
numeral 61 designates capacitance data describing an 
output load capacitance of each instance. Reference 
numeral 151 designates an eighth rectangular waveform 
height calculation section for extracting an output load 

15 capacitance of an instance, the instance being an object 
of EMI analysis, from the capacitance data 61 and 
calculating a peak value of a rectangular waveform serving 
as a model, from the c-h table library 150 and the 
thus -extracted output load capacitance. Reference 

20 numeral 152 designates an eighth rectangular waveform 
width calculation section for calculating the width of a 
rectangular waveform from the peak value calculated by the 
eighth rectangular waveform height calculation section 
151 and from the instantaneously-changing current, such 

25 that the area of the rectangular waveform becomes equal 
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to the instantaneously-changing electric current. An 
eighth rectangular waveform- shaping section 153 models 
the instantaneously-changing current calculation result 
12 as a rectangular waveform, the height of the rectangular 
5 waveform corresponding to the peak value of the 

rectangular waveform calculated by the eighth rectangular 
waveform height calculation section 151, and the width of 
the rectangular waveform corresponding to the width 
calculated by the eighth rectangular waveform width 

10 calculation section 152. 

FIG. 68 is a flowchart showing processing pertaining 
to the eighth rectangular waveform-shaping section 153, 
including processing pertaining to the eighth rectangular 
waveform height calculation section 151 and processing 

15 pertaining to the eighth rectangular waveform width 

calculation section 152. First, the c-h table library 150 
is read (step 1150) . The capacitance data 61 are read 
(step 1151) . 

Next, the subsequent processing operations are 
20 iterated on a per-event basis until the current waveform 
calculation operation is completed (step 1152) . Two 
points closest to "c"; (cl, hi) and (c2, h2) , are extracted 
from the c-h library table 150. Reference symbol "c" 
denotes a slew in an input waveform, the slew being 
25 extracted from the capacitance data 61 concerning an 
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instance which is an object of EMI analysis. 

In connection with the table shown in FIG. 66, in 
the case where c=12, the two points closest to c=12 are 
(cl, hi) = (1, 5) and (c2, h2) = (20, 20) (step 1153) . The 
5 height ^^h" of a rectangular waveform is calculated 

according to Eq . 12 (step 1154). In the table shown in FIG . 
66 , in the case where c=12 , h (12 ) =8 . This processing; i.e., 
processing pertaining to steps 1153 and 1154, corresponds 
to processing pertaining to the eighth rectangular 

10 waveform height calculation section 151. 

This processing corresponds to processing 
pertaining to the eighth rectangular waveform width 
calculation section 152, Next, "h" is added to I(t+x) 
until variable ^^x" changes from -W/2 to W/2 , and At is added 

15 to variable '''x" (steps 1156 and 1157). I (t+x) denotes 
total electric current flowing through all the cells at 
time 

"t+x." This processing corresponds to processing 
pertaining to the eighth rectangular waveform-shaping 

20 section 153. 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 
calculation operation, a rectangular waveform width 
calculation operation, and a rectangular waveform-shaping 

25 operation. The FFT processing section 9 subjects the 
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current waveform calculation result 8 to an FFT operation, 
to thereby determine the frequency characteristic of an 
EMI component. As a result, there can be obtained an EMI 
analysis result 10. 
5 (Seventeenth Embodiment) 

An electromagnetic interference (EMI) analysis 
method according to a seventeenth embodiment of the 
present invention will now be described. Among 
operations related to EMI analysis, the present embodiment 

10 describes particularly extraction of a characteristic of 
EMI and calculation of the amount of an electric current, 
in the case where a circuit element is a composite cell 
formed by combination of a plurality of logic elements. 
In a composite cell, an electric current flows every time 

15 each of logic elements constituting the circuit element 
is activated. A logic element provided in a circuit 
element is activated with a time lag unique to the logic 
element . 

Specifically, when a change arises in an input to 
20 a circuit element, an electric current flows many times 
within the circuit element. In a circuit element 
including logic elements in combination as shown in FIG. 
108, INl is changed with a lag of a gate PI provided in 
a cell (provided that the lag is Ins) from the time changes 
25 have arisen in terminals A and B. Further, Y is changed 
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with a lag of a gate P2 provided in the cell (provided that 
the lag is Ins) from the time a change has arisen in INl . 
In a case where a pattern of signal change shown in FIG. 
69 is imparted to the circuit element (cell) , there arises 
5 a current waveform such as that shown in FIG. 70. 

According to the conventional gate-level EMI 
analysis method, the total amount of electric current 
flowing through each circuit element is considered to be 
dissipated with an involved time lag until a signal appears 

10 at an output pin of the circuit element. For instance, 
in a case where a pattern of signal change shown in FIG. 
69 is imparted to a circuit shown in FIG. 73, the circuit 
employing a circuit element shown in FIG. 108, there is 
output current information such as that shown in FIG. 71. 

15 In contrast, according to the EMI analysis method of the 
present embodiment, an electric current, which would arise 
when a change occurs in a signal appearing at an input pin 
of a composite cell serving as an object, is registered 
in a library, through use of a peak value of the electric 

20 current and a time lag based on a change in an input signal . 
A logic simulator estimates an electric current, on the 
basis of the assumption that a peak value of an electric 
current includes the previously-described time lag in 
accordance with a change in a signal appearing at an input 

25 pin of a composite cell. 
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FIG. 72 shows the configuration of an EMI analysis 
method using the value of a peak current and the amount 
of electric current. 

A system to be used with the EMI analysis method shown 
5 in FIG. 72 comprises circuit connection information 
storage means 9001; pattern-of -signal-change storage 
means 9002; element current storage means 9003 for storing 
the amount of electric current flowing in a circuit 
element; total current storage means 9004 for storing a 

10 total amount of electric current flowing through the 

circuit element; and total current calculation means 9005 
for calculating the total amount of electric current 
flowing through a circuit element. 

Next will be described individual elements 

15 constituting the EMI analysis method shown in FIG. 72 and 
procedures for calculating the total amount of electric 
current shown in FIG. 74, from the circuit connection 
information shown in FIG. 73, the pattern of signal change 
shown in FIG. 69, and the amount of electric current 

20 flowing through a circuit element shown in FIG. 74. 

The circuit connection information storage means 
9001 stores information concerning the connection of a 
circuit to be analyzed, in the same manner as described 
in connection with the background art. Circuit 

25 connection information, such as that shown in FIG. 73, is 
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stored in advance in the circuit connection information 
storage means 9001. The pattern-of -signal-change 
storage 9002 stores the pattern of signal change to be 
applied to the circuit connection information, as in the 
5 case of a conventional pattern of signal change described 
in connection with the background art. As shown in FIG. 
69, a pattern of change in a signal input to an input port 
CK and a pattern of change in a signal input to an input 
port D, such as those shown in FIG. 73, are stored in the 

10 pattern-of -signal-change storage 9002 beforehand. 

The element current storage means 9003 stores, for 
each circuit element and in a form such as that shown in 
FIG. 75, pins 9301 in which electric current arises in an 
element after a change has arisen in an external input 

15 terminal of each element, time intervals 9302 (i.e., the 
time which elapses after a change has arisen in an external 
input terminal), total electric current 9303, and peak 
values 9304 . Reference numeral 9305 designates flow of 
an electric current in an ANDOR circuit shown in FIG. 108; 

20 specifically, a total amount of 1 mA including a peak 
current of 0.5 mA flows into the circuit element after 
lapse of 1 ns from when input A has arisen. After lapse 
of 2 ns from when the input A has changed, a total mount 
of 3 mA including a peak value of 1.5 mH flows into the 

25 circuit element. Reference numeral 9306 designates flow 
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of an electric current in a circuit element shown in FIG. 
108; namely, a total amount of 1 mA including apeak current 
of 0.5 xnA flows after lapse of 1 ns from when input B has 
changed. Further, after lapse of 2 ns from when the input 
5 B has changed, an electric current of 3 ma including a peak 
value of 1.5 ma flows into the circuit. Reference numeral 
9307 designates flow of an electric current in a circuit 
element shown in FIG. 108; specifically, a total amount 
of 3 mA including a peak value of 1.5 mA flows into the 
10 circuit element after lapse of 2 ns from when input C has 
changed . 

The total current calculation means 9005 operates 
according to a flowchart such as that shown in FIG. 76. 
According to the flowchart shown in FIG. 76 concerning a 

15 method of calculating an electric current, processing 
pertaining to step 1504 of the conventional method of 
calculating an electric current is replaced by processing 
(pertaining to step 9404) in which an electric current is 
calculated while a rectangular waveform of the electric 

20 current is assumed on the basis of the total amount of 
electric current and a peak value with a time interval 
provided between information pieces concerning the amount 
of current of an element, the information having been 
stored in the element current storage means in association 

25 with a change in an input signal. In step 9401 is read 



138 




circuit connection information shown in FIG. 73 stored in 
the circuit connection information storage means 9001. 
In step 9402 is read a pattern of signal change which is 
shown in FIG. 69 and is stored in the pattern-of- 
5 change-in-signal storage means 9002. In step 9403 is read 
element current information, which is shown in FIG. 75 and 
is stored in the element current information storage means 
9003. 

In step 9404, a pattern of signal change shown in 

10 FIG. 69 is imparted to the circuit connection information 
shown in FIG. 73 in step 9404 , thereby enabling propagation 
of a signal. In a case where a change arises in a signal 
appearing at an external input terminal of a circuit 
element, which element is described in the information 

15 concerning the amount of electric current flowing an 

element, a rectangular waveform of an electric current is 
produced, by means of averaging the amount of electric 
current flowing through the external input terminal by 
means of a peak value. The thus -produced rectangular 

20 waveform is added to the information concerning the total 
mount of electric current obtained at a time which lags 
a predetermined period of time from the time a change has 
arisen in the external input terminal. After completion 
of propagation of a signal, in step 9405 the thus- 

25 calculated information concerning the total amount of 
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electric current obtained at each time shown in FIG. 74 
is stored in the total current storage means 9004 . As a 
result, there is produced a current waveform, such as that 
shown in FIG. 74, which is more realistic than that which 
5 is obtained by the conventional current calculation method 
and shown in FIG. 71. 

(Eighteenth Embodiment) 

An electromagnetic interference analysis method 

10 according to an eighteenth embodiment of the present 
invention will now be described. 

Among operations related to EMI analysis, the 
present embodiment describes in particular extraction of 
a characteristic of EMI and calculation of the amount of 

15 an electric current, in the case where a circuit element 
is a register (signal holding) element synchronized with 
a clock signal . In a register element such as a flip- 
flop, an electric current flows at every rise or fall of 
a change which arises in a signal appearing at a clock pin. 

20 For example, a pattern of signal change shown in FIG. 78 
is input to a circuit shown in FIG. 77, an electric current 
shown in FIG. 79 flows through the circuit. As in the case 
of the seventeenth embodiment, an electric current flows 
every time each of the logic elements constituting the 

25 circuit element is activated. A conventional tool for 



140 




analyzing power consumption on a gate level represents a 
signal change only when data are input to a flip flop; for 
example, at the time of a rise in a clock signal pin. For 
example, in a case where a pattern of signal change shown 
5 in FIG. 78 is input to the circuit shown in FIG. 77, the 
circuit outputs current information such as that shown in 
FIG. 80. 

In the present embodiment, the timing at which an 
electric current flows through a flip-flop is determined 

10 by means of taking, as a reference, a rise and fall of a 
clock signal appearing at the flip-flop. Therefore, an 
electric current is estimated by means of taking, as a 
reference, a rise and fall of the clock signal. FIG. 72 
shows the configuration of an EMI analysis method using 

15 the value of a peak current and the amount of electric 
current. A system to be used with the EMI analysis method 
shown in FIG. 72 comprises circuit connection information 
storage means 9001; pattern-of -signal -change storage 
means 9002; element current storage means 9003 for storing 

20 the amount of electric current flowing in a circuit 

element; total current storage means 9004 for storing a 
total amount of electric current flowing through the 
circuit element; and total current calculation means 9005 
for calculating the total amount of electric current 

25 flowing through a circuit element. 
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[0203] 

Next will be described individual elements 
constituting the EMI analysis method shown in FIG. 72 and 
procedures for calculating the total amount of electric 
5 current shown in FIG. 82, from the circuit connection 
information shown in FIG. 77, the pattern of signal change 
shown in FIG. 78, and the amount of electric current 
flowing through a circuit element shown in FIG. 81. 

The circuit connection information storage means 
10 9001 stores information concerning the connection of a 
circuit to be analyzed, in the same manner as described 
in connection with the background art. Circuit 
connection information, such as that shown in FIG. 77, is 
stored in advance in the circuit connection information 
15 storage means 9001. 
[0205] 

The pattern-of -signal -change storage 9002 stores 
the pattern of signal change to be applied to the circuit 
connection information, as in the case of a conventional 

20 pattern of signal change described in connection with the 
background art. As shown in FIG. 78, a pattern of change 
in a signal input to an input port CK and a pattern of change 
in a signal input to an input port D, such as those shown 
in FIG. 77, are stored in the pattern-of -signal -change 

25 storage 9002 beforehand. 
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The element current storage means 9003 stores, for 
each circuit element and in a form such as that shown in 
FIG. 81, information 9901 concerning the total amount of 
electric current flowing in a circuit element at a rise 
5 and fall of a clock signal and information 9902 about a 
peak value of the electric current flowing at a rise and 
fall of the clock si.gnal . 

The total current calculation means 9005 operates 
according to a flowchart, such as that shown in FIG. 83. 

10 According to the flowchart relating to a method of 
calculating an electric current shown in FIG. 83, 
processing pertaining to step 1504 of the conventional 
current calculation method is replaced by processing of 
calculating a rectangular waveform of the electric current 

15 from a peak value and a total amount of electric current 
flowing at a rise and fall of a clock signal (step 10104) . 

In step 10101, there is read the circuit connection 
information which is shown in FIG. 77 and stored in the 
circuit connection information storage means 9001. In 

20 step 10102, there is read a pattern of signal change which 
is shown in FIG. 78 and is stored in the pattern-of- 
signal-change storage 9002. 

In step 10103, there is read the element current 
information which is shown in FIG. 114 and stored in the 

25 element current storage means 9003. In step 10104, the 
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pattern of signal shown in FIG. 78 is imparted to the 
circuit connection information shown in FIG. 77, thereby 
enabling propagation of a signal . In a case where a change 
arises in a signal appearing at an external input terminal 
5 of a circuit element, which element is described in the 
information concerning the amount of electric current 
flowing an element, a rectangular waveform of an electric 
current is produced by means of averaging an electric 
current corresponding to the electric current flowing 

10 through the external input terminal at a rise and fall of 
the signal by means of a peak value. The thus -produced 
rectangular waveforms are added to the information 
concerning the total mount of electric current obtained 
at a time when the change has arisen in the signal appearing 

15 at the external input terminal . 

After completion of propagation of a signal, in step 
10105 the thus-calculated information concerning the 
total amount of electric current obtained at each time 
shown in FIG. 82 is stored in the total current storage 

20 means 9004. As a result, there is produced a current 
waveform, such as that shown in FIG. 82, which is more 
realistic than that obtained by the conventional current 
calculation method and shown in FIG. 80. 
(Nineteenth Embodiment) 

25 An electromagnetic interference analysis method 
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according to a nineteenth embodiment of the present 
invention will now be described. 

Conventional gate-level power consumption analysis 
means can change total electric current according to the 
5 value of a signal input to a logic element. More 

specifically, the conventional power consumption analysis 
means enables computation of so-called state-dependent 
power consumption. However, the conventional means 
cannot effect state-dependent power consumption including 
10 secular change in electric current and secular change in 
peak value. 

Further, some conventional means determines a 
stationary value of another signal in the event of 
occurrence of a change in a certain signal . However, such 

15 means cannot determine a change in another signal. For 
instance, an electric current such as that shown in FIG. 
8 6 (an electric current may change in amount according to 
whether a change arises in both terminals A and B 
simultaneously or in only terminal A) , on the basis of the 

20 circuit connection information shown in FIG. 84 as well 
as a pattern of signal change shown in FIG. 85. In such 
a case, the conventional means can output only a waveform 
shown in FIG. 87 . 

In contrast, the EMI analysis method according to 

25 the present embodiment is characterized by holding secular 
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change in electric current and a peak value in the form 
of a state -dependent signal which takes into consideration 
changes in two or more signals. FIG. 72 shows the 
configuration of an embodiment of the EMI analysis method 
5 using the peak current and the amount of electric current, 
which are calculated in the fourteenth embodiment. 

A system to be used with the EMI analysis method shown 
in FIG. 72 comprises circuit connection information 
storage means 9001; pattern-of -signal -change storage 

10 9002; element current storage means 9003 for storing the 
amount of electric current flowing in a circuit element; 
total current storage means 9004 for storing a total amount 
of electric current flowing through the circuit element; 
and total current calculation means 9005 for calculating 

15 the total amount of electric current flowing through a 
circuit element. 

Next will be described individual elements 
constituting the EMI analysis method shown in FIG. 72 and 
procedures for calculating the total amount of electric 

20 current shown in FIG. 90, from the circuit connection 
information shown in FIG. 84, the pattern of signal change 
shown in FIG. 85, and the amount of electric current 
flowing through a circuit element shown in FIG. 88. 

The circuit connection information storage means 

25 9001 stores information concerning the connection of a 
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circuit to be analyzed, in the same manner as described 
in connection with the background art. Circuit 
connection information, such as that shown in FIG, 84, is 
stored in advance in the circuit connection information 
5 storage means 9001. 

The pattern-of -signal-change storage 9002 stores 
the pattern of signal change to be applied to the circuit 
connection information, as in the case of a conventional 
pattern of signal change described in connection with the 

10 background art. As shown in FIG. 85, a pattern of change 
in a signal input to an input port A and a pattern of change 
in a signal input to an input port B, such as those shown 
in FIG. 84, are stored in the pattern-of -signal -change 
storage 9002 beforehand. 

15 The element current storage means 9003 stores, for 

each circuit element and in a form such as that shown in 
FIG. 88, a changed pin concerning an electric current which 
flows in a circuit element in the event of a change arising 
in the circuit element when a change occurs in an external 

20 terminal of the circuit element, a change arising in the 
logic status of another external pin under the influence 
of the change in the external terminal, total electric 
current, and information concerning a peak value. 

25 For example, reference numeral 10601 designates a 
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circuit element: that is, a NAND circuit shown in FIG. 84, 
in which a total amount of 3 mA including a peak current 
of 1 mA flows into the circuit element, provided that an 
input signal A assumes a logic value of 1 and the logic 
5 value of an input signal B changes from 0 to 1 when a change 
arises in an output Y. 

The total current calculation means 9005 operates 
according to a flowchart, such as that shown in FIG. 89. 
According to the flowchart relating to a method of 

10 calculating an electric current shown in FIG. 89, 

processing pertaining to step 1504 of the conventional 
current calculation method is replaced by processing of 
calculating a rectangular waveform of the electric current 
from total electric current and a peak value at the time 

15 of occurrence of a change in consideration of the status 
of another terminal (step 10704) . 

In step 10701, there is read the circuit connection 
information which is shown in FIG. 84 and stored in the 
circuit connection information storage means 9001. In 

20 step 10702, there is read a pattern of signal change which 
is shown in FIG. 85 and is stored in the pattern-of- 
signal-change storage 9002. In step 10703, there is read 
the element current information which is shown in FIG. 88 
and stored in the element current storage means 9003. 

25 
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In step 10704, the pattern of signal shown in FIG. 
85 is imparted to the circuit connection information shown 
in FIG. 84, thereby enabling propagation of a signal. In 
a case where a change arises in a signal appearing at an 
5 external input terminal of a circuit element, which 

element is described in the information concerning the 
amount of electric current flowing in an element, a 
rectangular waveform of an electric current is produced 
by means of averaging the amount of electric current 

10 flowing through the external input terminal by means of 
a peak value in accordance with the status of another 
terminal. The thus-produced rectangular waveform is 
added to the information concerning the total mount of 
electric current obtained at a time when the change has 

15 arisen in the signal appearing at the external input 
terminal . 

After completion of propagation of a signal, in step 
10705 the thus-calculated information concerning the 
total amount of electric current obtained at each time 

20 shown in FIG. 90 is stored in the total current storage 
means 9004. As a result, there is produced a current 
waveform, such as that shown in FIG. 90, which is more 
realistic than that which is obtained by the conventional 
current calculation method and shown in FIG. 87. 

25 (Twentieth Embodiment) 



149 





FIG. 91 is a block diagram showing the overall 
flowchart of an EMI analysis method according to the 
present embodiment. The EMI analysis method according to 
the present embodiment is identical with the EMI analysis 
5 method according to the seventh embodiment described in 
connection with FIG. 28, except for addition of a cell 
characterization section G2 , replacement of the I-V table 
file AlO by an I-V table file (or an I-V function) G3 , and 
replacement of the voltage drop 
10 calculation/instantaneously-changing current correction 
section All by a line voltage drop 

calculation/instantaneously-changing current correction 
section G4 . 



15 section for producing information unique to each cell; G3 
denotes an I-V table file (or an I-V function) which 
represents, for each cell, the relationship between the 
amount of electric current and the degree of drop in line 
voltage; and G4 denotes a line voltage drop 

20 calculation/instantaneously-changing current correction 
section for calculating the degree of drop in line voltage 
and correcting the instantaneously-changing current 
calculated in step 11. 



Reference symbol G2 denotes a cell characterization 



[0226] 
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The operation of the cell characterization section 
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G2 will now be described. The cell characterization 
section G2 receives a library 3 and calculates in advance 
the amount of electric current at the time of a change 
arising in a line voltage of each cell, thereby producing 
5 an I-V table file (or an I-V function) for each cell. 

FIG. 92 is a block diagram showing the flow of 
processing of the voltage drop 

calculation/instantaneously-changing current correction 
section G4 according to the twentieth embodiment of the 

10 present invention. The function of the voltage drop 

calculation/instantaneously-changing current correction 
section G4 will now be described by reference to FIG, 92. 

Those portions of the EMI analysis method according 
to the present embodiment which are identical with those 

15 of the EMI analysis method according to the seventh 

embodiment are assigned the same reference numerals, and 
their repeated explanations are omitted. Further, since 
the data 12 have already been described in detail in 
connection with the first embodiment, the same reference 

20 numeral is assigned to the data 12, and repeated 

explanation is omitted. In contrast with the voltage drop 
calculation/instantaneously-changing current correction 
section All employed in the seventh embodiment, the 
voltage drop calculation/instantaneously-changing 

25 current correction section G4 is configured such that the 
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current value correction section B12 is replaced by step 
HI . 

There will now be described step HI. In step HI, 
there are entered the voltage drop information Bll 
5 prepared for each instance in step BIO according to the 
seventh embodiment and the I-V table file (or I-V function) 
which is produced for each cell by the cell 
characterization section G2 and represents the 
relationship between the amount of electric current and 

10 the degree of drop in line voltage. On the basis of a table 
file unique to each cell, the instantaneously-changing 
current information Bl , which is produced on the basis of 
the instantaneously-changing current result 12 produced 
in step 11, is corrected. An instantaneously-changing 

15 current calculation result A6 is produced for each 

instance in consideration of a drop in line voltage which 
is dependent on the resistance of a power line. 

As mentioned above, under the EMI analysis method 
20 according to the twentieth embodiment of the present 

invention, the relationship between the amount of electric 
current and the degree of drop in line voltage is 
characterized for each cell beforehand. In the seventh 
embodiment, EMI of a circuit element is analyzed on the 
25 basis of a uniform relationship between the amount of 
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electric current and the degree of drop in line voltage 
using the I-V table file {or I-V function) . In contrast 
with the EMI analysis method of the seventh embodiment, 
the EMI analysis method of the present embodiment enables 
5 dependence of EMI on the structure (or layout) of a cell 
to be reflected in EMI analysis. Therefore, the EMI 
analysis method of the present embodiment can produce a 
more accurate model of a current waveform, thereby 
reducing undesired higher harmonic components included in 

10 the frequency characteristic of a circuit to be analyzed, 
the characteristic being determined by means of subjecting 
the current waveform to FFT processing. 

Although the twentieth embodiment has described by 
reference to the first embodiment, the same result will 

15 be yielded even when the technical idea of the present 
embodiment is applied to any one of the fourteenth through 
nineteenth embodiments. 
(Twenty-first Embodiment) 

FIG. 93 is a block diagram showing the overall 

20 flowchart of an EMI analysis method according to a 

twenty-first embodiment of the present invention. The 
EMI analysis method according to the present embodiment 
is identical with the EMI analysis method according to the 
eighth embodiment described in connection with FIG. 36, 

25 except for addition of a cell characterization section II, 
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deletion of the C-V table file (or a C-V function) CI, 
replacement of the I-V table file AlO by an I -C-V table 
file (or an I-V function) 12, and replacement of the 
voltage drop calculation/instantaneously-changing 
5 current correction section C2 by a line voltage drop 
calculation/instantaneously-changing current correction 
section 13 . 

Those portions of the EMI analysis method according 
to the present embodiment which are identical with those 

10 of the EMI analysis method according to the first 

embodiment are assigned the same reference numerals, and 
their repeated explanations are omitted. Reference 
symbol II denotes a cell characterization section for 
producing information unique to each cell; 12 denotes an 

15 I-C-V table file (or an I-C-V function) which represents, 
for each cell, the relationship among the amount of 
electric current, the degree of drop in line voltage, and 
decoupling capacitance existing between a power line and 
ground; and 13 denotes a line voltage drop 

20 calculation/instantaneously-changing current correction 
section for calculating the degree of drop in line voltage 
and correcting the instantaneously-changing current 
calculation result 12 produced in step 11. 

The operation of the cell characterization section 

25 II will now be described. The cell characterization 
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section II receives a library 3, changes a line voltage 
of each cell, and calculates in advance the amount of 
electric current at the time of changing of the line 
voltage in consideration of decoupling capacitance 
5 existing between a power and ground, thereby producing an 
I-C-V table file (or an I-C-V function) 12 for each cell. 

FIG. 94 is a block diagram showing the flow of 
processing of the voltage drop 

calculation/instantaneously-changing current correction 

10 section 14 according to the twenty-first embodiment of the 
present invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 
section 13 will now be described by reference to FIG. 94. 
Those portions of the EMI analysis method according to the 

15 present embodiment which are identical with those of the 
EMI analysis method according to the eighth embodiment are 
assigned the same reference numerals, and their repeated 
explanations are omitted. Further, since the data 12 have 
already been described in detail in connection with the 

20 first embodiment, the same reference numeral is assigned 
to the data 12, and repeated explanation is omitted. 

In contrast with the voltage drop 
calculation/instantaneously-changing current correction 
section C2 employed in the eighth embodiment, the voltage 

25 drop calculation/instantaneously-changing current 
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correction section 13 is configured such that the current 
value correction section B12 is replaced by step Jl . 

There will now be described step Jl . In step 
Jl , there are entered the voltage drop information Bll 
5 prepared for each instance in step BIO according to the 
eighth embodiment, the I-C-V table file (or I-C-V 
function) 12 which is produced for each cell by the cell 
characterization section II and represents the 
relationship between the amount of electric current and 

10 the degree of drop in line voltage, and the capacitance 
information D2 concerning a capacitor provided in a 
segment interconnecting a power and ground, the 
information being prepared in step Dl in the eight 
embodiment. On the basis of a table file unique to each 

15 cell, the instantaneously-changing current information Bl , 
which is produced on the basis of the instantaneously- 
changing current result 12 produced in step 11, is 
corrected. An instantaneously-changing current 
calculation result A6 is produced for each instance, in 

20 consideration of fluctuations in line voltage which are 
dependent on decoupling capacitance as well as the 
resistance of a power line. 

As mentioned above, under the EMI analysis method 
according to the twenty-first embodiment of the present 

25 invention, the relationship among the amount of electric 
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current, the degree of drop in line vol tage , and decoupling 
capacitance existing between a power and ground is 
characterized for each cell beforehand. In the eighth 
embodiment, EMI of a circuit element is analyzed on the 
5 basis of a uniform relationship between the amount of 
electric current, the degree of drop in line voltage, and 
decoupling capacitance existing between a power and ground, 
using the I-V table file (or I-V function) AlO and the C-V 
table (or C-V function) CI. In contrast with the EMI 

10 analysis method of the eighth embodiment, the EMI analysis 
method of the present embodiment enables dependence of EMI 
on the structure (or layout) of a cell to be reflected in 
EMI analysis. Therefore, the EMI analysis method of the 
present embodiment can produce a more accurate model of 

15 a current waveform, thereby reducing undesired higher 
harmonic components included in the frequency 
characteristic of a circuit to be analyzed, the 
characteristic being determined by means of subjecting the 
current waveform to FFT processing. 

20 Although the twenty-first embodiment has described 

by reference to the first embodiment, the same result will 
be yielded even when the technical idea of the present 
embodiment is applied to any one of the fourteenth through 
nineteenth embodiments. 

25 (Twenty-second Embodiment) 
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FIG. 95 is a block diagram showing the overall 
flowchart of an EMI analysis method according to a 
twenty-second embodiment of the present invention. The 
EMI analysis method according to the present embodiment 
5 is identical with the EMI analysis method according to the 
seventh embodiment described in connection with FIG. 28, 
except for addition of a comparison-and-de termination 
section Kl . Those portions of the EMI analysis method 
according to the present embodiment which are identical 
10 with those of the EMI analysis method according to the 
first embodiment are assigned the same reference numerals , 
and their repeated explanations are omitted. Reference 
symbol Kl denotes a comparison-and-determination section 
for determining whether or not the instantaneously- 
15 changing current calculation result A6 — ^which takes into 
consideration a drop in line voltage produced by the 
voltage drop calculation/instantaneously-changing 
current correction section All — becomes equal to an input 
to the next processing step. 
20 Next will be described the operation of the 

comparison-and-determination section Kl . The EMI 
analysis method according to the present embodiment is 
characterized by means of iterating the following 
operations; that is, re-entry, into the voltage drop 
25 calculation/instantaneously-changing current correction 
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section All, of the instantaneously-changing current 
calculation result A6 which takes into consideration a 
drop in line voltage produced by the voltage drop 
calculation/instantaneously-changing current correction 
5 section All, thereby producing another 

instantaneously-changing current calculation result A6 
which newly takes into consideration a drop in line 
voltage . 

The comparison-and-de termination section Kl 

10 determines, through comparison, whether or not the number 
of times the foregoing operations has been iterated has 
reached a predetermined number of times. If a truth-value 
TRUE is yielded as a result of comparison, iteration of 
the operations is completed, and there is performed 

15 processing pertaining to the next step; that is, the 

operation of the FFT processing section 9, is performed. 

As mentioned above, under the EMI analysis 
method according to the twenty-second embodiment of the 
present invention, the instantaneously-changing electric 

20 current result which has been corrected and calculated in 
consideration of the influence of fluctuations in line 
voltage is fed back, to thereby iterate voltage drop 
calculation/instantaneously-changing current correction 
operations. Accordingly, variation in the amount of 

25 electric current , which would otherwise be caused by a drop 



159 



in line voltage, can be made stable. Therefore, the EMI 
analysis method of the present embodiment can produce a 
more accurate model of a current waveform, thereby 
reducing undesired higher harmonic components included in 
5 the frequency characteristic of a circuit to be analyzed, 
the characteristic being determined by means of subjecting 
the current waveform to FFT processing. 

In the present embodiment, the comparison-and- 
determination section Kl controls iteration of processing 

10 operations by means of determining whether or not a 
predetermined number of times has been reached. 
Alternatively, there may be employed a storage device for 
temporarily storing the instantaneously-changing current 
calculation result A6 which takes into consideration a 

15 drop in line voltage produced by the voltage drop 

calculation/instantaneously-changing current correction 
section All . The instantaneously-changing current 
calculation result A6 stored in the storage device is 
entered again into the voltage drop 

20 calculation/instantaneously-changing current correction 
section All, to thereby newly produce another 
instantaneously-changing current calculation result A6 . 
Comparison is made between the thus-produced another 
instantaneously-changing current calculation result A6 

25 and the instantaneously-changing current calculation 
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result A6 stored in the storage device. The foregoing 
operations are iterated until a difference falls within 
a predetermined range of allowance. 

Although the twenty-second embodiment has been 
5 described with reference to the seventh embodiment, the 
same result will be yielded even when the technical idea 
of the present embodiment is applied to any one of the 
eighth, ninth, twentieth, and twenty-first embodiments. 

10 (Twenty- third Embodiment) 

An EMI analysis method according to a tenth 
embodiment of the present invention will now be described. 
In the fourteenth through twenty-second embodiments, an 
instantaneously-changing current of each cell is 

15 calculated from event information concerning each cell, 
the information being output from a logic simulator, by 
means of characterizing a peak current of a cell according 
to the type of cell and the amount of instantaneously- 
changing current for each event. The thus-calculated 

20 instantaneously-changing current is modeled as a 

rectangular waveform such that a peak value becomes equal 
to a peak current of the cell, thereby determining change 
in electric current. However, so long as the waveform of 
an instantaneous electric current is modeled as a 

25 triangular wave, a more realistic change in electric 
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current can be determined. 

By way of example, there will now be described 
an EMI analysis method comprising the steps of: 

modeling an instantaneously-changing electric 
5 current as a triangular wave whose height becomes equal 
to a characterized peak current, the width of the 
triangular wave being calculated such that the 
instantaneously-changing current becomes equal to the 
width of the triangular wave, to thereby calculate a more 
10 realistic change in electric current; and 

subjecting the result of calculation of an electric 
waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
of interest. 

15 

FIG. 96 shows an instantaneous electric current and 
a modeled triangular wave according to the present 
embodiment. In the drawing, ^^i" denotes an 
instantaneously-changing current; and "^^t" denotes a 

20 current processing time. Further, ^^h" denotes the height 
of a triangular wave serving as a model of an 
instantaneously-changing current. The height of a 
triangular waveform is extracted from a library in which 
the heights of triangular waveforms of cells are 

25 characterized. Reference symbol "w" denotes the width of 
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a triangular waveform and is calculated on the basis of 
^^i" and "h." An electric current at time t+x is expressed 
by Eq.5 described previously. Further, an electric 
current at time t-x is also expressed by Eq.5. 
5 The overall flow of processing of the EMI analysis 

method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1). 

FIG. 97 is a detailed block diagram showing the 

10 current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 

15 information being output from the logic simulator 4 . 

Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an ins tantaneously- 

20 changing current (hereinafter referred to as an 

^^instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 

25 As can be seen from FIG. 3, according to the 
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instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 6 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
5 identical with the current wave calculation operation of 
the first embodiment. Reference numeral 130 designates 
a peak current library including peak currents of cells. 
Reference numeral 160 denotes a triangular waveform height 
calculation section for calculating the height of a 

10 triangular wave from the peak current library 130. 

Reference numeral 161 designates a triangular waveform 
width calculation section for calculating the width of a 
triangular waveform from the height extracted by the 
triangular waveform height calculation section 160 and the 

15 instantaneously-changing current, such that the area of 
a triangular wave becomes equal to an instantaneously- 
changing current. Reference numeral 162 denotes a 
triangular-waveform-shaping section for modeling the 
instantaneously-changing current result 12 as a 

20 triangular wave having a constant width, the height of the 
triangular waveform corresponding to the peak value of the 
triangular waveform calculated by the rectangular 
waveform height calculation section 160, and the width of 
the rectangular waveform corresponding to the width 

25 calculated by the rectangular waveform width calculation 
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section 161. 

FIG. 98 is a flowchart of processing pertaining to 
the triangular— waveform-shaping section 162, including 
processing pertaining to the triangular waveform height 

5 calculation section 160 and processing pertaining to the 
triangular waveform width calculation section 161 . First, 
the peak current library 130 is read (step 1160) . The 
following processing is iterated on a per-event basis 
until calculation of a current waveform is completed (step 

10 1161) . From the peak current library 130, a peak current 
of an instance which is an object of EMI analysis is 
extracted, and the thus-extracted peak current is taken 
as the height "h" of a triangular waveform (step 1162) . 
This processing corresponds to processing of the 

15 triangular waveform height calculation section 160. 

Provided that ^^w" is taken as the width of a triangular 
waveform, the area of a triangular wave is Wxh/2 and is 
equal to the instantaneously-changing current of a cell 
of the event. Therefore, the width "w" of a triangular 

20 waveform is calculated from 2xi/h (step 1163) . Reference 
symbol "i" denotes an instantaneously-changing current of 
a cell of an event, the event being an object of EMI analysis . 
This processing corresponds to processing pertaining to 
the triangular waveform width calculation section 161. 

25 Until variable "x" changes from 0 to W/2 , h(c, i) expressed 
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by Eq . 5 is added to I ( t+x) and I ( t-x) . Further , At i s added 
to variable ^^x'' (steps 1164 and 1165). Here, I (t+x) 
denotes total electric current flowing through all the 
cells at time t+x, and I (t-x) denotes total electric 
5 current flowing through all the cells at time t-x. This 
processing corresponds to processing pertaining to the 
triangular-waveform-shaping section 162. 

The current wave calculation result 8 can be obtained 
10 by means of the flow of a triangular waveform height 
calculation operation, a triangular waveform width 
calculation operation, and a triangular-waveform-shaping 
operation. The FFT processing section 9 subjects the 
current waveform calculation result 8 to an FFT operation, 
15 to thereby determine the frequency characteristic of an 
EMI component. As a result, there can be obtained an EMI 
analysis result 10. 

(Twenty-fourth Embodiment) 

20 An EMI analysis method according to a twenty-fourth 

embodiment of the present invention will now be described. 
In the fourteenth through twenty-second embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 

25 the information being output from a logic simulator, by 
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means of characterizing a peak current of a cell according 
to the type of cell and the amount of instantaneously- 
changing current for each event. The thus-calculated 
instantaneously-changing current is modeled as a 
5 rectangular waveform such that a peak value becomes equal 
to a peak current of the cell, thereby determining change 
in electric current. However, so long as the waveform of 
an instantaneous electric current is modeled as a Gauss 
function wave, a more realistic change in electric current 

10 can be determined. By way of example, there will now be 
described an EMI analysis method comprising the steps of: 

modeling an instantaneously-changing electric 
current as a multi-order-function waveform whose height 
becomes equal to a characterized peak current, the width 

15 of the multi-order-function wave being calculated such 
that the instantaneously-changing current becomes equal 
thereto, to thereby calculate a more realistic change in 
electric current; and 

subjecting the result of calculation of an electric 

20 waveform to FFT processing, thereby determining the 

frequency characteristic of an EMI component of a circuit 
of interest. 

FIG. 99 shows an instantaneous electric current and 
25 a modeled mul ti -order-function according to the present 
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embodiment. In the drawing, '"i" denotes an 
instantaneously-changing current; and '^t" denotes a 
current processing time. Further, ^^h" denotes the height 
of a multi-order-function wave serving as a model of an 

5 instantaneously-changing current. The height of a 

triangular waveform is extracted from a library in which 
the heights of multi-order-function waveforms of cells are 
characterized. Reference symbol "U" denotes the width of 
a multi-order-function wave and is calculated on the basis 

10 of "h," which has been calculated previously, and "i . " 
Provided that the width of the multi-order-function 
waveform is taken as 2xw, an electric current at time t+x 
is expressed by Eq.6. The area of the multi-order- 
function waveform is expressed by Eq.8. The area of the 

15 multi-order-function waveform is equal to the value of 
"i." Here, "a" is expressed by Eq.9. Since the height 
^''h" of a multi-order-function waveform is determined by 
Eq.6 when x=0 , the height ^^h'' is expressed by Eq.7. From 
Eq. 7 and Eq. 9, the width ^^w" is expressed by Eq. 14 provided 

20 below. 

w=^... (14) 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
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by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 100 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 

5 present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 

10 Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 

15 "instantaneously-changing current calculation result 

12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 

instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 

20 instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 6 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 

25 the first embodiment. 
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Reference numeral 130 designates a peak current 
library including peak currents of cells which are 
characterized according to the type of cell. Reference 
5 numeral 170 denotes a multi -order-function waveform 

height calculation section for calculating the height of 
a multi-order-function waveform from the peak current 
library 130. Reference numeral 171 designates a 
multi-order-function waveform width calculation section 
10 for calculating the width of a multi-order-function 

waveform from the height extracted by the multi -order- 
function waveform height calculation section 17 0 and the 
instantaneously-changing current, such that the area of 
a multi-order-function waveform becomes equal to an 
15 instantaneously-changing current . Reference numeral 172 
denotes a multi-order-function waveform-shaping section 
for modeling the instantaneously-changing current result 
12 into a multi-order-function waveform having a constant 
width, the height of the multi-order-function waveform 
20 corresponding to the peak value of the multi-order- 
function waveform calculated by the multi-order-function 
waveform height calculation section 170, and the width of 
the multi -order-function waveform corresponding to the 
width calculated by the mul ti -order-function waveform 
25 width calculation section 171. 
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FIG. 101 is a flowchart of processing pertaining to 
the multi -order-function waveform-shaping section 172, 
including processing pertaining to the multi-order- 
5 function waveform height calculation section 170 and 
processing pertaining to the multi-order-function 
waveform width calculation section 171. First, the peak 
current library 130 is read (step 1170) . The following 
processing is iterated on a per-event basis until 

10 calculation of a current waveform is completed (step 1171) . 
From the peak current library 130, a peak current of an 
instance which is an object of EMI analysis is extracted, 
and the thus -extracted peak current is taken as the height 
■^^h" of a multi-order-function waveform (step 1172) . This 

15 processing corresponds to processing of the multi- 
order-function waveform height calculation section 170. 
The width ^^w" of a multi -order-function waveform is 
calculated by Eq.l4 (step 1173). Reference symbol ^'i" 
denotes an instantaneously-changing current of a cell of 

20 an event, the event being an object of EMI analysis. This 
processing corresponds to processing pertaining to the 
multi-order-function waveform width calculation section 
171. Until variable ^^x" changes from -w to w, '^^y" 
expressed by Eq. 6 is added to I (t+x) and I (t-x) . Further, 

25 At is added to variable ^^x" (steps 1174 and 1175) . Here, 
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I (t+x) denotes total electric current flowing through all 
the cells at time t+x, and I (t-x) denotes total electric 
current flowing through all the cells at time t-x. This 
processing corresponds to processing pertaining to the 
5 multi -order-function waveform-shaping section 172. 
[0255] 

The current wave calculation result 8 can be obtained 
by means of the flow of a multi-order-function waveform 
height calculation operation, a multi -order-function 

10 waveform width calculation operation, and a multi- 
order-function waveform- shaping operation. The FFT 
processing section 9 subjects the current waveform 
calculation result 8 to an FFT operation, to thereby 
determine the frequency characteristic of an EMI component. 

15 As a result, there can be obtained an EMI analysis result 
10. 

(Twenty-fifth Embodiment) 

An EMI analysis method according to a twenty-fifth 

20 embodiment of the present invention will now be described. 
In the fourteenth through twenty-second embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 
the information being output from a logic simulator, by 

25 means of characterizing a peak current of a cell according 



172 




to the type of cell and the amount of instantaneously- 
changing current for each event. The thus -calculated 
instantaneously-changing current is modeled as a 
rectangular waveform such that a peak value becomes equal 
5 to a peak current of the cell, thereby determining change 
in electric current. However, so long as the waveform of 
an instantaneous electric current is modeled as a Gauss 
function waveform, a more realistic change in electric 
current can be determined. 

10 By way of example, there will now be described an 

EMI analysis method comprising the steps of : 

modeling, for each type of cell, an 
instantaneously-changing electric current as a Gauss 
function waveform whose height and width are characterized, 

15 to thereby calculate a more realistic change in electric 
current; and 

subjecting the result of calculation of an electric 
waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
20 of interest. 

FIG. 50 shows an instantaneous electric current and 
a modeled Gauss function according to the present 
embodiment. In the drawing, '^^i'' denotes an 
25 instantaneously-changing current; ^^t" denotes a current 
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processing time; and ^^h'' denotes the height of a Gauss 
function wave serving as a model. The height ^^h" is equal 
to a peak current characterized according to the type of 
cell. Reference symbol "w" denotes a width of a Gauss 

5 function serving a model, and the width ^^w" corresponds 
to a value characterized according to the type of cell. 
In the present embodiment, a Gauss function serving as a 
model is expressed as Eq.ll. This function is called a 
Gaussian function. Let x=0 , Eq.ll provides ^^h,'' which 

10 represents the height of a Gaussian function waveform. 
From Eq.ll, letx=3xW, there is derived a value of O.OOOlxh 
(i.e., 0.01% of the height of a Gaussian function wave) . 
This value is taken as a limit on the width of a Gaussian 
function serving as a model, and the width of the Gaussian 

15 function is taken as ^^6xw." 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
20 embodiment (see FIG. 1) . 

FIG. 102 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
25 calculating a momentary change in the current of each cell 
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from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
5 conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
^^instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
10 calculation section 11. FIG. 3 shows an example 

instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 6 function. 
15 Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. A Gaussian function waveform 
height calculation section 180 extracts a peak value of 
20 a Gaussian function serving as a model, from the peak 
current library 130 including peak currents of cells 
characterized according to the type of cell. Reference 
numeral 181 designates a wave width library includes wave 
widths of cells characterized according to the type of cell . 
25 Reference numeral 182 designates a Gaussian function 
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waveform width calculation section for extracting, from 
the wave width library 181, the width of a Gaussian 
function waveform serving as a model. 

5 Reference numeral 183 designates a Gaussian function 

waveform- shaping section for modeling the 
instantaneously-changing current result 12 as a Gaussian 
function wave having a peak value calculated by the 
Gaussian function waveform height calculation section 180 
10 and a width calculated by the Gaussian function waveform 
width calculation section 182. 

FIG. 103 shows the flowchart of processing 
pertaining to the Gaussian function waveform-shaping 

15 section 183, including processing pertaining to the 

Gaussian function waveform height calculation section 180 
and processing pertaining to the Gaussian function 
waveform width calculation section 182. First, the peak 
current library 130 is read (step 1180) , and the wave width 

20 library 181 is read (step 1181) . 

The processing provided below is iterated on a 
per-event basis until calculation of a current waveform 
is completed (step 1182) . From the peak current library 
25 130, a peak current of an instance which is an object of 
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EMI analysis is extracted, and the thus-extracted peak 
current is taken as the height "h" of a Gaussian function 
waveform (step 1183) . This processing corresponds to 
processing of the Gaussian function waveform height 
5 calculation section 180. The width '^w" of a Gaussian 
function waveform which is an object of EMI analysis is 
extracted from the wave width library 181, and the 
thus-extracted width is taken as the width "w" of a 
Gaussian function waveform (step 1184) . This processing 

10 corresponds to processing of the Gaussian function 

waveform width calculation section 182. Until variable 
^^x" changes from -3xw to 3xw, '''y" expressed by Eq.ll is 
added to I(t+x). Further, At is added to variable ^^x" 
(steps 1185 and 1186). Here, I(t+x) denotes total 

15 electric current flowing through all the cells at time t+x . 
This processing corresponds to processing pertaining to 
the Gaussian function waveform-shaping section 183. 

The current wave calculation result 8 can be obtained 
20 by means of the flow of a Gaussian function waveform height 
calculation operation, a Gaussian function waveform width 
calculation operation, and a Gaussian function 
waveform- shaping operation. The FFT processing section 
9 subjects the current waveform calculation result 8 to 
25 an FFT operation, to thereby determine the frequency 
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characteristic of an EMI component. As a result, there 
can be obtained an EMI analysis result 10. 

(Twenty-sixth Embodiment) 

5 An EMI analysis method according to a twenty-sixth 

embodiment of the present invention will now be described - 
The EMI analysis method according to the present 
embodiment is characterized by modeling an 
instantaneously-changing current while the information 

10 concerning a current waveform of each circuit element (or 
cell) which is described in connection with the fourteenth 
through twenty-second embodiments is separated into an 
element dependent on a short circuit current (also called 
"short circuit-current-dependent element") and another 

15 element dependent on a charge current (also called 

"charge-current-dependent element") . In contrast with 
the thirteenth embodiment, the present embodiment effects 
modeling of each of the circuit elements (cells) , thereby 
enabling formation of a more accurate model of electric 

20 current. 

For example, in a case where the circuit connection 
information shown in FIG. 110 and the pattern of signal 
change shown in FIG. 12 are imparted to a circuit element, 
25 an electric current shown in FIG. 104 flows through the 
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circuit element. In such a case, the conventional EMI 
analysis method enables formation of only a waveform such 
as that shown in FIG. 14. Even in the case of the 
thirteenth embodiment, there is produced only a waveform 
5 such as that shown in FIG. 58. In contrast, the EMI 
analysis method according to the present embodiment 
enables formation of a more correct electric current . FIG . 
72 shows the configuration of an EMI analysis method 
according to the present embodiment, the method using the 
10 value of a peak current and the amount of electric current 
described in connection with the fourteenth embodiment. 

A system to be used with the EMI analysis method shown 
in FIG. 72 comprises circuit connection information 

15 storage means 9001; pattern-of -signal-change storage 

9002/ element current storage means 9003 for storing the 
amount of electric current flowing in a circuit element; 
total current storage means 9004 for storing a total amount 
of electric current flowing through the circuit element; 

20 and total current calculation means 9005 for calculating 
the total amount of electric current flowing through a 
circuit element. 

Next will be described individual elements 
25 constituting the EMI analysis method shown in FIG. 72 and 
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procedures for calculating the total amount of electric 
current shown in FIG. 106, from the circuit connection 
information shown in FIG. 110 , the pattern of signal change 
shown in FIG. 12, and the amount of electric current 
5 flowing through a circuit element shown in FIG. 105. 

The circuit connection information storage means 
9001 stores information concerning the connection of a 
circuit to be analyzed in the same manner as described in 
10 connection with the background art. Circuit connection 
information, such as that shown in FIG. 110, is stored in 
advance in the circuit connection information storage 
means 9001 . 

15 The pattern-of -signal -change storage 9002 stores 

the pattern of signal change to be applied to the circuit 
connection information, as in the case of a conventional 
pattern of signal change described in connection with the 
background art. As shown in FIG. Ill, a pattern of change 

20 in a signal input to an input port CK and a pattern of change 
in a signal input to an input port D, such as those shown 
in FIG. 110, are stored in the pattern-of -signal-change 
storage 9002 beforehand. 

25 The element current storage means 9003 stores^, for 
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each circuit element and in a form such as that shown in 
FIG. 105, the total amount of electric current flowing in 
a circuit element (cell) , a peak value of the electric 
current, the total amount of charge current, and a peak 
5 current. 

The total current calculation means 9005 operates 
according to a flowchart, such as that shown in FIG. 107. 
According to the flowchart relating to a method of 
10 calculating an electric current shown in FIG. 107, 

processing pertaining to step 1504 of the conventional 
current calculation method is replaced by processing of 
calculating a rectangular waveform of the short circuit 
current from a peak value and the total amount of electric 
15 current, calculating a rectangular waveform of the charge 
current from a peak value and the total amount of electric 
current, and merging the thus -calculated rectangular 
waveforms into a single rectangular waveform (step 12504) . 
In step 12501, there is read the circuit connection 
20 information which is stored in the circuit connection 
information storage means 9001 and shown in FIG. 110. 

In step 12502, there is read a pattern of signal 
change which is shown in FIG. Ill and is stored in the 
25 pattern-of-signal-change storage 9002. In step 12503, 
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there is read the element current information which is 
shown in FIG. 105 and stored in the element current storage 
means 9003. In step 12504, the pattern of signal shown 
in FIG. Ill is imparted to the circuit connection 
5 information shown in FIG. 110, thereby enabling 

propagation of a signal. In a case where a change arises 
in a signal appearing at an external input terminal of a 
circuit element, which element is described in the 
information concerning the amount of electric current 

10 flowing in an element, a rectangular waveform of an 

electric current is produced by means of averaging a short 
circuit current component of the electric current by a peak 
value of the short circuit current component, and a 
rectangular waveform of an electric current is produced 

15 by means of averaging a charge current component of the 
electric current by a peak value of the charge current 
component. The thus-formed rectangular waveforms are 
added to the information concerning the total mount of 
electric current obtained at a time when the change has 

20 arisen in the signal appearing at the external input 

terminal. After completion of propagation of a signal, 
in step 12505 the thus-calculated information concerning 
the total amount of electric current obtained at each time 
shown in FIG. 106 is stored in the total current storage 

25 means 9004. As a result, there is produced a current 
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waveform, such as that shown in FIG. 106, which is more 
realistic than that which is obtained by the conventional 
current calculation method and shown in FIG, 14 . 

5 (Twenty- seventh Embodiment) 

Next will be described an EMI analysis method 
according to a twenty- seventh embodiment of the present 
invention. The tenth embodiment has described that a more 
realistic change in current can be materialized by means 

10 of calculating an instantaneously-changing current of 
each cell from event information concerning the cell, the 
information being output from a logic simulator, and 
modeling the thus-calculated instantaneously-changing 
current as a triangular wave. The EMI analysis method 

15 according to the present embodiment is characterized by 
calculating a current waveform in consideration of output 
slew information; that is, the time which elapses from the 
time a change arises in a signal appearing at an output 
terminal of a cell until the change has been canceled. 

20 By way of example, there will now be described an 

EMI analysis method comprising the steps of: 

modeling an instantaneously-changing current as a 
triangular waveform by means of calculating the width of 
a triangular waveform for each event information in 

25 consideration of slew information concerning an output 
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terminal of a cell and the height of the triangular 
waveform on the basis of the thus-calculated width such 
that the area of the triangular waveform becomes equal to 
the amount of current of each event; 
5 calculating a more realistic change in electric 

current by means of using an existing library in its 
present form; and 

subjecting the result of calculation of an electric 
waveform to FFT processing, thereby determining the 
10 frequency characteristic of an EMI component of a circuit 
of interest. 

FIG. 117 shows a logic value (see FIG. 117A) and a 
modeled triangular waveform {see FIG. 117B) , which are to 

15 be used for executing the EMI analysis method according 
to the present embodiment. Reference symbol ^^t" denotes 
a current processing time, and "W" denotes the width of 
a triangular waveform serving as a model of an 
instantaneously-changing current. The width ^''W" is 

20 determined by means of calculating an instantaneously- 
changing current of each cell from event information of 
each cell, in consideration of output slew information. 
Further, ^^h" denotes the height of a triangular waveform 
and is calculated from "i" and ''^W." Here, h=2xl/W. 

25 FIG. 118 shows the configuration of a system for 
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effecting the EMI analysis method according to the present 
invention . The system used for effecting the EMI analysis 
method shown in FIG. 118 comprises pattern-of -signal - 
change storage means 1171; element current storage means 

5 1172 for storing the amount of electric current flowing 
in a circuit element; output slew information storage 
means 1173; and current waveform calculation means 1175 
for calculating a current waveform on the basis of data 
output from the pattern-of -signal-change storage means 
10 1171, the element current storage means 1172, the output 
slew information storage means 1173, and circuit 
connection information storage means 1174. The system 
outputs a calculation result of a current waveform 117 6 
(hereinafter referred to simply as a "current waveform 

15 calculation result 1176") . 

FIG. 119 is a flowchart of a triangular- 
waveform- shaping operation, including a triangular 
waveform height calculation operation and a triangular 

20 waveform width calculation operation. The following 

processing is iterated on a per-event basis after reading 
of an output slew (step 1190) until calculation of a 
current waveform is completed (step 1191) . A base of the 
triangular waveform is extracted as an output slew while 

25 a point at which a change arises in an event is taken as 
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a center (step 1192). An isosceles triangle whose area 
is equal to the amount of electric current flowing through 
the event is taken as a current waveform of the event. The 
area of a triangular wave is Wxh/2 and assumes a value "i . " 
5 The height ^^h" of a triangular wave is calculated from 
2xi/W (step 1193) . Reference symbol ^^i" denotes an 
instantaneously-changing current of a cell of an event 
which is an object of EMI analysis. This processing 
corresponds to processing pertaining to the triangular 

10 waveform height calculation section 100 according to the 
tenth embodiment. Until variable "x" changes from 0 to 
W/2, (2xh/W)«x+h is added to I (t+x) and I (t-x) . Further, 
At is added to variable ^^x" (steps 1194, 1195, and 1196) . 
Here, I (t+x) denotes total electric current flowing 

15 through all the cells at time t+x, and I (t-x) denotes total 
electric current flowing through all the cells at time t-x. 
This processing corresponds to processing pertaining to 
the triangular-waveform-shaping section 101. 

20 The current wave calculation result 7 6 can be 

obtained by means of the flow of a triangular waveform 
height calculation operation and a triangular- 
waveform-shaping operation. The FFT processing section 
9 subjects the current wave calculation result 76 to an 

25 FFT operation, to thereby determine the frequency 
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characteristic of an EMI component. As a result, there 
can be obtained an EMI analysis result 10. 

The EMI analysis method enables estimation of an 
5 electric current with sufficient accuracy from only the 
information which is currently used by a power consumption 
analysis tool and used for verifying a delay. 

(Twenty-eighth Embodiment) 

10 Next will be described an EMI analysis method 

according to a twenty-eighth embodiment of the present 
invention. The EMI analysis method according to the 
present embodiment is characterized by means of 
calculating a current waveform, in consideration of 

15 whether an event is a rise or a fall and in addition to 
the feature of the EMI analysis method according to the 
twenty- seventh embodiment. FIG. 120 shows a logic value 
(see FIG. 120A) and a modeled triangular wavef orm (see FIG. 
120B) , which are to be used for effecting an EMI analysis 

20 method according to the present embodiment. Reference 
symbol '^t" denotes a current processing time, and ''^W" 
denotes the width of a triangular waveform serving as a 
model of an instantaneously-changing current. The width 
"W" is determined by means of calculating an 

25 instantaneously-changing current of each cell from event 
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information of each cell, in consideration of output slew 
information- There is employed rise-dependent 
coefficient/f all-dependent coefficient storage means 
1207, and current waveform calculation means 1205 
5 calculates a current waveform in accordance with whether 
an event is a rise or a fall. 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 

10 by the same block, diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 121 shows the configuration of a system for 
effecting the EMI analysis method according to the present 
invention . The system used for effecting the EMI analysis 

15 method shown in FIG. 121 comprises pattern-of -signal - 
change storage means 1201; element current storage means 
1202 for storing the amount of electric current flowing 
in a circuit element; output slew information storage 
means 1203; and current waveform calculation means 1205 

20 for calculating a current waveform on the basis of data 
output from the pattern-of -signal-change storage means 
1201, the element current storage means 1202, the output 
slew information storage means 1203, and circuit 
connection information storage means 1204 . The system 

25 outputs a calculation result of a current waveform 1206 
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(hereinafter referred to simply as a ^^current waveform 
calculation result 1206") . 

In the present embodiment, a current waveform is 
5 calculated in consideration of whether an event is a rise 
or a fall, by reference to the rise/fall-dependent 
coefficient storage means 1207, in addition to the 
operations of the EMI analysis method described in 
connection with the twenty- seventh embodiment. A more 
10 realistic change in current is calculated by use of an 
existing library in its present form, and by subjecting 
a calculation result of a current wave to FFT processing, 
to thereby determine the frequency characteristic of an 
EMI component and enable EMI analysis. 

15 

(Twenty-ninth Embodiment) 

Next will be described an EMI analysis method 
according to a twenty-ninth embodiment of the present 
invention. According to the present method, there is 

20 estimated the amount of electric current on the basis of 
event information, the number of cells determined 
according to a SPICE netlist provided in a cell, an output 
load capacitance, input slew information, the 
relationship between an input slew and the amount of 

25 current, and a time at which a change has arisen in an input 
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signal, or a delay time. Information concerning the 
amount of" electric current which flows in a single event 
is calculated from an output load capacitance and a table 
representing the relationship between an input slew and 
5 the amount of electric current. FIG. 122 shows a logic 
value (see FIG. 122A) and a modeled triangular waveform 
(see FIG. 122B) , which are to be used for effecting an EMI 
analysis method according to the present embodiment. 
Reference symbol "t^" designates a delay time; ^^t" denotes 

10 a current processing time; and ^^W" denotes the width of 
a triangular waveform serving as a model of an 
instantaneously-changing current. The width '■^Vf" is 
determined by means of calculating an instantaneously- 
changing current of each cell from event information of 

15 each cell, on the basis of output slew information taking 
into consideration a delay time td caused by multistage 
connection. Here, there is calculated the amount of 
current independent of an output load (i.e., the amount 
of non-dependent current) beforehand, and the thus- 

20 calculated amount of current is stored. A base of the 
triangular waveform is taken as an output slew while a 
point at which a change arises in an output is taken as 
a center. Further, there are employed, as output loads, 
an angular waveform whose area is equal to the amount of 

25 non-dependent current or to the amount of current flowing 
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when 

electric ciarrent-capacitance is 0 and a triangular 
waveforms equal in number to the number of cells within 
the delay time. As shown in FIG. 123, the EMI analysis 
5 method according to the present embodiment is provided 
with a load capacitance/slew-dependent element current 
information table 1226; element current calculation means 
1227 for calculating the amount of electric current 
flowing through a cell from the load capacitance/slew- 

10 dependent element current information table 1226; and 
non-load element current storage means 1228 for storing 
the amount of electric current flowing through a non-load 
element calculated by the element current calculation 
means 1227. By use of these means, a current waveform is 

15 calculated. The systemwhich is shown in FIG. 123 and used 
for effecting the EMI analysis method according to the 
present invention, comprises the load 

capacitance/slew-dependent element current information 
table 122 6; the element current calculation means 1227 for 

20 calculating the amount of electric current flowing through 
a cell from the load capacitance/slew-dependent element 
current information table 1226; the non-load element 
current storage means 1228 for storing the amount of 
electric current flowing through a non-load element 

25 calculated by the element current calculation means 1227; 
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pattern-of-signal-change storage means 1221; element 
current storage means 1222 for storing the amount of 
electric current flowing in a circuit element at the time 
of an event, the electric current being calculated by the 
5 element current calculation means 1227; output slew 
information storage means 1223; and current waveform 
calculation means 1225 for calculating a current waveform 
on the basis of data output from the pattern-of- 
signal-change storage means 1221, the element current 

10 storage means 1222, the output slew information storage 
means 1223, circuit connection information storage means 
1224, number-of- stages storage means 1229 for storing the 
number of stages in a composite cell, and the non-load 
element current storage means 1228 . The system outputs 

15 a calculation result of a current waveform 1206 

(hereinafter referred to simply as a ^^current waveform 
calculation result 1206") . 

In the twenty-ninth embodiment, there are stored a 
20 plurality of triangular waveforms which are equal in 

number to cells, such that the total area of a triangular 
waveform becomes equal to the amount of electric current 
obtained when a capacitance assumes a value of 0, during 
a period of time from the time when a change has arisen 
25 in an input signal until a change arises in an event. A 
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more realistic change in current is calculated from an 
existing library in its present form, and the calculation 
result of current waveform is subjected to FFT processing, 
thereby determining the frequency characteristic of an EMI 
5 component arising in a circuit to be analyzed and enabling 
EMI analysis. 



The present invention evaluates electromagnetic 
10 interference arising in an LSI through a simulation, by 
means of accurately analyzing a change in a power current, 
which change can be said to primarily account for 
electromagnetic interference. In order to solve the 
problem, a highly-accurate power current analysis 
15 function is added to a gate-level simulation, thereby 
preventing consumption of much time, which would otherwise 
be required by the conventional EMI analysis method using 
transistor-level current analysis. 
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What is claimed is: 

1. An electromagnetic interference analysis method 
for analyzing the amount of electromagnetic interference 
arising in an LSI by means of performing a logic simulation, 
5 the method comprising: 

an instantaneous current calculation step of 
calculating the amount of instantaneous electric current 
from event information, the information being produced 
when a change arises in a signal and including the instance 

10 name of each cell in which the change has arisen, the name 
of the signal, a time at which the change has arisen, and 
transition information; 

a modeling step of modeling the instantaneous 
electric current according to a predetermined rule; and 

15 an FFT processing step of subjecting to fast Fourier 

processing (hereinafter referred to as "FFT processing") 
the information concerning a change in electric current, 
the information being calculated through a modeling step. 

2 . The electromagnetic interference analysis method 

20 as defined in claim 1, wherein the modeling step includes 
an averaging step of averaging the instantaneous current 
over a certain discrete width, and the FFT processing step 
includes a step of subjecting to FFT processing 
information concerning a change in current, the 

25 information being produced by the averaging step. 
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3. The electromagnetic interference analysis method 
as defined in claim 1, wherein the modeling step includes 
a rectangular waveform modeling step of modeling the 
instantaneous current as a rectangular waveform whose 
5 height is calculated for information for each event such 
that the area of the rectangular waveform becomes equal 
to the amount of electric current of each event, and the 
FFT processing step includes a step of subjecting to FFT 
processing information concerning a change in current, the 
10 information being calculated in the rectangular waveform 
modeling step. 

4 . The electromagnetic interference analysis method 
as defined in claim 1, wherein the modeling step includes 
a geometrically-similar rectangular waveform modeling 
15 step of modeling the instantaneous current as a 

geometrically-similar rectangular waveform whose height 
and width are calculated such that the area of the 
rectangular waveform becomes equal to the amount of 
electric current of each event, and the FFT processing step 
20 includes a step of subjecting to FFT processing 
information concerning a change in current, the 
information being calculated in the geometrically-similar 
rectangular waveform modeling step. 

5. The electromagnetic interference analysis method 
25 as defined in claim 1, wherein the modeling step includes 
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a rectangular waveform modeling step of calculating the 
instantaneous electric current for each event information, 
and a step of modeling the instantaneous current as a 
rectangular waveform through use of the amount of electric 
5 current and a table representing the relationship between 
the width and height of a rectangular waveform, to thereby 
subject to FFT processing the information concerning a 
change in electric current calculated in the rectangular 
waveform modeling step. 

10 6. The electromagnetic interference analysis method 

as defined in claim 1, wherein the modeling step includes 
a step of calculating the instantaneous electric current 
for information for each event, and a rectangular waveform 
modeling step of modeling the instantaneous current as a 

15 rectangular waveform through use of a slew in input 

waveform and a table representing the relationship between 
the width and height of a rectangular waveform, to thereby 
subject to FFT processing the information concerning a 
change in electric current calculated in the rectangular 

20 waveform modeling step. 

7 . The electromagnetic interference analysis method 
as defined in claim 1, wherein the modeling step includes 
a step of calculating the instantaneous electric current 
for information for each event, and a rectangular waveform 

25 modeling step of modeling the instantaneous current as a 
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rectangular waveform through use of an output load 
capacitance and a table representing the relationship 
between the width and height of a rectangular waveform, 
to thereby subject to FFT processing the information 
5 concerning a change in electric current calculated in the 
rectangular waveform modeling step. 

8 . The electromagnetic interference analysis method 
as defined in claim 1, wherein the modeling step includes 
a step of calculating a drop in voltage from the amount 

10 of electric current flowing in each cell and the resistance 
of a power and correcting the amount of instantaneous 
electric current of each cell for each event, on the basis 
of the relationship between the drop in voltage and the 
amount of instantaneous electric current. 

15 9. The electromagnetic interference analysis method 

as defined in claim 1, wherein the modeling step includes 
a step of calculating a drop in voltage from the amount 
of electric current flowing in each cell, the resistance 
of a power line, and the capacitance of an on-chip 

20 capacitor, and correcting the amount of instantaneous 
electric current of each cell for each event, on the basis 
of the relationship between the drop in voltage and the 
amount of instantaneous electric current. 

10. The electromagnetic interference analysis 

25 method as defined in claim 1, wherein the modeling step 
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includes a step of transiently analyzing a power RC of each 
cell and a cell power source, accurately calculating a drop 
in voltage, and a correction step of correcting the amount 
of instantaneous electric current of each cell for each 
5 event, on the basis of the relationship between the drop 
in voltage and the amount of instantaneous electric 
current . 

11. The electromagnetic interference analysis 
method as defined in claim 1, wherein the modeling step 

10 includes a triangular waveform modeling step of modeling 
the instantaneous current as a triangular waveform which 
has a given width and whose height is calculated for each 
event information such that the amount of instantaneous 
electric current becomes equal to the area of the 

15 triangular waveform, and the FFT processing step includes 
a step of subjecting to FFT processing information 
concerning a change in current, the information being 
calculated in the triangular waveform modeling step. 

12. The electromagnetic interference analysis 

20 method as defined in claim 1, wherein the modeling step 
includes a multi-order-function waveform modeling step of 
modeling the instantaneous current as a multi-order- 
function waveform, and the FFT processing step includes 
a step of subjecting to FFT processing information 

25 concerning a change in current, the information being 
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calculated in the multi-order-function waveform modeling 
step . 

13. The electromagnetic interference analysis 
method as defined in claim 1, wherein the modeling step 

5 includes an exponential function waveform modeling step 
of modeling the instantaneous current as an 
exponential -function waveform, and the FFT processing 
step includes a step of subjecting to FFT processing 
information concerning a change in current, the 
10 information being calculated in the exponential -function 
waveform modeling step. 

14 . The electromagnetic interference analysis 
method as defined in claim 1, wherein the modeling step 
includes a step of modeling the amount of instantaneous 

15 electric current while separating the same into a short 
circuit electric current component and a charge current 
component . 

15. The electromagnetic interference analysis 
method as defined in claim 1, wherein the modeling step 

20 includes a calculation step of calculating the height of 
a rectangular waveform from a library in which peak 
currents of cells are characterized according to the type 
of cell, and a rectangular waveform modeling step of 
modeling the amount of instantaneous electric current as 

25 a rectangular waveform, the peak current calculated in the 
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calculation step being taken as the height of the 
rectangular waveform and the area of the triangular 
waveform being equal to the amount of electric current of 
each event, and the FFT processing step includes a step 
5 of subjecting to FFT processing information concerning a 
change in current, the information being calculated in the 
rectangular waveform modeling step. 

16. The electromagnetic interference analysis 
method as defined in claim 15 , wherein the calculation step 

10 includes a step of calculating a peak current for each cell 
from information concerning a slew in the cell, by 
reference to a library in which the relationship between 
a slew in input waveform and a peak current is 
characterized in the form of a table according to the type 

15 of cell. 

17. The electromagnetic interference analysis 
method as defined in claim 15 , wherein the calculation step 
includes a step of calculating a peak current for each cell 
from information concerning a load capacitance of a cell, 

20 by reference to a library in which the relationship between 
a load capacitance and a peak current is characterized in 
the form of a table according to the type of cell . 

18. The electromagnetic interference analysis 
method as defined in claim 15 , wherein the calculation step 

25 includes a step of setting a plurality of peak currents 
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for a composite cell and calculating the heights of a 
plurality of rectangular waveforms through use of a 
characterized library, and the rectangular waveform 
modeling step corresponds to a step of modeling the amount 
of electric current into a plurality of rectangular 
waveforms . 

19. The electromagnetic interference analysis 
method as defined in claim 15 , wherein the calculation step 
includes a step of setting a plurality of peak currents 
for each of the rise and fall of a flip-flop (FF) cell and 
calculating the heights of a plurality of rectangular 
waveforms through use of a characterized library, and the 
rectangular waveform modeling step corresponds to a step 
of modeling the amount of electric current into a plurality 
of rectangular waveforms . 

20. The electromagnetic interference analysis 
method as defined in claim 15 , wherein the calculation 
step includes a step of calculating the height of a 
rectangular waveform through use of a library in which peak 
currents are characterized, in consideration of the state 
of an input signal . 

21. The electromagnetic interference analysis 
method as defined in claim 15, wherein the modeling step 
includes 

a step of calculating a drop in voltage from the 
201 




amount of electric current determined according to the 
type of cell and from the resistance of a power line; and 
a correction step of characterizing, for each cell, 
the relationship between a drop in voltage and the amount 
5 of instantaneous electric current in the form of a table, 
to thereby correct the amount of instantaneous electric 
current for each event of the cell . 

22. The electromagnetic interference analysis 
method as defined in claim 15, wherein the modeling step 
10 includes 

a step of calculating a drop in voltage from the 
amount of electric current determined according to the 
type of cell, the resistance of a power line, and the 
capacitance of an on-chip capacitor; and 

15 a correction step of characterizing, for each cell, 

the relationship between a drop in voltage and the amount 
of instantaneous electric current in the form of a table, 
to thereby correct the amount of instantaneous electric 
current for each event of the cell . 

20 23. The electromagnetic interference analysis 

method as defined in Claim 10, wherein the correction step 
includes a step of iterating several times calculation of 
a drop in voltage and correction of a current waveform, 
24 . The electromagnetic interference analysis 

25 method as defined in claim 15 , wherein the calculation step 
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includes a step of modeling the amount of instantaneous 
electric current while separating the same into a short 
circuit electric current component and a charge current 
component . 

5 25. The electromagnetic interference analysis 

method as defined in claim 1, wherein the modeling step 
includes 

a triangular waveform modeling step of modeling the 
instantaneous current as a triangular waveform whose width 

10 is calculated for each event information in consideration 
of slew information (i.e., an output slew) for an output 
terminal of a cell for each event information such that 
the area of the triangular waveform becomes equal to the 
amount of electric current of each event, the height of 

15 the triangular waveform being calculated on the basis of 
the width, and 

the FFT processing step includes a step of subjecting 
to FFT processing information concerning a change in 
current, the information being calculated in the 

20 triangular waveform modeling step. 

26. The electromagnetic interference analysis 
method as defined in claim 1, wherein the modeling step 
includes 

a triangular height calculation step of calculating 
25 the height of a triangular waveform such that the area of 
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the triangular waveform becomes equal to the amount of 
electric current of each event, by means of multiplying 
the amount of instantaneous electric current by a 
coefficient corresponding to the state of an event of a 
5 cell, in consideration of whether the event of the cell 
is in a rising state or a falling state. 

27. The electromagnetic interference analysis 
method as defined in claim 1, wherein the modeling step 
includes 

10 a step of calculating the amount of instantaneous 

electric current for each event information in the case 

of a composite cell; and 

a triangular waveform modeling step of modeling the 

amount of instantaneous electric current as a plurality 
15 of triangular waveforms which are equal in number to the 

stages provided in the composite cell, through use of a 

table representing the relationship between the width and 

height of a triangular waveform; and 

the FFT processing step includes a step of subjecting 
20 to FFT processing information concerning a change in 

current, the information being calculated in the 

triangular waveform modeling step. 

28. An electromagnetic interference analysis system 
for analyzing the amount of electromagnetic interference 

25 arising in an LSI by means of performing a logic simulation , 
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the system comprising: 

a logic simulator; 

computation means which is connected to the logic 
simulator and calculates the amount of instantaneous 
electric current from event information, the information 
being produced when a change arises in a signal and 
including the instance name of each cell in which the 
change has arisen, the name of the signal, a time at which 
the change has arisen, and transition information; 

modeling means for modeling the instantaneous 
electric current according to a predetermined rule; and 
fast Fourier (FFT) conversion means for subjecting 
to fast Fourier processing the information concerning a 
change in electric current, the information being 
calculated through a modeling step , thereby analyzing the 
amount of electromagnetic interference arising in an LSI 
on the basis of a signal output from the FFT conversion 
means . 
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ABSTRACT 

In a gate-level logic simulation, a change in electric 
current is calculated from event information 5 output from 
a logic simulator 4 through use of a current waveform 
5 calculation section 7. The thus-calculated change in 
current is subjected to FFT processing through use of an 
FFT processing section 9, thereby determining a frequency 
characteristic of EMI and enabling EMI analysis. 



206 



:50376l0033 



1 1¥ 7^ 9 B 
M: 1/94 




^mm^= 5037610033 



mma ¥fi5ci 1^ 7^ 9a 

M: 2/ 94 




03 mmmT^tntamm 

m^im it [As] 



0 



5 



|:-^^ 5037610033 



mti^B 1^ 7^ 9B 
M: 3/ 94 



S1010 







lb = 0 











lb = lb + i(t) 







0 



101 1 


lave 


= Ib/T 








1012 










1013 


I(x) = 


: lave 








X = X 


+ At 



mmn ¥5^1 i¥ 7^ 9b 
M: 4/ 94 



n.mM^= 5037610033 



®5 msitiB^ifmaii^mc^jsejf^^i) 



T = 5 CD^ 



# • 

iittiB 1 1^ 7^ 9 B 
mm.^-^= 5037610033 M: 5/94 

gl6 FFTffi^fg^ 



dB 



0 




^7 ^r^ri&=E^)i^mmmm2) 




c 



t - W/2 t t + W/2 



:5037610033 



^ffiB 1^ 7^ 9B 
M: 6/ 94 








21 


















8 



5037610033 



M: 7/ 94 



119 l|l£D$gJf^;j^^i^MS:7P-li] 









h = i/W 







X = -W/2 -> W/2 







I(t+x) = I(t+x) + h 







x = x+ At 



5037610033 



^tBB ¥J5Jc1 l^P 9B 
M: 8/ 94 



W = 5 CDB# 



Fi<J- II 



111 $gM^^;i'(OTf^S3) 



i =Wxh = - 
h = V (i X K) 
h/w = K = — ^ 



K 



t-w/2 t l + w/2 



5037610033 



¥J55c 1 1 ^ 7 ^ 9 B 
M: 9/ 94 




^= 5037610033 



mmn 1^ 7^ 9B 

M: 10/ 94 



h = X K) 



w = 


i/h 






X = -W/2 W/2 






I(t+x) = I(t+x) + h 






X= X 


+ At 




mthB 1 1^ 7^ 9 0 
mm^-^= 5037610033 M: 11/ 94 























At 1 





















0 5 10 



1 



t - w/2 t t + w/2 



^Mm^= 5037610033 



M: 12/ 94 



11116 i-h^—^J'^U 




0 50 70 lOOi 



^mM^= 5037610033 



mttiB 11^ 7^ 9 H 
M: 13/ 94 





^a#^= 5037610033 



mmn 1^ 7^ 9B 

M : 14/ 94 



(il,hl)(i2,h2)tti±i 







, h2-hl V . 1 






w = i/h 







X - -W/2 W/2 



I(t+x) = I(t+x) + h 



mM.^^= 5037610033 



mmu 1^ 7^ 9B 

M: 15/ 94 




i = 70.s = i 

sl = 0,s2 = 2,il = 50,i2 = 100 
h(sl.il) = 5,h(sl,i2) = 14 
h(s2,il) = 6,h(s2.i2) = 20 

r h(sl.il)(s2-s) ^ h(s2,il)(s-sl) Y ^ r h(sl.i2)(s2-s) ^ h(s2.i2)(s-s]) V i-il ^ 
^'•'^ [ (s2-sl) (s2-sl) Jli2-ilJ^ (s2-sl) (s2-sl) Jii2-ilJ 

h(1.70) = 10.1 



^-^= 5037610033 



M: 16/ 94 



gi2i m^immi^mm 




¥fi5c 1 1 ^ 7 ^ 9 B 
M: 17/ 94 



1^= 5037610033 



g]22 m4(D^?^;^SIfi^M^n— E 



x = -W/2— W/2 ] 1056 



h(sl,il),h(sl,i2) 4^^, T053 
h(s2.il) , h(s2.i2) ' 



I(t+x) = I(t+x) + h 



^|' h(sl,il)(s2-s) ^ h(s2,il)(s-sl) Y i2-i 
"(^ (s2-s]) (s2-sl) Jii2-il, 

Yh(sl. 
t (s2- 



il.i2)(s2-s) ^ h(s2,i2)(s-sl) Y i-il > 
(s2-sl) ^ (s2-sl) Jli2-il J 



E123 m3fi-»^tBfgm(ll»Si5) 



6.9 








10.1 



5 t 



^tBB 1^ 7^ 9B 
5037610033 M: 18/ 94 



i]24i-c-h-T— ^S^-ZU 





0 


50 


100 


0 


0 


5 


20 


20 


0 


8 


25 




= 70,c = loco^-a- 



cl = 0,c2 = 20,11 = 50,52 = 100 
h(cl,il) = 5.h(cl,i2) = 20 
h(c2,il) = 8,h(c2,i2) = 25 

_/^ h(cl,il)(c2-c) ^ h(c2,il)(c 
~[ (c2-cl) (c2-cl) 
h(10,70) = 12.9 



h(s.i) = 



fa(cl,>2)(c2-c) ^ h(c2.i2Xc-c 



"IS] 



5037610033 



7^ 9 B 
[: 19/ 94 



11125 m.^Li^'mntii^^ 




I 







64 1 















mMM^= 5037610033 



^fttiB 1^ 7^ 9H 
M: 20/ 94 



(1126 m50)^jf^;SllMS:7o— g 



i-c-h-7 — 















h(c2,il) , h(c2,i2) 



r h(cl,il)(c2-c) ^ h(c2,il)(c-cl) Y i2-i ^ 

\ (c2-cl) (c2-cl) Jli2-il J 

i2)(c2-c) h(c2,i2)(c-cl) Y i-iQ 

:2-cl) (c2-cl) Ai2-ilJ 



(c2-, 



x = -W/2->W/2 ] 1066 







I(t+x) = I(t+x) + h 







C= 1 



:5037610033 



^tBB 7^ 9B 
M: 22/ 94 



Al 0 





1 2 




B1 












B6 



:5037610033 



m 9 B 

I: 23/ 94 



-f>X^?>X1 -i'>X^z>X2 -f'>X^?>X3 



v/DD- l S gl h- T Sg2 

4 



12 



13 



'OX^^T.I -r>X^>7^2 'i-^T^^T.Z 

sfio sfli sti2 



5ND[i Sg7 
R1 



Sg8 



Sg9 



<R4 <R5 < 



1= 5037610033 



^ffiS 1^ 9B 
M: 24/ 94 



Isl 



Is2 



Is3 



>Is4 



n 



Isl =11+12+13 

Is2=I2+I3 

Is3=I3 

Is4=I1 

Is5=I2 

Is6=I3 



j>Is5 |>Is6 



12 



13 



-r>X^>X1 -<>X^>X2 ^>;=5.^>X3 



AVsl 

N/dd-AAAa- 



AVs2 AVs3 
Vs4 ^AVs5 ^AVs6 



AVsl^RI *Is1 
AVs2 = R2 * Is2 
AVs3 = R3 * Is3 
AVs4 = R4*Is4 
AVs5 = R5 *Is5 
AVs6 = R6 * Is5 



n 




12 




13 



-OX^>X2 



:5037610033 



^ms 1^ 7^ 9 B 
M: 25/ 94 



AVs2 AVs3 
Vs4 >AVs5 >AVs6 



n 




12 




13 


AVI 




AV2 




AV3 



-OX-J^^XI '1'>X^>X2 -<>X^>X3 



'f>7.^>X3CDmEl^T 



AVI = AVsl + AVs4 
AV2= AVsl + AVs2+ AVs5 
AV3= AVsl + AVs2+ AVs3+ AVs6 



^^^^= 5037610033 



mniB 1^ 7^ 9B 
M: 28/ 94 



Sg2 h -T Sg3 



£g5 



NC- | Sg7 




13 



v/dd-AAAa 



AVs21 AVs31 

-AAA/ cAAA/- 



3ND-^AAA 




AVs11=f(C1. AVsl) 
AVs21 =f(C2. AVs2) 
AVs31 =f(C3. AVsS) 
AVs41 = AVs4 
A Vs51 = AVsS 
AVs61 = AVs6 



'r>X^>X2CDl|I±(^T 



AVI = AVsl 1 + AVs41 

AV2= AVsl 1 + AVs21 + AVs51 

AV3= AVsl 14- AVs21 + AVs31 + AVs61 



^mM^= 5037610033 



imS 11^ 7^ 9 B 
M: 29/ 94 




^g^-^= 5037610033 



^tiiB 1^ 7^ 9B 
M: 31/ 94 



y 




t - W/2 t t + W/2 



5037610033 



itma ¥J55c1 1^ 7^ 9 B 
M: 32/ 94 





mi 















r 


100 




101 






1 


r 



^g^^= 5 037610033 



mma 1^ 7^ 9B 

M: 33/ 94 









h = 2xi/W 







X = 0 W/2 1 1 1 02 



I(t+x) = I(t+x) - X + h 

I(t-x) = I(t-x) - -^^x + h 



m&B ¥fi£l 1^ 7^ 9B 
5037610033 M: 34/ 94 



40 










At 


20 













W = 5 CD^ 



0 5 10 




W = — ^ 

y = a (x - W?(x + Wf 

15 . 
^= I6xw^' 
4 W=2.5, j=100CDB# 
h =a X W a = 0.96 
h = 37.5 



W = 2.5 . i 
a = 0.48 
h= 18.75 



= 50 (D^ 



^^^= 50 3 7610033 



1^ 7^ 9B 
M: 35/ 94 



mmm 


m 




r 





mn.m^^= 5037610033 



1^ 7^ 9B 
M: 36/ 94 









16 xW^ 





















I(t+x)= I(t+x)+a (X -W) 2 (x+W) 2 










x = x+ At 










Re-. 




y = h X e 



i 


h 


0 


0 


50 


10 


100 


40 




i =50, 100(7? ig^ 
(il,hl)=(50,5) 
(j2,h2)=( 100,20) 



h(50)-10 
h(100)=40 



gg^-^= 5 0 3 7 6 1 0 0 3 3 



i 


w 


0 


0 


20 


1 


100 


1-5 




0 20 50 



^^11^ 7 ^ 9 B 
M: 38/ 94 



i = 50.100 0?^^ 

(il.wl)=(20,l) 

(i2.w2)=(100,1.5) 



_ W2-Wl ^; 

" i2-il 



w(50)=1.3 
w(100)=1.5 



(i-il)+wl 



mmm^= 5037610033 



jtma 1^ 7^ 9B 

M: 39/ 94 



ms3 m^dmmtamm 



1 ^ 

mm^Lum ^ ^ 


















mm 






^tHH 1^ 7^ 9H 
^MM^= 5037610033 M : 40/ 94 





11 20 








i - w -J — 


1121 






1 122 








(il,hl)(i2,h2)ttili 


1 123 








, h2-hl,. ,1 


1 124 








(il,wl)(i2,w2)ttai 


1 125 








W2-Wl , 

w- .^ .^ (i-il)+wl 


1 126 






1 1 27 


I x = -3Xw-»3Xw 






I(t+x) = I(t+x) + h X e ^ w ^ 


1 128 








[^^^^^ x = x+At 







j:^= 5037610033 



littiH 1^ 7^ 9B 
M: 42/ 94 



7306 




2iA~-^^^401 

i=3ir>^7402 



IliiiA: 



;,^7404 



5 0 3 7 6 1 0 0 3 3 



g: 43/ 94 



3 
2 
1 



7601 



Ins 

3ns 



7602 



¥J55c 1 1 ^ 7 ^ 9 a 
n : 44/ 94 



5037610033 



Fi9^ to 



7705^ 

7706^ 
7707 



7701^ 
7702 



7703W 



7704^ 



5 037610033 



7^ 9 B 
: 45/ 94 




5037610033 



Uti^B ^mil^ 7^ 9B 
M: 46/ 94 













w = i/h 







x = -W/2-W/2 lii34 







I(t+X) = 


I(t+x) + h 






x = x+ At 



:5037610033 



Fi^. fo3 



m&B ¥55cll^ 7^ 9B 
M: 47/ 94 



E]63s-hT Zfjiy 

h 



s 


h 


0 


0 


1 


5 


2 


20 




s= 1.2 COii-^ 

(sl,hl)=(l,5) 

(s2,h2)=(2.20) 



I1I64 m^ii^^UtiiBM 




^j^^^ = 5037610033 



s - h f—^ )V=y^^-y^)M3i-^^ 



m-i''^mmm}\^--:^\ 1142 







(sl,hl)(s2,h2)tttii 






h=-i^(s-sl)+hl 
s2-sl 






w = 


= i/h 







x = -W/2-W/2 ^ 







I(t+x) = I(t+x) + h 






x = x+ At 



^g#-§-= 50 3 7610033 



¥fiJcl 1^ 7M 9B 
M: 49/ 94 



mee c - h ^ 



c 


h 


0 


0 


10 


5 


20 


20 



2C- 



(cl,hl)=(10,5) 
(c2.h2)=(20,20) 

h(c)=-^f^(c-cl)+hl 



h(12)=8 



:B 037610033 





r 






151 




r 






152 




r 






153 1 










r 



^Igg:^ = 5037610033 



M: 51/ 94 









^m-j — ^ 









X = -W/2 W/2 







(cl,hl)(c2,h2)ati[S 






c2-cl 


(:c-cl)+hl 






w 


= i/h 











I(t+x) = 


I(t+x) + h 






X = > 


+ At 



g-^r= 5037610033 



1. 5 
0. 5 



Ins 2ns 



0 3 7610033 



M: 53/ 94 



F.4-. n I 



1. 5 
0. 5 



Ins 2ns 



9003 



9004- 



9002 9001 



^Mm^ = 5037610033 



M : 54/ 94 



NANDNOR 



A ^ 

B o- 



-c> Y 



1. 5 
0. 5 



Ins 2ns 



m 



JSHiB ¥fiSl 1^ 9H 
5037610033 H: 55/ 94 



NANDNOR 



9301 9302 9303 9304 




9305 
9306 
9307 



:5037610033 



iitBS 1^ 7^ 9B 
M: 56/ 94 



9404 



9405 



9401^ 
9402 



X 



R*. 77 



FF 



CK' 





> 




> 

> 


[ 





-O Q 



50 3 7610033 



CK 



D 0 



i2 
i1 



mm 



5037610033 



M: 58/ 94 



rise 15mA 7. 5mA \ 

fall 5mA ^ 2. 5mA i 



9901 9902 



5Q3761Q033 



M: 60/ 94 



10101 ^ 



10102 



10104 



10105 



1 01 Q^-^ m^^mmmm^^^^ 



NAND 




o Y 



1= 5037610033 



lims 1^ 9B 
M: 61/ 94 



A 
B 



11 



mma w-^i 7m 9 a 

M: 62/ 94 



f= 5037610033 



► 

mm ^ 



liSfng^ NAND Yi^^ 
i 'A=l, B=0^1B#cr)m5ttigM ^-^-^ SmA i 
i A= 1 , B=0-* 1 B#cD t° - ^ 1 mk ; 

T=l, B=1^0iFCDSjjfti^.M 7.^.1 3mA 
A=1, B=1->0B#(Dt°-^<g 1mA 
B=1, A=0-»1Bt(Dm}?Stli« -^^4 3mA 
B=l, A=0-^ia#(Dt°-i>fi IfflA 
B=1,A=W0B#(DS}tlSM 3mA 
B=1, A=1-^0B#(Dt°-^'(g r ^ -' 1mA 
A=0-^1, B=0->1B#C9S;y1tlfe>m ' ■ 4mA 
A=0-*1, 8=0^1 B#CDt°-^'(il A--" 1, 5mA 
A=WO, B=1-^0P#Oli>MtM '-'^ 4mA 
A=WO, B^WOP^cDt'-^'fit , - 1. 5mA 



:Kn3 7S10033 



atBH ¥h£11* 7H 9B 
W: 63/ 94 



10701^ 
10702 



v4 



X 



X 



10704 



10705 



ItHS 1^ 7^ 9B 
M : 64/ 94 



1= 5037610033 



^tbH 1 1^ 7^ 9 B 
M: 65/ 94 



f4, ^1 




msm^ = 5 037610033 

F;4- '^'2, 








AULA'S 











1^ 7^ 9B 

^g^-^=5 0 3 7610033 





50 3 7610033 



^tBB ¥fi^ll^^7M 9B 
M: 69/ 94 




:B 037610033 



M: 70/ 94 



Pie, ^6 



:0 




t - w/2 t t + w/2 



5037610033 



iiaiB 1^ 7^ 9B 
M: 71/ 94 



1197 mSfJSff^^tii^aS 










^3 






r 




mi 




i 











:5037610033 



398 ^f^}^mm^myn—m 



1 160 

















w =2 X i/h 






^^^""""^ x = 0 


w/2 



1163 



I(t+x) = I(t+x) - x + h 

w 

I(t-x) = I(t-x) - ^ x + h 



§:^= 5037610033 



^tiiB 1^ 7^ 90 
M: 73/ 94 




y = a (x - wHx + w)^ 



1:^= 5 037610033 



ItHB 1 1^ 7^ 9 H 
M : 74/ 94 



Fi0 too 




^j$^^= 5 0 3 7610033 

B9 1^1 



M: 75/ 94 



0101 m^'^m^mmmmz?o-m 



I(t + x) = I(t+x) + a(x-w)-(x + w)^ 



1171 













15i 








w-^ w 



5037610033 



1 1^ 7^ 9 B 
M: 76/ 94 



Fc^, lOZ 



m 1 02 m^i^Mnta^m 




^^g^^= 5 037610033 

R9. lOS 



iimS 1^ 7^ 90 
M: 77/ 94 



^^^^ x=-3Xw-»3Xw '"^^"'^ 






I(t + x) = I(t + 


x) + hxe"^*-^ 



JSMJH ¥S8 1 1* 7^ 9 B 
mm^9= 5037610033 . ^- ^'^ 

7U^;/y7P':/:/FF/Q 

SilS)»M'* * 3mA 

^^-^-vm^^e-^M / 4 2mA 

y u :/:7 p y ^ff/ck 

/\-^y7 7BUF/Y 

miim}«M "^"-^ 2mA 

miim>)ite-^{iiK^MmA 
^'i'-vmj/te-^'ft , - omA 



: 5 0 3 7 6 1 0 0 3 3 



¥fi5cl 1^ 9B 
M : 80/ 94 



3 

2 
1 



h nil In nil Ih Ih In 



mm 



5037610033 



7^ 9B 

M: 81/ 94 



12501^ 
12502 



12504 



12505 



i 



t.^= 5037610033 



SmS ¥fiSl 1^ 1^ 9B 
a : 82/ 94 



Fi3- {OS 
NANDNOR 



A 
B 

C 




1003 



1002 



1001 



1005 



1004- 



^= 5037610033 



fid uo 



M: 83/ 94 



~ D Q " 
-> CK 



CLK 



CLK 

D 
Q1 



D1 
Q 



mm ^! 



mtUB ¥mi i¥ IS 9B 

^mm^= 5037610033 M: 84/ 94 

<xMimm}^m 2tnA::r::: ^ 1 3oz 



^^yZTm ^^,^mz 



fie. us 



3 
2 
1 



mm 



^g^^= 50 3 7610033 



mtHE 1^ 7fl 9H 
M: 85/ 94 



1504 



1505 



1501^ 
1502 



5037610033 



Fv9- 115 



CLK 



iSttiB 1 1^ 7^ 9 B 
M: 86/ 94 



Q1 

D1 
Q 



:5Q37610033 



^iiia 1 1^ 7 H 9 0 
M: 87/ 94 





5037610033 

(a) 



¥J5jcl 1^ 7^ 9 B 
M: 88/ 94 



(b) 



5037610033 



ItBH ¥J52 1 1^ 7^ 9 0 
M: 89/94 



e\9 m 



1171 



1172 1173 



D 

m 
I- 

I 




1175 



1176 



o 

M 
O 
M 
O 
O 



/ 



1-^= 5037610033 



^ttSH ¥mi 1^ m 9H 
M: 90/ 94 



Re- lie 



1 190 















h=2 X l/w 






x = 0 


^ w/2 


1 


I(t+x) = I(t+x) 
I(t-x) = I(t-x) 


2xh ^, 

- X + ll 

w 

- X + h 

w 


1 


L^^^^^ X = X 


+ At 



11 92 



1 193 



1 194 



1 195 



1 196 



5037610033 



MffiB ^mi 1^ 7^ 9B 
M: 91/ 94 



(a) 

1 



(b) 



l;j1tMl?<2xriseBt{^it 



i/tib^XUir— w 

m^iMl X 2 X fallBtfigSf 

/ai^ixuzL— w 





5037610033 



¥J5g 1 1^ 7^ 9 B 
M: 92/ 94 



1201 



1202 



1203 



1204 



1207 



rise/falHa^#f^i^I 



1205 



1206 



I^r^ 5037610033 



1^ 7^ 9H 
M: 93/ 94 



(a) 



0 



(b) 



SMtd 



x2/t±i±iXij3.— 



/(iss-ai^xu^— i2) 




!-§•= 5037610033 



M: 94/ 



1221 



1222 



1224 



iBIt^lS 



IBII^IS 



1223 



1225 
1226 




[Description of Reference Numerals] 

1 net list 

2 test vector 

3 library 

4 logic simulator 

5 event information 

6 simulation result 

7 current waveform calculation section 

8 current waveform calculation result 

9 FFT processing section 

10 result of EMI analysis 

11 instantaneously-changing current calculation section 

12 calculation result of instantaneously-changing 
current 

13 current waveform averaging section 

20 first rectangular waveform height calculation section 

21 first rectangular waveform- shaping section 

30 second rectangular waveform height calculation 
section 

31 second rectangular waveform width calculation section 

32 second rectangular waveform-shaping section 

40 i-h table 

41 third rectangular waveform height calculation section 

42 third rectangular waveform width calculation section 

43 third rectangular waveform-shaping section 

50 i-s-h table 

51 Slew data 

52 fourth rectangular waveform height calculation 
section 

53 fourth rectangular waveform width calculation section 

54 fourth rectangular waveform- shaping section 

60 i-c-h table 

61 capacitance data 

62 fifth rectangular waveform height calculation section 

63 fifth rectangular waveform width calculation section 

64 fifth rectangular waveform- shaping section 
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100 triangular waveform height calculation section 

101 triangular waveform- shaping section 

120 Gaussian function waveform height calculation 
section 

121 i-w table 

122 Gaussian function waveform width calculation section 

123 Gaussian function waveform- shaping section 

130 peak current library 

131 sixth rectangular waveform height calculation 
section 

132 sixth rectangular waveform width calculation section 

133 sixth rectangular waveform shaping section 

140 s-h table 

141 seventh rectangular waveform height calculation 
section 

142 seventh rectangular waveform width calculation 
section 

143 seventh rectangular waveform shaping section 

150 c-h table 

151 eighth rectangular waveform height calculation 
section 

152 eighth rectangular waveform width calculation 
section 

153 eighth rectangular waveform shaping section 

160 triangular waveform height calculation section 

161 triangular waveform width calculation section 

162 triangular waveform shaping section 

170 multi -order-function waveform height calculation 
section 

171 multi-order-function waveform width calculation 
section 

172 multi-order-function waveform shaping section 

180 Gaussian function waveform height calculation 
section 

181 width library 

182 Gaussian function waveform width calculation section 
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183 Gaussian function waveform shaping section 
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[DESIGNATION OF DOCUMENT] ABSTRACT OF THE DISCLOSURE 
[ABSTRACT] 

[PROBLEM] EMI analysis involving a transistor-scale simulation 
requires consumption of much processing time, and analysis of 
overall EMI arising in an LSI has been unrealistic. Further, 
a gate-scale logic simulation cannot correctly replicate a 
change in electric current. Hence, the gate -scale logic 
simulation cannot be used for EMI analysis. 

[Means of Solution] In a gate-scale logic simulation, a change 
in electric current is calculated from event infoirmation 5 
output from a logic simulator 4 through use of a current waveform 
calculation section 7 . The thus-calculated change in current 
is subjected to FFT processing through use of an FFT processing 
section 9, thereby determining a frequency characteristic of 
EMI and enabling EMI analysis. 

[Selected Drawing] FIG. 18 
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[DESIGNATION OF DOCUMENT] 
[FIG. 1] 

BLOCK DIAGRAM SHOWING OVERALL FLOW OF PROCESSING OF EMI ANALYSIS 
METHOD ACCORDING TO THE PRESENT INVENTION 

1 NET LIST 

2 TEST VECTOR 

3 LIBRARY 

4 LOGIC SIMULATOR 

5 EVENT INFORMATION 

6 SIMULATION RESULT 

7 CURRENT WAVEFORM CALCULATION SECTION 

8 RESULT OF CURRENT WAVEFORM CALCULATION 

9 FFT PROCESSING SECTION 

10 RESULT OF EMI ANALYSIS 

[FIG. 2] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(FIRST EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

13 CURRENT WAVEFORM AVERAGING SECTION 

[FIG. 3] 

CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING ELECTRIC 
CURRENT 

AMOUNT OF CURRENT 
[FIG. 4] 

FLOWCHART FOR AVERAGING CURRENT WAVEFORM 
STEP SlOlO AVERAGING LOOP 

[FIG. 5] 

CALCULATION RESULT OF CURRENT WAVEFORM (FIRST EMBODIMENT) 
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AT T=5 



[FIG. 6] 

RESULT OF FFT ANALYSIS 

INTENSITY 

FREQUENCY 

[FIG. 7] 

RECTANGULAR WAVEFORM MODEL (SECOND EMBODIMENT) 
W IS CONSTANT 

[FIG. 8] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
{SECOND EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CBIANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

20 FIRST RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

21 FIRST RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 9] 

FLOWCHART OF FIRST RECTANGULAR WAVEFORM SHAPING OPERATION 
1020 RECTANGULAR WAVEFORM SHAPING LOOP 

[FIG. 10] 

CALCULATION RESULT OF CURRENT WAVEFORM (SECOND EMBODIMENT) 
AT W = 5 

[FIG. 11] 

RECTANGULAR WAVEFORM MODEL (THIRD EMBODIMENT) 
h/w = K = CONSTANT 

[FIG. 12] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
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(THIRD EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

30 SECOND RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

31 SECOND RECTANGULAR WAVEFORM WIDTH CALCULATION SECTION 

32 SECOND RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 13] 

FLOWCHART OF SECOND RECTANGULAR WAVEFORM SHAPING OPERATION 
1030 RECTANGULAR WAVEFORM SHAPING LOOP 

[FIG. 14] 

CALCULATION RESULT OF CURRENT WAVEFORM (THIRD EMBODIMENT) 
AT K=4 

[FIG. 15] 

RECTANGULAR WAVEFORM MODEL (FOURTH, FIFTH, AND SIXTH 
EMBODIMENTS) 

[FIG. 16] 
i-h TABLE 

IN THE CASE OF i=70 
[FIG. 17] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(FOURTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

40 i-h TABLE 

41 THIRD RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 



101 



42 THIRD RECTANGULAR WAVEFORM WIDTH CALCULATION SECTION 

43 THIRD RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 18] 

FLOWCHART OF THIRD RECTANGULAR WAVEFORM SHAPING OPERATION 

1040 READ i-h TABLE 

1041 RECTANGULAR WAVEFORM SHAPING LOOP 

1042 EXTRACT (il, hi) (i2, h2) 

[FIG. 19] 

CALCULATION RESULT OF CURRENT WAVEFORM (FOURTH EMBODIMENT) 

[FIG. 20] 
i-S-h TABLE 

IN THE CASE OF i=70, s=l 
[FIG. 21] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(FIFTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY -CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

50 i-s-h TABLE 

51 SLEW DATA 

52 FOURTH RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

53 FOURTH RECTANGULAR WAVEFORM WIDTH CALCULATION SECTION 

54 FOURTH RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 22] 

FLOWCHART OF FOURTH RECTANGULAR WAVEFORM SHAPING OPERATION 

1050 READ i-s-h TABLE 

1051 READ SLEW DATA 

1052 RECTANGULAR WAVEFORM SHAPING LOOP 

1053 EXTRACT (si, il) , (si, i2) , (s2, il) , (s2, 12) 
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[FIG. 23] 

CALCULATION RESULT OF CUKEIENT WAVEFORM (FIFTH EMBODIMENT) 
IN THE CASE OF i=70, C=10 

[FIG. 24] 
i-c-h TABLE 

[FIG. 25] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(SIXTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY -CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

60 i-c-h TABLE 

61 CAPACITANCE DATA 

62 FIFTH RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

63 FIFTH RECTANGULAR WAVEFORM WIDTH CALCULATION SECTION 

64 FIFTH RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 26] 

FLOWCHART OF FIFTH RECTANGULAR WAVEFORM SHAPING OPERATION 

1060 READ i-c-h TABLE 

1061 READ CAPACITANCE DATA 

1062 RECTANGULAR WAVEFORM SHAPING LOOP 

1063 EXTRACT 

[FIG. 27] 

CALCULATION RESULT OF CURRENT WAVEFORM (SIXTH EMBODIMENT) 

[FIG. 28] 

1 NET LIST 

2 TEST VECTOR 

3 LIBRARY 
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4 LOGIC SIMULATOR 

5 EVENT INFORMATION 

6 SIMULATION RESULT 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

9 FFT PROCESSING SECTION 

10 RESULT OF EMI ANALYSIS 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

13 CURRENT WAVEFORM AVERAGING SECTION 

A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
A9 LAYOUT DATA 

AlO I-V TABLE OR I-V FUNCTION 

All LINE VOLTAGE DROP CALCULATION/ INSTANTANEOUSLY-CHANGING 
CURRENT CORRECTION SECTION 

[FIG. 29] 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
A9 LAYOUT DATA 

AlO I-V TABLE OR I-V FUNCTION 

Bl INFORMATION CONCERNING INSTANTANEOUSLY-CHANGING CURRENT 

FOR EACH INSTANCE 

B2 EXTRACTION OF POWER NET 

B3 CALCULATION OF CURRENT FLOWING IN EACH SEGMENT 

B4 INFORMATION CONCERNING EACH SEGMENT IN POWER LINE 

B5 INFORMATION CONCERNING CURRENT PER EACH SEGMENT 

B6 CALCULATION OF RESISTANCE OF EACH SEGMENT 

B7 CALCULATION OF VOLTAGE DROP IN EACH SEGMENT 

B8 INFORMATION CONCERNING RESISTANCE OF EACH SEGMENT 

B9 INFORMATION CONCERNING VOLTAGE DROP IN EACH SEGMENT 

BIO CALCULATION OF VOLTAGE DROP OF EACH INSTANCE 

Bll INFORMATION CONCERNING VOLTAGE DROP IN EACH INSTANCE 

B12 CORRECTION OF ELECTRIC CURRENT 
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[FIG. 30] 

INSTANCE 1, INSTANCE 2, INSTANCE 3 
[FIG. 31] 

INSTANCE 1, INSTANCE 2, INSTANCE 3 
[FIG. 32] 

INSTANCE 1, INSTANCE 2, INSTANCE 3 
[FIG. 33] 

INSTANCE 1, INSTANCE 2, INSTANCE 3 
[FIG. 34] 

INSTANCE 1, INSTANCE 1, INSTANCE 3 
[FIG. 35] 

INSTANCE 1, INSTANCE 2, INSTANCE 3 
VOLTAGE DROP IN INSTANCE 1 
VOLTAGE DROP IN INSTANCE 2 
VOLTAGE DROP IN INSTANCE 3 

[FIG. 36] 

1 NET LIST 

2 TEST VECTOR 

3 LIBRARY 

4 LOGIC SIMULATOR 

5 EVENT INFORMATION 

6 SIMULATION RESULT 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

9 FFT PROCESSING SECTION 

10 RESULT OF EMI ANALYSIS 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

13 CURRENT WAVEFORM AVERAGING SECTION 

A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
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A9 LAYOUT DATA 

AlO I-V TABLE OR I-V FUNCTION 
CI C-V TABLE OR C-V FUNCTION 

C2 LINE VOLTAGE DROP CALCULATION/INSTANTANEOUSLY-CHANGING 
CURRENT CORRECTION SECTION 

[FIG. 37] 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
A9 LAYOUT DATA 

AlO I-V TABLE OR I-V FUNCTION 
CI C-V TABLE OR C-V FUNCTION 

Bl INFORMATION CONCERNING INSTANTANEOUSLY-CHANGING CURRENT 

FOR EACH INSTANCE 

B2 EXTRACTION OF POWER NET 

B3 CALCULATION OF CURRENT FLOWING IN EACH SEQUENT 

B4 INFORMATION CONCERNING EACH SEGMENT IN POWER LINE 

B5 INFORMATION CONCERNING CURRENT PER EACH SEGMENT 

B6 CALCULATION OF RESISTANCE OF EACH SEGMENT 

B8 INFORMATION CONCERNING RESISTANCE OF EACH SEGMENT 

B9 INFORMATION CONCERNING VOLTAGE DROP IN EACH SEGMENT 

BIO CALCULATION OF VOLTAGE DROP OF EACH INSTANCE 

Bll INFORMATION CONCERNING VOLTAGE DROP IN EACH INSTANCE 

B12 CORRECTION OF ELECTRIC CURRENT 

Dl CALCULATION OF CAPACITANCE PROVIDED BETWEEN SEGMENTS 
D2 INFORMATION CONCERNING CAPACITANCE PROVIDED BETWEEN 
SEGMENTS 

D3 CALCULATION OF VOLTAGE DROP IN EACH SEGMENT 
[FIG. 38] 

INSTANCE 1, INSTANCE 2, INSTANCE 3 
[FIG. 39] 

INSTANCE 1, INSTANCE 2, INSTANCE 3 
VOLTAGE DROP IN INSTANCE 1 
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VOLTAGE DROP IN INSTANCE 2 
VOLTAGE DROP IN INSTANCE 3 

[FIG. 40] 

1 NET LIST 

2 TEST VECTOR 

3 LIBRARY 

4 LOGIC SIMULATOR 

5 EVENT INFORMATION 

6 SIMULATION RESULT 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

9 FFT PROCESSING SECTION 

10 RESULT OF EMI ANALYSIS 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

13 CURRENT WAVEFORM AVERAGING SECTION 

A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
A9 LAYOUT DATA 

AlO I-V TABLE OR I-V FUNCTION 

El LINE VOLTAGE DROP CALCULATION/ INSTANTANEOUSLY-CHANGING 
CURRENT CORRECTION SECTION 

[FIG. 41] 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
A9 LAYOUT DATA 

AlO I-V TABLE OR I-V FUNCTION 

Bl INFORMATION CONCERNING INSTANTANEOUSLY-CHANGING CURRENT 

FOR EACH INSTANCE 

B2 EXTRACTION OF POWER NET 

B4 INFORMATION CONCERNING EACH SEGMENT IN POWER LINE 
B6 CALCULATION OF RESISTANCE OF EACH SEGMENT 
B8 INFORMATION CONCERNING RESISTANCE OF EACH SEGMENT 
Bll INFORMATION CONCERNING VOLTAGE DROP IN EACH INSTANCE 
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B12 CORRECTION OF ELECTRIC CURRENT 

Dl CALCULATION OF CAPACITANCE PROVIDED BETWEEN SEGMENTS 
D2 INFORMATION CONCERNING CAPACITANCE PROVIDED BETWEEN 
SEGMENTS 

Fl TRANSIENT ANALYSIS 
[FIG. 42] 

TRIANGULAR WAVEFORM MODEL (TENTH EMBODIMENT) 
W IS CONSTANT 

[FIG. 43] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(TENTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

100 TRIANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

101 TRIANGULAR WAVEFORM SHAPING SECTION 

[FIG. 44] 

FLOWCHART OF FIFTH TRIANGULAR WAVEFORM SHAPING OPERATION 
1100 TRIANGULAR WAVEFORM SHAPING LOOP 

[FIG. 45] 

CALCULATION RESULT OF CURRENT WAVEFORM (TENTH EMBODIMENT) 
AT W=5 

[FIG. 46] 

MULTI -ORDER-FUNCTION MODEL (ELEVENTH EMBODIMENT) 
W IS CONSTANT 

IN THE CASE OF W=2 . 5 , i=100 
IN THE CASE OF W=2 . 5 , i=50 

[FIG. 47] 
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BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(ELEVENTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

110 MULTI -ORDER-FUNCTION WAVEFORM HEIGHT CALCULATION SECTION 

111 MULTI -ORDER- FUNCTION WAVEFORM SHAPING SECTION 

[FIG. 48] 

FLOWCHART OF MULTI-ORDER- FUNCTION WAVEFORM SHAPING OPERATION 
1110 MULTI -ORDER-FUNCTION WAVEFORM SHAPING LOOP 

[FIG. 49] 

CALCULATION RESULT OF CURRENT WAVEFORM (ELEVENTH EMBODIMENT) 
AT W=2.5 

[FIG. 50] 

GAUSSIAN FUNCTION MODEL (TWELFTH AND TWENTY-FIFTH EMBODIMENTS) 

[FIG. 51] 
i-h TABLE 

IN THE CASE OF i=50, 100 

[FIG. 52] 
i-w TABLE 

IN THE CASE OF i=50, 100 
[FIG. 53] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(TWELFTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 
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12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
40 i-h TABLE 

121 i-w TABLE 

120 GAUSSIAN FUNCTION WAVEFORM HEIGHT CALCULATION SECTION 

122 GAUSSIAN FUNCTION WAVEFORM WIDTH CALCULATION SECTION 

123 GAUSSIAN FUNCTION WAVEFORM SHAPING SECTION 

[FIG. 54] 

FLOWCHART OF GAUSSIAN FUNCTION WAVEFORM SHAPING 

1120 READ i-h TABLE 

1121 READ i-w TABLE 

1122 PERFORM GAUSSIAN FUNCTION WAVEFORM SHAPING LOOP 

1123 EXTRACT 
1125 EXTRACT 

[FIG. 55] 

CALCULATION RESULT OF CURRENT WAVEFORM (TWELFTH EMBODIMENT) 
[FIG. 56] 

7301 CIRCUIT CONNECTION INFORMATION STORAGE MEANS 

7302 PATTERN-OF-SIGNAL-CHANGE STORAGE MEANS 

7303 ELEMENT CURRENT STORAGE MEANS 

7304 TOTAL CURRENT STORAGE MEANS 

7 605 TOTAL CURRENT CALCULATION MEANS 

7306 MEANS FOR STORING BASIC INFORMATION FOR ESTIMATING 
ELECTRIC CURRENT 

[FIG. 57] 
FLIP-FLIP FF 

THE AMOUNT OF SHORT CIRCUIT CURRENT FLOWING WHEN Y IS CHANGED 
THE AMOUNT OF CHARGE CURRENT FLOWING WHEN Y IS CHANGED 
THE AMOUNT OF SHORT CIRCUIT CURRENT FLOWING WHEN CK IS CHANGED 
THE AMOUNT OF CHARGE CURRENT FLOWING WHEN CK IS CHANGED 

BUFFER BUF 

THE AMOUNT OF SHORT CIRCUIT CURRENT FLOWING WHEN Y IS CHANGED 
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THE AMOUNT OF CHARGE CURRENT FLOWING WHEN Y IS CHANGED 

[FIG. 58] 

CURRENT 

TIME 

[FIG. 59] 

7601 DISCRETE WIDTH USED FOR AVERAGING SHORT CIRCUIT CURRENT 
COMPONENT 

7602 DISCRETE WIDTH USED FOR AVERAGING CHARGE CURRENT COMPONENT 
[FIG. 60] 

START ACTIVATION OF TOTAL CURRENT CALCULATION MEANS 

7701 READ CIRCUIT CONNECTION INFORMATION 

7702 READ PATTERN OF SIGNAL CHANGE 

7703 READ INFORMATION CONCERNING THE AMOUNT OF ELECTRIC CURRENT 
FLOWING IN ELEMENT 

7704 READ BASIC INFORMATION FOR ESTIMATING ELECTRIC CURRENT 

7705 IMPART A PATTER OF SIGNAL CHANGE TO CIRCUIT CONNECTION 
INFORMATION, TO THEREBY ENABLE PROPAGATION OF A SIGNAL. IN A 
CASE WHERE A CHANGE HAS ARISEN IN A SIGNAL APPEARING AT AN 
EXTERNAL INPUT TERMINAL OF A CIRCUIT ELEMENT, THE ELEMENT BEING 
DESCRIBED IN THE INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING THROUGH AN ELEMENT, ADD A SHORT CIRCUIT CURRENT 
COMPONENT AND A CHARGE CURRENT COMPONENT INCLUDED IN THE 
ELECTRIC CURRENT FLOWING IN THE EXTERNAL TERMINAL, AT 
RESPECTIVE POINTS IN TIME, TO THE SHORT CIRCUIT CURRENT 
COMPONENT AND THE CHARGE CURRENT COMPONENT INCLUDED IN THE 
INFORMATION CONCERNING THE TOTAL AMOUNT OF ELECTRIC CURRENT, 
RESPECTIVELY. 

7706 AVERAGE SHORT CIRCUIT CURRENT COMPONENT INCLUDED IN TOTAL 
CURRENT INFORMATION BY DISCRETE WIDTH THEREOF, AND AVERAGE 
CHARGE CURRENT COMPONENT INCLUDED IN THE SAME BY DISCRETE WIDTH 
THEREOF, ADD AVERAGING RESULTS TOGETHER AS TOTAL AMOUNT OF 
CURRENT AT RESPECTIVE TIME 

7707 STORE THE TOTAL AMOUNT OF CURRENT INTO TOTAL CURRENT 
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STORAGE MEANS 

END OF OPERATION OF TOTAL CURRENT STORAGE MEANS 
[FIG. 61] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(FOURTEENTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

130 PEAK CURRENT LIBRARY 

131 SIXTH RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

132 SIXTH RECTANGULAR WAVEFORM WIDTH CALCULATION SECTION 

133 SIXTH RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 62] 

FLOWCHART OF SIXTH RECTANGULAR WAVEFORM SHAPING 

1130 READ PEAK CURRENT LIBRARY 

1131 PERFORM RECTANGULAR WAVEFORM SHAPING LOOP 

1132 EXTRACT HEIGHT OF RECTANGULAR WAVEFORM OF OBJECT INSTANCE 
FROM PEAK CURRENT LIBRARY 

[FIG. 63] 
s-h TABLE 

IN THE CASE OF s=l . 2 
[FIG. 64] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(FIFTEENTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
51 SLEW DATA 
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140 s-h TABLE LIBRARY 

141 SEVENTH RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

142 SEVENTH RECTANGULAR WAVEFORM WIDTH CALCULATION SECTION 

143 SEVENTH RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 65] 

FLOWCHART OF SEVENTH RECTANGULAR WAVEFORM SHAPING 

1140 READ s-h TABLE 

1141 READ SLEW DATA 

1142 PERFORM RECTANGULAR WAVEFORM SHAPING LOOP 

1143 EXTRACT 

[FIG. 66] 
c-h TABLE 

IN THE CASE OF c=1.2 
[FIG. 67] 

BLOCK DIAGRAM SHOWING A CURRENT WAVEFORM CALCULATION SECTION 
(SIXTEENTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
61 CAPACITANCE DATA 

150 c-h TABLE LIBRARY 

151 EIGHTH RECTANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

152 EIGHTH RECTANGULAR WAVEFORM WIDTH CALCULATION SECTION 

153 EIGHTH RECTANGULAR WAVEFORM SHAPING SECTION 

[FIG. 68] 

FLOWCHART OF EIGHTH RECTANGULAR WAVEFORM SHAPING 

1150 READ c-h TABLE 

1151 READ CAPACITANCE DATA 

1152 PERFORM RECTANGULAR WAVEFORM SHAPING LOOP 
1143 EXTRACT 
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[FIG. 69] 
TIME 

[FIG. 70] 

CURRENT 

TIME 

[FIG. 71] 

CURRENT 

TIME 

[FIG. 72] 

9001 CIRCUIT CONNECTION INFORMATION STORAGE MEANS 

9002 PATTERN-OF-SIGNAIi-CHANGE STORAGE MEANS 

9003 ELEMENT CURRENT STORAGE MEANS 

9004 TOTAL CURRENT STORAGE MEANS 

9005 TOTAL CURRENT CALCULATION MEANS 

[FIG. 74] 
CURRENT 
TIME 

[FIG. 75] 
CHANGED PINS 
TIME 

TOTAL CURRENT 
PEAK CURRENT 

[FIG. 76] 

START ACTIVATION OF TOTAL CURRENT CALCULATION MEANS 

9401 READ CIRCUIT CONNECTION INFORMATION 

9402 READ PATTERN OF SIGNAL CHANGE 

94 03 READ INFORMATION CONCERNING THE AMOUNT OF ELECTRIC CURRENT 
FLOWING IN ELEMENT 

9404 IMPART A PATTER OF SIGNAL CHANGE TO CIRCUIT CONNECTION 
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INFORMATION, TO THEREBY ENABLE PROPAGATION OF A SIGNAL. IN A 
CASE WHERE A CHANGE HAS ARISEN IN A SIGNAL APPEARING AT AN 
EXTERNAL INPUT TERMINAL OF A CIRCUIT ELEMENT, THE ELEMENT BEING 
DESCRIBED IN THE INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING THROUGH AN ELEMENT, A RECTANGULAR WAVEFORM OF 
ELECTRIC CURRENT IS PRODUCED, BY MEANS OF AVERAGING THE AMOUNT 
OF ELECTRIC CURRENT FLOWING THROUGH THE EXTERNAL INPUT TERMINAL 
BY A PEAK VALUE. THE RECTANGULAR WAVEFORM IS ADDED TO THE 
INFORMATION CONCERNING THE TOTAL AMOUNT OF ELECTRIC CURRENT 
OBTAINED AT A TIME WHICH LAGS A PREDETERMINED PERIOD OF TIME 
FROM THE TIME A CHANGE HAS ARISEN IN THE EXTERNAL INPUT TERMINAL 
9405 STORE, INTO TOTAL CURRENT STORAGE MEANS, TOTAL AMOUNT OF 
CURRENT OBTAINED AT EACH TIME 

END OF OPERATION OF TOTAL CURRENT STORAGE MEANS 

[FIG. 78] 
LOGIC VALUE 
TIME 

[FIG. 79] 
CURRENT 
TIME 

[FIG. 80] 

CURRENT 

TIME 

[FIG. 81] 
EDGE 
CURRENT 
PEAK VALUE 

[FIG. 82] 

CURRENT 

TIME 
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[FIG. 83] 

START ACTIVATION OF TOTAL CURRENT CALCULATION MEANS 

10101 READ CIRCUIT CONNECTION INFORMATION 

10102 READ PATTERN OF SIGNAL CHANGE 

10103 READ INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING IN ELEMENT 

10104 IMPART A PATTER OF SIGNAL CHANGE TO CIRCUIT CONNECTION 
INFORMATION, TO THEREBY ENABLE PROPAGATION OF A SIGNAL. IN A 
CASE WHERE A CHANGE HAS ARISEN IN A SIGNAL APPEARING AT AN 
EXTERNAL INPUT TERMINAL OF A CIRCUIT ELEMENT, THE ELEMENT BEING 
DESCRIBED IN THE INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING THROUGH AN ELEMENT, A RECTANGULAR WAVEFORM OF 
AN ELECTRIC CURRENT IS PRODUCED BY MEANS OF AVERAGING THE AMOUNT 
OF ELECTRIC CURRENT FLOWING THROUGH THE EXTERNAL INPUT TERMINAL 
AT A RISE AND FALL OF THE SIGNAL BY MEANS OF A PEAK VALUE . THE 
RECTANGULAR WAVEFORMS ARE ADDED TO THE INFORMATION CONCERNING 
THE TOTAL AMOUNT OF ELECTRIC CURRENT OBTAINED AT A TIME WHEN 
THE CHANGE HAS ARISEN IN THE SIGNAL APPEARING AT THE EXTERNAL 
INPUT TERMINAL. 

10105 STORE, INTO TOTAL CURRENT STORAGE MEANS, TOTAL AMOUNT OF 
CURRENT OBTAINED AT EACH TIME 

END OF OPERATION OF TOTAL CURRENT STORAGE MEANS 

[FIG. 85] 
TIME 

[FIG. 86] 
CURRENT 
TIME 

[FIG. 87] 

CURRENT 

TIME 

[FIG. 88] 
NOR, Y- TERMINAL 
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TOTAL CURRENT AT LOGIC VALUE OF 1 
PEAK CURRENT AT LOGIC VALUE OF 1 
TOTAL CURRENT AT LOGIC VALUE OF 0 
PEAK CURRENT AT LOGIC VALUE OF 0 
TOTAL CURRENT AT LOGIC VALUE OF 1 
PEAK CURRENT AT LOGIC VALUE OF 1 
TOTAL CURRENT AT LOGIC VALUE OF 0 
PEAK CURRENT AT LOGIC VALUE OF 0 
TOTAL CURRENT AT LOGIC VALUE OF 1 
PEAK CURRENT AT LOGIC VALUE OF 1 
TOTAL CURRENT AT LOGIC VALUE OF 0 
PEAK CURRENT AT LOGIC VALUE OF 0 

[FIG. 89] 

START ACTIVATION OF TOTAL CURRENT CALCULATION MEANS 

10701 READ CIRCUIT CONNECTION INFORMATION 

10702 READ PATTERN OF SIGNAL CHANGE 

10703 READ INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING IN ELEMENT 

10704 IMPART A PATTER OF SIGNAL CHANGE TO CIRCUIT CONNECTION 
INFORMATION, TO THEREBY ENABLE PROPAGATION OF A SIGNAL. IN A 
CASE WHERE A CHANGE HAS ARISEN IN A SIGNAL APPEARING AT AN 
EXTERNAL INPUT TERMINAL OF A CIRCUIT ELEMENT, THE ELEMENT BEING 
DESCRIBED IN THE INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING THROUGH AN ELEMENT, A RECTANGULAR WAVEFORM OF 
AN ELECTRIC CURRENT IS PRODUCED BY MEANS OF AVERAGING THE AMOUNT 
OF ELECTRIC CURRENT FLOWING THROUGH THE EXTERNAL INPUT TERMINAL 
BY MEANS OF A PEAK VALUE IN ACCORDANCE WITH THE STATUS OF ANOTHER 
TERMINAL . THE RECTANGULAR WAVEFORM IS ADDED TO THE INFORMATION 
CONCERNING THE TOTAL AMOUNT OF ELECTRIC CURRENT OBTAINED AT A 
TIME WHEN THE CHANGE HAS ARISEN IN THE SIGNAL APPEARING AT THE 
EXTERNAL INPUT TERMINAL. 

10705 STORE, INTO TOTAL CURRENT STORAGE MEANS, TOTAL AMOUNT OF 
CURRENT OBTAINED AT EACH TIME 

END OF OPERATION OF TOTAL CURRENT STORAGE MEANS 
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[FIG. 90] 

CURRENT 

TIME 

[FIG. 91] 

1 NET LIST 

2 TEST VECTOR 

3 LIBRARY 

4 LOGIC SIMULATOR 

5 EVENT INFORMATION 

6 SIMULATION RESULT 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

9 FFT PROCESSING SECTION 

10 RESULT OF EMI ANALYSIS 

11 INSTANTANEOUSLY -CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

13 CURRENT WAVEFORM AVERAGING SECTION 

A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 

A9 LAYOUT DATA 

G2 CHARACTERIZATION 

G3 I-V TABLE OR I-V FUNCTION FOR EACH CELL 

G4 LINE VOLTAGE DROP CALCULATION/INSTANTANEOUSLY-CHANGING 
CURRENT CORRECTION SECTION 

[FIG. 92] 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
A6 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
A9 LAYOUT DATA 

Bl INFORMATION CONCERNING INSTANTANEOUSLY-CHANGING CURRENT 

FOR EACH INSTANCE 

B2 EXTRACTION OF POWER NET 

B3 CALCULATION OF ELECTRIC CURRENT FLOWING IN EACH SEGMENT 

B4 INFORMATION CONCERNING EACH SEGMENT IN POWER LINE 

B5 INFORMATION CONCERNING ELECTRIC CURRENT OF EACH SEGMENT 
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B6 CALCULATION OF RESISTANCE OF EACH SE ORIENT 

B7 CALCULATION OF VOLTAGE DROP IN EACH SEGMENT 

B8 INFORMATION CONCERNING RESISTANCE OF EACH SEGMENT 

B9 INFORMATION CONCERNING VOLTAGE DROP IN EACH SEGMENT 

BIO CALCULATION OF VOLTAGE DROP IN EACH INSTANCE 

Bll INFORMATION CONCERNING VOLTAGE DROP IN EACH INSTANCE 

HI CORRECTION OF ELECTRIC CURRENT 

G3 I-V TABLE OR I-V FUNCTION FOR EACH CELL 

[FIG. 93] 

1 NET LIST 

2 TEST VECTOR 

3 LIBRARY 

4 LOGIC SIMULATOR 

5 EVENT INFORMATION 

6 SIMULATION RESULT 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

9 FFT PROCESSING SECTION 

10 RESULT OF EMI ANALYSIS 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 

13 CURRENT WAVEFORM AVERAGING SECTION 
A9 LAYOUT DATA 

11 CHARACTERIZATION 

12 I-C-H TABLE OR I-C-H FUNCTION FOR EACH CELL 

13 LINE VOLTAGE DROP CALCULATION/INSTANTANEOUSLY-CHANGING 
CURRENT CORRECTION SECTION 

[FIG. 94] 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
A6 CALCULATION RESULT OF INSTANTANEOUSLY -CHANGING CURRENT 
TAKING INTO CONSIDERATION A DROP IN LINE VOLTAGE 
A9 LAYOUT DATA 

Bl INFORMATION CONCERNING INSTANTANEOUSLY-CHANGING CURRENT 

FOR EACH INSTANCE 

B2 EXTRACTION OF POWER NET 
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B3 CALCULATION OF ELECTRIC CURRENT FLOWING IN EACH SEGMENT 

B4 INFORMATION CONCERNING EACH SEGMENT IN POWER LINE 

B5 INFORMATION CONCERNING ELECTRIC CURRENT OF EACH SEGMENT 

B6 CALCULATION OF RESISTANCE OF EACH SEGMENT 

B8 INFORMATION CONCERNING RESISTANCE OF EACH SEGMENT 

B9 INFORMATION CONCERNING VOLTAGE DROP IN EACH SEGMENT 

BIO CALCULATION OF VOLTAGE DROP IN EACH INSTANCE 

Ell INFORMATION CONCERNING VOLTAGE DROP IN EACH INSTANCE 

Jl CORRECTION OF ELECTRIC CURRENT 

12 I-C-V TABLE OR I-C-V FUNCTION FOR EACH CELL 

Dl CALCULATION OF CAPACITANCE PROVIDED BETWEEN SEC3XIENTS 

D2 INFORMATION CONCERNING A CAPACITANCE PROVIDED BETWEEN 

SEGMENTS 

[FIG. 95] 

1 NET LIST 

2 TEST VECTOR 

3 LIBRARY 

4 LOGIC SIMULATOR 

5 EVENT INFORMATION 

6 SIMULATION RESULT 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

9 FFT PROCESSING SECTION 

10 RESULT OF EMI ANALYSIS 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OP INSTANTANEOUSLY-CHANGING CURRENT 

13 CURRENT WAVEFORM AVERAGING SECTION 
A9 LAYOUT DATA 

AlO I-V TABLE OR I-V FUNCTION FOR EACH CELL 

All LINE VOLTAGE DROP CALCULATION/INS TANTANEOUSLY-CPIANGING 
CURRENT CORRECTION SECTION 
Kl COMPARISON /DETERMINATION 

[FIG. 96] 

TRIANGULAR WAVEFORM MODEL (TWENTY-THIRD EMBODIMENT) 
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[FIG. 97] 

BLOCK DIAGRAM SHOWING CALCULATION OF A CURRENT WAVEFORM 
(TWENTY-THIRD EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
130 PEAK CURRENT LIBRARY 

160 TRIANGULAR WAVEFORM HEIGHT CALCULATION SECTION 

161 TRIANGULAR WAVEFORM WIDTH CALCULATION SECTION 

162 TRIANGULAR WAVEFORM SHAPING SECTION 

[FIG. 98] 

FLOWCHART OF TRIANGULAR WAVEFORM SHAPING 
1150 READ PEAK CURRENT LIBRARY 

1161 PERFORM TRIANGULAR WAVEFORM SHAPING LOOP 

1162 EXTRACT HEIGHT OF TRIANGULAR WAVEFORM OF OBJECT INSTANCE 
FROM PEAK CURRENT LIBRARY 

[FIG. 99] 

MULTI -ORDER-FUNCTION MODEL (TWENTY -FOURTH EMBODIMENT) 
[FIG. 100] 

BLOCK DIAGRAM SHOWING CALCULATION OF A CURRENT WAVEFORM 
(TWENTY -FOURTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
130 PEAK CURRENT LIBRARY 

170 MULTI -ORDER- FUNCTION WAVEFORM HEIGHT CALCULATION SECTION 

171 MULTI -ORDER-FUNCTION WAVEFORM WIDTH CALCULATION SECTION 

172 MULTI -ORDER- FUNCTION WAVEFORM SHAPING SECTION 
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[FIG. 101] 

FLOWCHART OF MULTI -ORDER- FUNCTION WAVEFORM SHAPING 

1170 READ PEAK CURRENT LIBRARY 

1171 PERFORM MULTI -ORDER- FUNCTION WAVEFORM SHAPING LOOP 

1172 EXTRACT HEIGHT OF TRIANGULAR WAVEFORM OF OBJECT INSTANCE 
FROM PEAK CURRENT LIBRARY 

[FIG. 102] 

BLOCK DIAGRAM SHOWING CALCULATION OF A CURRENT WAVEFORM 
(TWENTY-FIFTH EMBODIMENT) 
3 LIBRARY 

5 EVENT INFORMATION 

8 CALCULATION RESULT OF CURRENT WAVEFORM 

11 INSTANTANEOUSLY-CHANGING CURRENT CALCULATION SECTION 

12 CALCULATION RESULT OF INSTANTANEOUSLY-CHANGING CURRENT 
130 PEAK CURRENT LIBRARY 

181 WIDTH LIBRARY 

180 GAUSSIAN FUNCTION WAVEFORM HEIGHT CALCULATION SECTION 

182 GAUSSIAN FUNCTION WAVEFORM WIDTH CALCULATION SECTION 

183 GAUSSIAN FUNCTION WAVEFORM SHAPING SECTION 

[FIG. 103] 

FLOWCHART OF GAUSSIAN FUNCTION WAVEFORM SHAPING 

1180 READ PEAK CURRENT LIBRARY 

1181 READ WIDTH LIBRARY 

1182 PERFORM GAUSSIAN FUNCTION WAVEFORM SHAPING LOOP 

1183 EXTRACT HEIGHT OF TRIANGULAR WAVEFORM OF OBJECT INSTANCE 
FROM PEAK CURRENT LIBRARY 

1184 EXTRACT WIDTH OF OBJECT INSTANCE FROM WIDTH LIBARRY 

[FIG. 104] 

CURRENT 

TIME 

[FIG. 105] 
FLIP-FLOP FF/Q 
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TOTAL SHORT CIRCUIT CURRENT 

PEAK VALUE OF SHORT CIRCUIT CURRENT 

TOTAL CHARGE CURRENT 

PEAK VALUE OF CHARGE CURRENT 

FLIP-FLOP FF/CK 

TOTAL SHORT CIRCUIT CURRENT 

PEAK VALUE OF SHORT CIRCUIT CURRENT 

TOTAL CHARGE CURRENT 

PEAK VALUE OF CHARGE CURRENT 

BUFFER BUF/Y 

TOTAL SHORT CIRCUIT CURRENT 

PEAK VALUE OF SHORT CIRCUIT CURRENT 

TOTAL CHARGE CURRENT 

PEAK VALUE OF CHARGE CURRENT 

[FIG. 106] 
CURRENT 
TIME 

[FIG. 107] 

START ACTIVATION OF TOTAL CURRENT CALCULATION MEANS 

12501 READ CIRCUIT CONNECTION INFORMATION 

12502 READ PATTERN OF SIGNAL CHANGE 

12503 READ INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING IN ELEMENT 

12504 IMPART A PATTER OF SIGNAL CHANGE TO CIRCUIT CONNECTION 
INFORMATION, TO THEREBY ENABLE PROPAGATION OF A SIGNAL. IN A 
CASE WHERE A CHANGE HAS ARISEN IN A SIGNAL APPEARING AT AN 
EXTERNAL INPUT TERMINAL OF A CIRCUIT ELEMENT, THE ELEMENT BEING 
DESCRIBED IN THE INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING THROUGH AN ELEMENT, A RECTANGULAR WAVEFORM OF 
AN ELECTRIC CURRENT IS PRODUCED BY MEANS OF AVERAGING SHORT 
CIRCUIT CURRENT OF THE ELECTRIC CURRENT BY MEANS OF A PEAK VALUE 
THEREOF. A RECTANGULAR WAVEFORM OF AN ELECTRIC CURRENT IS 
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PRODUCED BY MEANS OF AVERAGING CHARGE CURRENT OF THE ELECTRIC 
CURRENT BY MEANS OF A PEAK CURRENT THEREOF. THE RECTANGULAR 
WAVEFORMS ARE ADDED TO THE INFORMATION CONCERNING THE TOTAL 
AMOUNT OF ELECTRIC CURRENT OBTAINED AT A TIME WHEN THE CHANGE 
HAS ARISEN IN THE SIGNAL APPEARING AT THE EXTERNAL INPUT 
TERMINAL . 

12505 STORE , INTO TOTAL CURRENT STORAGE MEANS , TOTAL AMOUNT OF 
CURRENT OBTAINED AT EACH TIME 

END OF OPERATION OF TOTAL CURRENT STORAGE MEANS 
[FIG. 109] 

1001 CIRCUIT CONNECTION INFORMATION STORAGE MEANS 

1002 PATTERN-OF-SIGNAL-CHANGE STORAGE MEANS 

1003 ELEMENT CURRENT STORAGE MEANS 

1004 TOTAL CURRENT STORAGE MEANS 

1005 TOTAL CURRENT CALCULATION MEANS 

[FIG. Ill] 
TIME 

[FIG. 112] 
FLIP-FLOP FF 

THE AMOUNT OF CURRENT FLOWING WHEN Y IS CHANGED 
THE AMOUNT OF CURRENT FLOWING WHEN CK IS CHANGED 

BUFFER BUF 

THE AMOUNT OF CURRENT FLOWING WHEN Y IS CHANGED 

[FIG. 113] 
TOTAL CURRENT 
TIME 

[FIG. 114] 

START ACTIVATION OF TOTAL CURRENT CALCULATION MEANS 

1501 READ CIRCUIT CONNECTION INFORMATION 

1502 READ PATTERN OF SIGNAL CHANGE 
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1503 READ INFORMATION CONCERNING THE AMOUNT OF ELECTRIC CURRENT 
FLOWING IN ELEMENT 

1504 IMPART A PATTER OF SIGNAL CHANGE TO CIRCUIT CONNECTION 
INFORMATION, TO THEREBY ENABLE PROPAGATION OF A SIGNAL. IN A 
CASE WHERE A CHANGE HAS ARISEN IN A SIGNAL APPEARING AT AN 
EXTERNAL INPUT TERMINAL OF A CIRCUIT ELEMENT, THE ELEMENT BEING 
DESCRIBED IN THE INFORMATION CONCERNING THE AMOUNT OF ELECTRIC 
CURRENT FLOWING THROUGH AN ELEMENT, THE AMOUNT OF CURRENT 
APPEARING AT THE CIRCUIT ELEMENT IS ADDED TO THE INFORMATION 
CONCERNING THE TOTAL AMOUNT OF ELECTRIC CURRENT OBTAINED AT EACH 
TIME. 

1505 STORE, INTO TOTAL CURRENT STORAGE MEANS, TOTAL AMOUNT OF 
CURRENT OBTAINED AT EACH TIME 

END OF OPERATION OF TOTAL CURRENT STORAGE MEANS 

[FIG. 115] 
TIME 

[FIG. 116] 

01 LAYOUT DATA 

02 LPE RULE 

03 LPE PROCESSING 

04 NET LIST 

05 TEST PATTERN 

06 SWITCH-SCALE CIRCUIT SIMULATION 

07 CURRENT WAVEFORM OF EACH TRANSISTOR 

08 MODELING OF CURRENT SOURCE 

09 MODELING OF CUFIFUENT SOURCE ELEMENT 

010 LPE PROCESSING OF POWER LINE 

011 POWER LINE NET LIST 

012 TRANSIENT ANALYSIS SIMULATION 

013 RESULT OF CURRENT WAVEFORM 

014 FFT PROCESSING 

015 RESULT OF EMI ANALYSIS 

016 WIRE -LEAD FRAME IMPEDANCE 
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[FIG. 117A] 
LOGIC VALUE 
TIME 

[FIG. 117B] 
CURE^NT 

h = CURRENT (I) X 2 /OUTPUT SLEW (W) 

OUTPUT SLEW (W) 

TIME 

[FIG. 118] 

1171 PATTERN-OF-SIGNAL-CHANGE STORAGE MEANS 

1172 ELEMENT CURRENT STORAGE MEANS 

1173 OUTPUT SLEW INFORMATION STORAGE MEANS 

1174 CIRCUIT CONNECTION INFORMATION STORAGE MEANS 

1175 CURRENT WAVEFORM CALCULATION MEANS 

1176 CALCULATION RESULT OF CURRENT WAVEFORM 

[FIG. 119] 

FLOWCHART OF CALCULATION OF CURRENT WAVEFORM 

1190 READ OUTPUT SLEW 

1191 PERFORM TRIANGULAR WAVEFORM SHAPING LOOP 

1192 EXTRACT BASE OF TRIANGULAR WAVE FROM OBJECT INSTANCE FROM 
OUTPUT SLEW 

[FIG. 120A] 
LOGIC VALUE 
OUTPUT 
TIME 

[FIG. 120B] 

CURRENT 

TIME 

OUTPUT SLEW 

CURRENT (I) X 2 X RISE -DEPENDENT COEFFICIENT/OUTPUT SLEW (W) 
CURRENT (I) X 2 X FALL-DEPENDENT COEFFICIENT/OUTPUT SLEW (W) 
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FALL-DEPENDENT COEFFICIENT = 2 - RISE -DE PENDENT COEFFICIENT 



[FIG. 121] 

1201 PATTERN-OF-SIGNAL-CHANGE STORAGE MEANS 

1202 ELEMENT CURRENT STORAGE MEANS 

1203 OUTPUT SLEW INFORMATION STORAGE MEANS 

1204 CIRCUIT CONNECTION INFORMATION STORAGE MEANS 

1205 CURRENT WAVEFORM CALCULATION MEANS 

1206 CURRENT WAVEFORM STORAGE MEANS 

1207 RISE -DE PENDENT COEFFICIENT/FALL-DEPENDENT COEFFICIENT 
STORAGE MEANS 

[FIG. 122A] 
LOGIC VALUE, DELAY 
INPUT 
OUTPUT 
TIME 

[FIG. 122B] 

CURRENT 

TIME 

(CURRENT - CURRENT WHICH IS NOT DEPENDENT ON OUTPUT LOAD)X 
2/OUTPUT SLEW 

CURRENT WHICH IS NOT DEPENDENT ON OUTPUT LOAD X 2 / (DELAY - 
OUTPUT SLEW/2) 

(DELAY - OUTPUT SLEW/2 )/ (NUMBER OF STAGES -1) 
OUTPUT SLEW 
TIME 

[FIG. 123] 

SLEW/ CAPACITANCE CALCULATION MEANS 

OUTPUT LOAD CAPACITANCE INFORMATION STORAGE MEANS 
INPUT SLEW INFORMATION STORAGE MEANS 

1221 PATTERN-OF-SIGNAL-CHANGE STORAGE MEANS 

1222 ELEMENT CURRENT STORAGE MEANS 

1223 OUTPUT SLEW INFORMATION STORAGE MEANS 
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1224 CIRCUIT CONNECTION INFORMATION STORAGE MEANS 

1225 CURRENT WAVEFORM CALCULATION MEANS 

1226 CURRENT WAVEFORM STORAGE MEANS 

1227 ELEMENT CURRENT CALCULATION MEANS 

1228 NON-LOAD ELEMENT CURRENT STORAGE MEANS 

1229 NUMBER- OF -STAGES STORAGE MEANS 




Attorney's Docket No. 32796 

COMBINED DECLAR4TION AND POWER OF ATTORNEY 
IN ORIGINAL APPLICATION 

(Sole or Joint - Foreign) 
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